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在 非 平 静 大 气 中 飞行 器 的 运动 方程 

关于 风 对 飞行 器 适 动 的 影响 … aaa 
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有 成 切 变 时 的 和 运动 万 稳 … IER 


飞行 器 运动 的 计算 和 分 析 方法 
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作为 线性 系统 的 飞行 器 的 动态 特性 分 析 *… 


飞行 器 平 衡 状态 的 计算 … 


附录 四 抑 数 在 飞行 器 和 运动 方程 中 的 应 用 
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1. 物理 矢量 用 黑体 字母 表示 。 例 如 
Yvo E 
2. W 
让) 当 写 出 双 部 元 素 时 ， 用 E ”J 表示 。 和 例如 
eosa sing D Uad 
| ~ sing coag 0 | D | 
0 0 1 Ved 
入 ) 当 用 字母 代表 整个 短 阵 时 ， 用 黑体 守 母 表示 。 并 且 ， 列 阵 往 往 用 小 写 宇 母 表示 ， 
例如 r,; 其 它 阵 用 大 写字 母 表示 ， 例 如 Bi 
， it) 笠 对 称 考 阵 ， 用 上 标 ss 表 示 ， 例 如 
0 


tx sr 一 Gs Dy 
| zonnal 加 > 0 Tx | 
O; “Oy Ox 0 


iwp HHE- IRR, 
HERRE EITEK 





ds 


Ki Mate A EE A AKT SS Ara EA, a "Real SIE H. Mm 
EES, SIN GRE) 的 计算 以 及 稳定 性 和 操纵 性 的 分 析 是 十 分 必要 的 ， 击 这 些 工 
作 都 必须 以 改行 器 的 运动 方程 为 出 发 点 。 

不 同类 型 的 飞行 器 ， 由 于 飞行 条 体 的 不 同 ， 描 述 基 运动 的 方程 也 不 相同 。 例如 ， 对 于 研 
帘 常 规 飞机 或 战术 导 强 来 阅 ， 了 可 以 把 大 昌 当 作 平坦 的 、 不 旋转 的 ;而 对 于 研究 超 禹 速 的 飞行 
器 《如 弹 道 式 火箭 或 人 造 卫星 运载 火箭 》 米 说 ， 大 地 的 球 撒 和 施 转 就 成 为 不 可 忽略 的 闪 素 。 
在 麻 洲 大 气 中 飞行 器 航运 动 方程 就 比 在 平 贱 大 气 中 的 运动 方程 复杂 得 多 。 进 行 飞行 器 航 秋 这 
算 时 往往 委 用 全 基 的 运动 方程 ;而 进行 飞行 器 的 稳定 性 和 操纵 性 分 析 以 及 控制 系统 设计 时 ， 
则 常常 条 用 线性 化 的 运动 方程 。 紫 外， 即使 在 同样 的 飞行 条 件 下 ， 运 动 方 程 也 可 以 刻写 成 儿 
种 不 同 的 形式 ， 这 冲 能 是 有 由 于 研究 才 的 传统 习 从 、 或 者 其它 特点 决定 的 。 所 以 在 航空 航天 各 
技 文献 中 出 现 多 种 多 样 的 飞行 器 运动 方程 ， 有 关 的 科技 人 员 应 深入 理解 这 些 方程 的 实质 和 特 
点 ， 帮 根据 实际 情况 如 以 选择 。 

此 外 ， 任 何 文献 〈 包 括 本 蕊 ) 中 列举 的 飞行 器 运动 方程 都 不 可 能 是 包罗 万 稼 的 。 利 技 人 
员 必 须 具 备 根据 具体 的 特殊 情况 建立 运动 方程 的 能 肪 ， 玫 能 解决 名 种 各 样 的 有 类 飞行 蜂 运 动 
的 问题 上 

本 全 就 是 为 了 满 赴 上 述 的 需要 而 写成 的 关于 飞行 器 运动 方程 的 一 裔 专 募 。 一 方面 ， 系 统 
地 欣 民 帮 诊 迹 各 种 飞行 器 在 各 种 基体 条 你 下 的 过 动 方 种 ， 另 -万 耐 ， 讲 述 建立 飞行 器 运动 广 
程 的 -… 般 原理 和 方法 。 

EBRAR mem. ai, ëng HR: ei UAN 
mm DUR RI a lm MAE. EIERE ARIS (IP EEN SE de E 
旋转 ) 建 立 飞行 器 的 运动 方 程 ,这 些 方程 适合 于 迷 度 不 人 高 的 飞机 和 战术 导弹 。 第 ew E 
到 地 球 的 曲率 和 旋转 建立 超 高 速 飞行 器 的 运动 方程 ， 而 且 对 于 特殊 类 型 的 飞行 器 一 一 有 弹道 式 
火 第 或 人 造 卫 品 运 载 火 箭 - 一 -的 运动 方程 作 了 专门 的 描述 。 第 四 章 介 绍 飞 行 器 运动 方程 的 线 
性 化 形式 ， 这 样 的 常 系数 线性 微分 方程 组 适宜 于 飞行 器 的 稳定 性 和 操纵 作品 质 的 分 析 以 及 控 
制 系统 的 设计 。 第 无 章 讲 光 在 非 平静 大 气 中 飞行 器 的 运动 方程 ， 基 这些 方程 可 以 研究 飞行 器 
对 风 切 变 和 大 气 系 流 的 响应 。 第 大 章 介绍 飞行 器 运动 的 计算 和 分 析 方 法 ， 包 括 ， 作为 线性 条 
统 的 飞行 器 的 动态 特性 分 析 , 作为 非 线 狂 系统 的 飞行 器 运动 微分 方程 的 数值 求解 ,《 行 器 平衡 
状态 的 计算 。 

本 书 可 供 航 空 航天 高 等 院 科 的 高 年级 学 生 、 人 研究 生 和 数 师 以 及 航空 航天 部 门 的 科 按 信 员 
作为 系统 学 习 飞行 器 运动 方 径 用 的 数 材 , 亦 可 作为 科技 工作 者 的 实用 参考 资料 。 此 外 .对 十 综 
合 大 学 力学 条 的 师 生 也 有 一 定 的 参考 价值 。 





SS HAAR 


§ 1 坐标 变换 的 一 般 方法 和 性 质 


LI 坐标 变换 抢 阵 的 定义 


设 : Y RIRE r 在 坐标 系 $，( 即 Ooxpypzp)》 DIS, (BIC, ze Neie l 中 的 分 量 分 别 
2 (xpo Yp? Zp) 和 CBR Han za) {地图 1~13》 o 它们 之 al 的 RE CL nt wb 


ZS, 

Xy bu bus Dis Xp 

| Ya L bar bss bas] | De ! (1-1) 
Z4 ba bag Pss Ze ， 

其 中 的 5 , 即 是 对 应 各 轴 之 问 的 方向 余 嘴 ， 


bi rus en ët bya = cos( Ypo Xal bis = cos py xg) 





bat = cos(x p> yea); bas = cos{ Ha: ët bss= pos ze: Del (1-2) 
bau 一 cosf xps Zg)» bess=ecos(yp; Zei bas=cos(2p; Zei 
引用 符号 
Xp Xa 
,| Re | d Ye (1-3) 
Ze Ze 
bu Aus bus 
Bs=| ba baz bas Get 
ban Daa bas 
WA C-D) 可 以 写成 图 1-1 内 标 肥 Sp 和 SS。 
ra= Bir, det 


MEB gg ERAS gl ke 5: EERE. PRERA IEE PLATEN a ds 
变换 矩阵 的 问题 。 

因为 在 飞行 力学 的 六 多 数 坐 标 变换 问题 中 ， 只 注意 坐标 轴 的 指 淘 ， 而 不 关心 坐标 原点 的 
位 置 ， 所 以 以 后 将 主要 研究 具有 共同 原点 、 不 同 指向 的 两 个 坐标 系 之 间 的 关系 。 


1-2 坐标 变换 矩阵 的 基本 性 质 
若 仍 然 利 用 式 {1-2) WHER SARRE., ME 
Ze buy ban ba Za 
[> -| bus be Daa | È 
Ze bis Das bss Ëg 


或 r= Btr; (1-86) 


EI ba ba 


其 中 B= bis Bas Dar (1—7) 
bag brs Dag 
比较 式 (1-4) 十 (1-7) ， 可 见 
Bs=(B3)” 或 日 一 (日 人 t 1-8) 


AR eT RERE. Gsm, B EBIT AH EHE. 
GR, BA (1-5) wm 
rp =(B3)tr, 
HERA (1-6) 比较， 得 到 
Bs=(BY)"! 或 Bi-CBri) WE 
ZERE, B5B FARE. 
(DÄ (1-8) d (1-9) 可 见 


EHNEN (1-10) 
HAIA EREEREER AARRE RR ER CzzmbEpnA ek, 
3 1 当 i 二 时 
Zeche) a ais int 
, f WEI: 
È babash | AISAN | 
CH O remm 


ak, Eir ERRER EiT (或 列 ) 的 二 元 素 的 平方 之 和 等 于 1; 不 同行 【或 列 ) 的 
元 素 彼此 乘积 之 和 等 十 0。 根据 解析 几何 的 定理 ， 也 能 二 楼 推出 上 述 结 论 。 
Aa (1-8), (1-9), (1-10) i (1-11) 反 轴 了 坐标 变换 紫 阵 的 重 机 性 质 。 


1-5 坐标 系 的 基 元 变换 

为 了 骨 完 三 维 坐标 系 变 换 的 -… 般 关 系 式 ， 首 先 来 苗 虑 三 维 坐 标 系 变换 的 最 简单 情况 《或 
者 基 元 情况 ) ， 即 有 一 个 轴 征 此 重合 的 情况 。 

设 ， 有 两 个 坐标 系 .和 :SS$。， 其 中 的 轴 z 和 > 是 重合 的 ， 而 xp ypx y HRZ 
fjas RERA S t tt fa WEARS ,的 位 置 ( 见 图 1-2a》。 在 此 情况 下 ， 


显然 有 
Xa cosg siog Ô Xp 
Ya D -sing eoa Ô [> {1-12} 
. Ze 0 0 1 Ze 
或 者 写成 
r=B.,(a)r, 
cosg sing Q 
其 中 dE -sing cos Ü | {1-13) 


0 0 1 


是 在 上 述 条 件 下 由 上 坐标 对 SS, 玫 $ 的 变换 矩阵 。 
同样 地 ， 如 果 坐 标 系 S$* 和 ,S ,的 轴 y5 与 yo 重合 ,并且 当 5, 绕 轴 y, 转 过 角 a 时 就 达到 5。 的 
wR ( 见 图 1-25》, 旭 由 S ,到 5 。 的 变换 卸 阵 是 


A e 


vos 0 -sing 





dE 0 1 0 | (1-14) 
sing D cos 
ikii. WRA Si, 03.00 hp w BI LS E E, El Ha ol A, DÉI 
位 置 ( 见 图 1-2c》， 则 由 SS; 到 .5 pn 
Py - 
ër 
ka 
Ze | Ee 
DN AN 
ta) iÈ) Lei 
iz Seefe 
1 0 0 
d 0 cosd sind (1-15) 
D -sind rosg 


式 (1-13), (1-14) 和 (1-15) 所 表 
TZAT IRER, iO E EE. 


1-4 坐标 系 变换 的 一 般 情 况 

任意 两 个 空间 坐标 系 之 间 的 关系 可 以 
用 -2H (ZF) Euler 角 来 表示 。 

RE: HERRAS AS, HEREA 
过 如 下 三 式 转动 来 实 现 的 〈《 见 图 1-3) ， 
Ok, ERRO, ype SE, EE 8 
i GEIR) ， 成 为 Ox'y'z,; ERR Za 
VREA ESA) ， EA Oxyz: 
最 后 绕 轴 x。 转 过 角 三 〈 自 转角 ) ， 到 达 
Dxoyoyeo 角 上 .7 尖 组 成 一 组 Euler Ho 

根据 ia-3 节 所 讲 泪 的 原理 ， 可 以 写 出 


x’ Xp Xa x’ Xa Xa 
eei Geo) EJee 
pOg ze" Ze z’ 


Zb Z 


把 以 上 三 式 合 起 来 ， 在 





图 1-3 eum EES 


ro=B.(8)B,(n)B., Er, 
Pris, ale, 的 坐标 变换 矩阵 为 


B8=8-(s)B,(2)7B:(C) (1 .16) 
HAB. CS), B O), BCE) 的 表述 式 可 以 按 公式 (1-13), (ti, GI SI. H 
Eai, p, Za, EB, d U-16) RRAZ TAM BSAS 对 应 的 旋转 租 序 
是 相克 的 。 

IE HEY TARH ds br kën mu, ptaz, dek, 经 过 怎样 的 ASS 
AAR S, MER: AIEN S, EE 
BFERZ AGERRE, LAE REREN OA EAA E 元 旋转 中 
的 侠 座 l 

DHR EAI ERA A E RE mm, eiA e AMARRA ` 
组 ， 抽 在 天 给， 这 是 因为 ， 从 一 个 坐标 又 到 另 ARAE, ARA T Ah A jafi MA 
序 进 行 : 

Xo $F TR, Wi 
的 一 组 (二 个 ) Euler DER. Häipe e E RERE EA ët rel, Di te E 
ali ZS gie Dn, ang (SPA 《1-1)、(1-2)) 。 

企 实际 问题 中 ， 遇 个 些 标 采 之 间 的 相互 方位 ， 往 往 AR 用 确定 mm ZB tele DEN, 
这 样 的 入 择 常 常 是 由 传统 习惯 决定 的 ,或 者 由 物理 意义 决定 的 ， 或 者 由 测量 方式 决定 的 。 例 
如 总 机 机 体 相 对 于 地 而 的 方位 ， 用 篇 怠 角 多、 入 仰角 # 和 滨 竺 角 ? 来 确定 〈 华 见 第 一 章 82) ， 
mA EG ARE M 


1-5 “坐标 变换 矩阵 的 传递 
让 用 -个 全 标 系 $  ，Sy 和 S。， 且 已 知 变 换 年 阵 B8 和 Bg， 如 何 求 出 B8 


因为 rz Biren, re Bir, 
E r,=B;Bšr, 
mH ie ii 
Steel (1-17) 


这 就 是 洗 标 变换 的 传递 公式 。 
1-6 ”物理 矩阵 的 转换 


彼得 ， 有 两 个 临 理 名 是 上 和 YY ， 它 们 在 任意 至 标 系 中 的 分 量 由 如 下 关系 式 联系 : 
YA (1-18) 
och 13-20, EJER N E T ERRAR R o 
在 坐标 村 SS p Po uA y 的 分 量 的 关系 是 


die (1-19) 
在 夺标 系 S i 
ve A (1-20) 
WA. TEUAN EE PREA KIE BrEA KEER RR IDERE: 
ec Bir, u= Bgu (1-21) 


fRA AC (1-19) t: 


= p 
Bv, =å „Btu, 








d o, 利用 式 《1-9) ,得 到 

Yo=~BrA, B+u, 
把 它 与 式 (1-20) 此 较 ， 就 得 到 

A Bi. Br (1-22) 

ETARE hetk S, PRERE A, RRHH pA S, nä, , KhA 

3x 34ER, HARAR (1-18) 所 表示 的 物理 意义 。 例 如 ， 后 面 将 看 到 ,刚体 的 动 直 利和 天 
EH 与 角速度 矢量 巴 是 通过 惯性 年 阵 1 seg (ELA (1-40) ) ， 因 此 可 以 按照 公式 
(1-22) 把 霸 性 矩阵 (及 共和 拇 一 个 元 素 ) 由 一 个 举 标 承 的 算 到 另 一 个 坐标 采 。 利 用 公式 (1-22) 
DS -个 便 子 是 气动 旋转 导数 的 换算 【 见 第 四 章 86) o 


$ 2 旋转 坐标 系 中 矢量 的 导数 
2-1 基本 情况 
设 ， Aak SCOcuz RP oC HRIAN., REO- R HNG 量 为 ws, 


tys vi, HI 
ozstdsl two jt oK 
其 中 i 、j 、k 大 x、y，z 轴 的 单位 无 最。 
设 有 一 个 可 变 的 和 关 重 af) ， 它 在 此 坐标 对 中 的 分 量 
Ha. Gys Qg: H 








a=a,i+a,j+a,k (1-23) 
HEARRE a mt 上 的 导数 : 
dä do, . da . 
dt Ddi it dc i+ 
一 E k++ a, d +a; zj +a: ak (1-24) 
理论 力学 已 指出 ， 当 一 个 刚体 绕 定 点 口 记 角 速 座 包 旋转 有 时， 刚体 上 在 何 一 点 也 的 束 府 
SI a SET (1-25) 


Ah r MORBA PRIR IRo 
现在 ， 把 单位 舌 景 i AERE eh Sim, 于是， 与 式 (1-25) 类 似 。 
可 以 号 出 
di : 
dr TOi 
HEER, H 
dj ， k 
est, Zack 
把 以 上 关系 式 代入 式 (1-24) 中 ， 得 到 


da de, de e de 


= de togi It gr Kt Oxorita,jta,k) 





了 最 后 得 到 旋转 坐标 系 中 估量 的 娃 数 的 公式 


da _ òa . (1-26) 

di Tg tO*a 
` , da — do, ， doy ， do: | 27) 
其 山 A, =- yi- Ti j= di k {1 27) 


5a707 称 为 在 活动 坐标 系 中 的 “ 租 对 学 数 ”， 相 当 二 站 在 此 活动 坐标 系 中 的 观 娩 者 所 看 到 的 
关 最 a 的 变化 率 。 而 dayvdt 则 称 为 “绝对 导数 ”， 相 当 于 站 在 固定 坐标 系 中 的 观 周 者 所 看 到 
的 矢 最 a 的 变化 奉 。 例 如 ， 若 aa 是 某 点 的 秋 径 ， 则 3ay767 代 表 该 点 的 相对 RE GA r DE 
标 系 ) ， 而 daydt 则 代表 该 点 的 绝对 速度 。 


2-2 更 一 般 情况 
现在 把 公式 《1-26) al ELE. DER Steil A DEI, ie 
AS (Oaz) TE EECHER EECH Ao: Æ RR aU) 可 以 分 别 在 
EZE AS :中 分 解 成 
AEN E taaki =S&asigtayzja taske 


根据 关 晤 导数 的 公式 (1-26) 和 (1-27) : A 


GEREENT 











Za düxi ， dayı E dazi 
EE di Vt ar jt gr ki 
其 中 
ĉa \ dän, days ; daz 
Lol: di eT gp Er k: 
FAHRE a Ep S S a PEE o 
由 此 得 到 
óa ña ， 
(ar) =h 8 ),+ (0:- oxa (1728) 


2 dch kR & ETARA AHRR. N o- 0) SS 


§ 3 Rss 
大 家 知道 ， 矢 量 & 和 的 及 莘 积 a xb 是 一 个 笑 基 ， 丰 坐标 际 x8z 中 它 的 分 最 等 下 
apb aby 
a Ds axb: 
Gry Ost, 
ane kan, TASE 


D -az Ge KÉ 
el ER 0 -axs | by | (1—29) 


Tay Ux D 


EREHE MEME a A EPE (sherw-symmeiric mairix) 。 为 了 今后 局 
FPE., RAin THER S 


二 是 公式 (1-29) WHEE 


aa yes y Bx 
xl a, | | b, (1-297) 
ER D. 
dr axb=a"*b {1-29"} 


另外 ， 点 乘积 &，5 是 一 个 标量 ， 等 干 


8 axb +a, by tab 
Aur Seit, SW 


b, 
ab 一 [os。 a, aa] | by | (1-30) 
b, 
或 ah 一 a7b (1-307) 
公式 (1-29) 和 (1-30) 就 症 矢 最 的 乘积 的 矩阵 表达 式 。 


§ 4 山体 的 一 般 运 动 方程 
4-1 在 活动 泽 标 系 中 列 写 的 刚 体 质心 动力 学 方程 
在 亩 论 力 学 中 知道 ， 刚 体质 心 的 运动 方程 为 


dy 

ITT a 

Ehm ARIRE: ARER OHTI AREER; F 为 作 肯 在 刚体 上 的 外 
JERE. 

设 ， 有 一 个 活动 如 参 才 坐标 系 O 〇 xyz， 它 的 砍 转 衣 速 度 为 mm。 在 该 坐 标 At. oko E 


KÉN Gys (en 即 


-=F (1-31) 


OSv,i tæ, jtk 【1-32) 
其 中 i 、j、kk 为 党 x、y、z 辅 的 单 仁 人 疾苦 。 
同样 ， 质 心 速度 v 和 力 上 也 在 该 活动 坐标 系 中 分 解 如 地: 








v=v itv ]jtv,k 
Fast | (138) 
根据 在 旋转 坐 妹 系 中 估量 的 导数 的 公式 〈1-26) ， 方 程 【1-31)》 可 以 写成 
å 
nl g + EA) )=F (1-34) 
6 do, duy ， ad- 
dä IS Er EN 





利用 公式 《1I-23289) ， 则 方程 【1-34) 可 以 写成 如 下 于 和 阵 形式 : 


d Uy D =o, Oy Ue Fy 
A7 | 用 | " ? 97 | J- ru | (1236) 
Ua 一 I Fa 


ER EP d Ve 


ROO BiA aA T : 


m(- ER "Foot, -oo )= E, 

duy (ER 

Wo i 一 
m( T; t oUr Eet, 1-37) 
(F v )=F 
m Fr F Uy Oy 


这 就 是 在 任意 活动 做 着 坐标 条 中 列 写 的 刚体 质心 动力 学 方程 。 在 以 后 的 两 章 中 将 把 方程 
(1-37) 应 用 到 不 同 的 坐标 系 。 





4-2 刚体 统 质 心 转动 的 动力 学 方程 


理 诊 力 学 中 得 知 ， 刚 体 甸 动 的 动力 学 方程 是 


dH 
-5 M (1-38) 


KM EHN ER EAD AREG DEAA BR. 
用 是 刚体 对 质心 的 动量 给 ， 其 一 般 表 达 式 是 
H=frx (ow xr)am {1-39) 
其 中 “四 是 刚体 的 转动 角 速 诬 ; 
r 是 由 质心 O 〇 到 质点 dm 的 矢 径 。 
若 以 质心 为 原点 ， Si RIEDER AE SO, Mai RH H 与 角速度 o 


的 分 量 有 如 下 美 承 : 
H, I; -Ia -dyr DI e 
H Se | >, | (140) 
FH, —- Tox” S d, D e 
其 中 T= f(y +2?)dm> 1,= har xdm, 1 [exe ty ite, 
I= Led I p= yedm, I 一 12x 人 1 


PAPERE ale Ee Ou: HP ale Rn, mAA (1-40) 布 迪 的 方形 阵 称 为 惯性 
阵 。 

把 庭 转 坐标 系 中 疾 量 的 导数 的 公式 (1-26》 应 用 到 刚体 旋转 动力 学 方程 (1-38) 中 ， 有 

OH ` 

Ât 

利 骨 公式 (1~29’) ， 把 上 式 写 成 组 阵 形式 


D 10 a 


+o xH=M (l-41) 





H, One, fa M, 
alee (a) | 
H, o. H. Ma 
JFJ 2 
He io, H,- o, H, M: 
dH 


Tt eHe- oHa M, (1-42) 


一 一 -一 -一 一 -一 


PH H p H Dk 101 A ER, ART, JN, Kai IEN 
FMI SSC A gp eh D Le bn SE AC än E, 


do, z 
Si a -y -I)a ploi- oi- 1,.( rn )- 


1 H 





-fay Del 一 中 Da E 


do, 


1 EE EE 





(1-43) 
-fye (2 -= g Úy jam, 


do, 


ege (i, -F OOY (tee 0 -7,.( di tOO }- 


一 了 (22 Oyz E 


SEAN, KIA h N E E E Pe RHE audio aëérgIRsCTIal E 
行 器 的 转动 运动 产生 明显 的 影响 ， 它 称 为 陀螺 效应 。 因 此 在 飞行 器 运动 方程 中 有 必要 计 政 这 
作为 最 一 盘 情 况 ， 假 设 旋 较 部 件 的 动量 和 矩 在 与 飞 行 器 固 连 的 坐标 系 中 县 有 二 个 分 量 ， 
WITZ ETS OT ERC {1-44} 
TRT KITARAA BIRER (1-40) PEDS mH- as Hrga H. 
ELR A AE (1-42) 中 ,于 是 飞行 器 转动 动力 学 方程 应 懂 正 成 为 





1 
H 








GC T Odi ex = Drd = M, ( 1-45) 
meg Hra t+ midry— e — | 


其 中 第 省 略 的 部 分 C0 MEFE (1-743) 的 左边 。 


a II ew 





BOE "PKW arm Kirech LE 


Aë ERHET Rik, MERZER {地球 】 的 曲率 和 旋转 的 条 件 十 建 焉 网 
性 飞行 疾 的 运动 方程 碟 。 对 丁 日 前 去 多数 类 型 的 飞行 器 ， 即 备 种 飞机 和 战术 哇 弹 求 说 ， 把 天 
地 春 作 菠 面 是 适 宜 的 。 所 以 本 章 建 立 的 运动 方程 有 广 阅 的 应 用 场所 。 考 虚 大 地 出 宗 和 放 转 冉 
的 芒 行 器 运动 方程 将 在 第 沁 音 进行 妍 究 。 

在 将 中 在 定义 从 标 系 和 运动 变量 时 劣 虑 风 的 存在 ， 世 在 建立 运动 方程 时 假 普 大 气 是 平 衣 
A (MRETI) 。 在 非 在 散 大 人 气 中 飞行 医 的 运动 方程 将 在 第 五 章 中 讲述 。 


$ 1 坐标 系 的 定义 


1-1 WA 

DUPEK EHIE AREE TERA (或 你 pr ph 
R) o MIURA S PEAR. DA -一 指出 。 

EREE A Aide KITA DULAS PA AROE E A hb E, HRE. MT 
v Sr, KIRA EOI TEER (an jeet) RERA AEE BRAE 
AR u WRR. UUR ANRE RAM E Sw Ko EMEBEREK 

YY 二 出 十 W 

4R ARESA Didj, KITA ASHER -G Ae ëm, SR, EER 

KIPRE, 


1-2 ”地面 坐 标 系 


地 而 CRH) 侍 标 系 Cpxpypzp【 见 图 ?1 或 2-2) 国定 于 大 地 。 它 的 原点 Qo 选取 人 在 地 
面 二 的 基点 ( 俩 如 起 飞 点 )。 轩 yp 错 恒 向 上 。 办 xpn 和 zp 在 水 乎 而 内 ,它们 的 指向 可 以 接 
上 县 体 条 件 规 以 选取 (没有 严格 限制 )。 
河 为 不 考虑 大 地 的 旋转 ， 所 尺 此 坐标 系 可 以 看 作 是 惯性 坐标 系 。 飞 行 散 的 位 置 和 将 夺 以 
Senioren, 
亚 把 藉 点 放 在 飞行 器 质心 mietit Eir THERA f xo. Hu, zos Wax HE 
构成 的 坐标 冬 Oxoyszs 称 为 飞行 器 奉 连 的 地 面 坐 标 系 〈 胸 图 21 或 2-2) 。 当 仅 关心 坐标 柔 的 
专 位 ， 而 不 甘心 原点 位 置 时 ， 此 举 标 条 可 以 简称 为 地 面 坐 标 系 。 





1-3 MEERA 


Hikäbkbeëcse ez, 〈 见 图 2-1 或 3-3)》 固 连 于 飞行 器 。 原 点 口 在 飞行 器 质心 。 纵 铀 x: 乎 
行 于 机 身 炙 线 或 平行 于 折 井 平均 气动 奖 ， 指 办 前 方 。 又 轴 zs CSS aire, weg"! 
Hx AEA. Wehe 重 雪 于 飞行 器 对 称 平 面 ， 指 向 布 方 。 


IER 





CEA EELER 


Ap 





48 KEM Xp 
i 


et Tat Ai] 
GEREENT E E EE 


1-4 气流 坐标 系 


RWBER Ox Yeza 〔 见 网 2-3) 的 原点 口 在 飞行 器 质心 。 辅 xv* 沿 飞行 器 窑 速 矢量 u 。 
äu, fc rage, era, 、 指 向 上 方 。 轴 >, BETEK ， 指 向 右 方 。 

通常 作用 在 飞行 器 上 的 宗 气 动力 是 在 此 华人 奈 对 内 分 解 成 阻力 Q@ 《党 x。 负 向 ) 、 升 力 了 和 
WD ZPJ 


1-5 MERR 


Bäbb ROX zi (HRZ 31 的 原点 O 在 飞行 RELO. fh x ii KIT F E REY. 
轴 y; 在 包含 矢量 * HEREHA, wA Fx WELF. Me BE TE 重 平 面 xy (D 
市 是 水 平 的 ) ， 指 向 右 方 。 , 

建立 飞行 器 质心 运动 方程 时 ， 要 利用 此 泽 标 系 。 

SV, Edi Së, 与 气 流 坐 ten, SH (oni 时 是 不 … 致 
Di. HEA (w=0) DF. 


e 13 > 





On 
2 水 平面 xa 
1-6 半 机 体 坐 标 系 





DIE be Oasen 
( 匈 岁 2-37 的 原点 吕 在 飞行 器 
BC. Sé, WARRE Uu dr Ar 
Zens: tr. hys 也 在 对 
EPEHA; EFA rix 指向 上 
Zei IN Lëäu, Si. 。 轴 zx EE 
于 飞行 路 村 称 平面 ， 指 阿布 方 
(Rm bü: Si, 

有 时， 空气 动力 的 分 最大 在 
EHLI br A RR a 


H 
图 2-3 Ges, 、 牢 机 体 坐 标 系 和 机 体 坐 标 和 的 关 楷 


8 2 各 坐标 系 之 几 的 关系 及 变换 年 障 
PERRE ÈK RAEI e Bn Euer D 来 确定 ， 而 根 据 这 些 Euer 角 的 定义 就 可 以 
导出 坐标 条 之 间 的 变换 给 阵 。 
2-1 地 画 坐 标 系 与 机 体 坐 标 系 的 关系 
这 节 个 华 标 和 对 之 问 的 英 邓 用 以 下 三 个 角 认 来 确定 【 见 图 2-1》。 
e jd e 











at md, Mk, ERKE e KITEE op DD DE. Au, D EM 
Pë, Wehr, 

MA: Wiki SKEM azad EIDEN. x A EARE lr, 

滨 转 角 CERD v: VITRA EE y 与 包含 轴 x, MEER A E ie A 
PAAR IGATA, A8 E, "DE, 

这 三 个 角 y、 相 ,决定 了 飞行 器 在 空间 的 妆 态 。 

地 面 坐 标 条 Sa 可 以 佐 次 经 过 三 奖 转 动 ， 达到 瑟 机 体 坐 标 际 5S, 重合 (图 2-23) ek 
Xaysg24 绕 轴 y4 转 过 和 ， 成 为 x yarr DEER Eh éd ESO, 成 为 x:y 2 iehi efh x SE 
过 ?人鱼 ， 达 到 与 x,yizi 重 合 。 注 意 ! x*、y’、2' 厦 中 则 过 滤 轴 ,它们 并 不 构成 正 交 轴承 。 

DR, EA gih Sea, 可 上 太 写 出 由 ;到 S ,的 变换 什 阵 为 

Bi=B,(y)B, GHB, (4) {2-1} 
RRA Q-13), (1714), (1-15) HHH EERE 
了 0 0 
dE 0 COM siny | 
D —siny cosy 
cost sinf O 
dE -aint cost 0 
0 0 ] 
co% D -siny 
dm 0 1 0 | 
siny D cosi 
RAR (2-1) p, 8HRS a HOUE AS H nE A 
cos cosy : sing : — cosh sind | 
-vosy sint eos ? 


cosp rong | 
+ siuny sin : i +sinp cosy : (ESCH 
sa a 


cosy sind sivy | 





B; = 


. I ~ siny sinf sin ` 
- sinp cost H v F | 


teosy ang : ; 十 cosP cosy 


1 Weit, 
| sin sinf vos g : 


2-2 地面 坐 标 系 与 航 迹 举 标 系 的 关系 


EG E KEE SEET EE EG ECH 

friem (mn) pe MaRS (Dim Bai 在 小 平面 xz 上 的 提 影 线 
x“ 与 地 面 系 的 轴 xs 之 问 的 角度 。 当 速度 v Di LIRE, DE, 

Stand, fies OERE Y) 与 水 平面 xszz 之 间 的 角度。 当 速 E v IL 
CRI KITARE HO 时 ，8 为 | 上。 

角 闪 :和 8 决定 了 飞行 器 地 速 失 量 ¥ 在 空间 的 六 向 。 

地 面 谷 慰 系 Ss 可 以 通过 如 下 的 两 次 转动 达到 与 航 迹 坐标 对 5S ;重合 。 首 先 xyyazs 线 轴 ya 
Hf MAr yaz RMA hiza AO, AAH 

因此 ， 由 Ss 到 3 AE RERE 可 以 守成 








l B;= B.C0)B, (Cy.) . ZAE 
RAR (1-13) A (1-14) 与 出 B (OMB, (PARA (2-3) 中 ， Rit, PAHA RE 
br 25 200 e e A R E RAPE Au F a 
vos une, ` inf —eolsing, 
B+:= sinf eos, eos sinf siny, | HERR 
Sing, 0 Coag, 


2-3 Gäre AL ës SD e 


12 tb, R BÉ 5 Him pn A ETE Wl). 
Dë më, "me In, (ERK Rui D'Kkiräë P y A num, DR 
Ria Kiräm Im (nm k O WOH, 
pa Ti E pu, Mie ik Ru) Faer E Lon, DRUI, HI mm 
BR. FAGE u ing PIURE E A x: e ti AA F Woh, 
节气 流 坐标 系 Ox, Saz, MRR y IAA WD Out (也 就 是 小 机 Khi 
Oxz) o Rietz tte fH, ep WLP ERO Heo 
DE, RA iiair AEHL En ARAI RHR D 
BB toi (£) (2-5) 
HER EFI 
cosdeos f sina = positsin f? 
Bi 一 | -singes fl eose  singsinf? | (2-6) 


sinp 9 vun 


2-4 SBS bake SEL rs DOC 

Im LR, Take SSC EES ab be gp KE ko 5 AER CH 17 
问 的 电 何 其 系 已 经 介绍 了 -。 因 此， 很 容易 得 到 如 下 的 丙 个 变换 扎 竹 。 

tiA EAE pr A PP HL iE or E hy TE e ERE 


wap 0 -sinf 
au) 0 1 0 (2-7) 
sinf 0 coa 8 
出 年 机 体 坐 标 系 到 机 体 坐 标 系 的 变换 年 阵 
cosa sing Ü 
vue) — sind cosg 0 ! {2-8) 


0 0 1 


2-5 航 迹 坐标 么 与 气流 坐标 标的 关系 


在 IEN. DDR (Ww 志 0)》 的 情况 下 ， 这 此 个 坐 际 系 的 关系 非常 复 杂 ， 因 为 它们 之 
ben vg, d o, o, HERE {以 后 通过 冉 2-5 将 会 看 得 清楚 〉》。 如 果 要 得 到 变换 
SS, 281 


D 16 H 








B/{=B7B;B7 (2-9) 


kp, HATH. 

TERRIER, MEA {w=0) 的 情况 下 ， 它们 之 间 的 其 系 十 分 简单 。 畴 为 这 Wär 

Jj u ms, A pea Dh, B e a AARAA DR, DEI SE Z Din, 与 

BE: Sp ën, 2 Dir AER, bs, A Bim fär RRE Bea DEI POE E K E A) Ya o 
Mfz FEIS, ehr (W24) 。 

EEIE TO An ëtät ab AS tr ROX Yrer Zë a 转 
过 角 ?>-。， 训 与 气流 坐标 系 Ox, 重合 。 因 此 变换 第 
EZ 

1 D D 

B}=B.(v,)=| 0 ya siaa | (2-10) 

D = sinye Oyy 


(w= 0) 


2-6 ” 池 面 坐标 系 与 气流 坐标 系 的 关系 


.在 有 风 的 情况 下 ， 这 贞 个 人 准 标 条 的 关 邓 只 能 通过 机 
eink WEAR (1-17) ， 有 








图 2-4 sites 
Bs= GIE (2-11) 姜 举 标 系 的 美 系 


共 中 包 合 据 个 角 ， y y 下、 o, Bo RIERA REEE EI A o 
在 无 风 的 情况 F, TOET ETEN EAA EPR AREE. Di 


Di B4B4 {2-12} 
PA 2-10 和 (20RA, 88 


í Ou (eos: : sin fl : -cos 月 nps D 
~ OS, D rasg, H os? cosd ` tos Yy inĝ Aug. | 
Bs—| siny, sing, : : Tax: of | {2-13) 


和 H 

` iny D > liny, an Aug, 
:sinys Oak ` 

Jugen, sing, : E peost, cospe b 


(us — 0) 
其 二 只 包含 3 个 角度 : d dé, VARAR Sé. JEITZ. 


| - 
| sinp, sin cosp, 


2-7 各 个 坐标 系 的 相互 关系 的 汇总 


以 上 介绍 的 五 套 和 坐标 对 之 疗 的 相 志 关系 吓 以 用 一 -个 示意 图 ( 见 图 2-5) 清楚 而 形 轨 地 表 
PHR. HERTE E Ame (ORRA = 转 过 4 角 。 寄 和 白 小 方块 器 示 过渡 性 的 无 
名 党 标 系 。 粤 线 卉 示 仅 在 无 风 条 体 下 才 存 在 的 联系 ， 实 线 则 表示 在 任何 条 件 〈 有 风 或 无 风 》 
FENKE. F ARRP. 







Rr yl 


i 





m-i 部 华 标 系 和 相互 关 杂 示 意图 


§ 3 作用 在 飞行 器 上 的 力 和 力 算 
在 给 中 飞行 时 作用 在 飞行 器 上 的 力 包括 : 发 动机 推力 P 、 空 气动 力 R 和 R ip, "em 
明 这 些 尹 的 特性 以 及 它们 对 质心 的 力矩 。 
人 于 飞行 器 沿 地 丽 、 水 面 滑 跑 或 沿 发 射 集运 动 时 受到 的 约束 力 ， 则 不 予 讨论 。 
3-1 发 动机 推力 P 


在 正常 情况 《 非 故 障 情况 ) FF 发 动机 的 总 扒 力 在 飞行 器 对 称 平面 Ox 内 。 设 推力 作用 
线 与 机 体 织 轴 x, 怕 成 角 g,， 夺 对 十 质心 OO 有 偏心 焉 e,。# ,和 e。 的 正 向 在 图 2-6 中 作 了 规定 。 
因此 ， 发 动机 扒 力 在 机 体 坐 标 系 中 的 分 虹 为 : 


Pai 5 eosyps Daer Dame: IN (2-14) 
FED D HEAL e iP a Rb, 
Mp=0, My,s=0, Aur Pep (2-15) 


HR. EARRA T HEDA 
EJES RAM tH. KK 
ERM RH RERIT HH e 

一 般 地 说 ， 发 动机 推力 忆 随 
痊 习 行商 诬 4、 窑 速 4 和 油门 于 
度 s5r 而 变化 。 写 成 

P=p(u, ks 6) (2-18) Wor San 

ARARAU AHLERS GER SEI. MERIETE (E e SR, 





3-2 空气 动力 RR 


通常 把 空气 动力 R 在 气流 维 标 邓 中 分 解 成 三 个 分 量 ( 见 图 2-3) : AUO ( 沿 轴 x i 
Di, AAY ( 沿 轩 yo》， 侧 力 Z《 洛 轴 z,》。 因 而 


D 18 D 


D. -Q 
| Ryg -| Y | (2-17) 
Fag Z 

实 气 动力 分 量 可 以 者 示 成 如 下 形式 ; 


n l n2 | 
Q =e, z Pu S | 
Y =e, pw 

Scy a PH S {2-18} 


Z =C S 





pu’ S 


Ahs KERSANA TAEL KTR SAA aA R, ATERA KITA 
(PTR EA TA) , SRAI RER: 4 D KRHA; oe Loire 
TE: cv，cv、c> 分 别称 为 阻力 系数 、 升 力 系 数 和 侧 力 系数 ， 它 们 都 取决 于 许多 因素 。 对 于 一 
定 的 飞行 器 外 形 来 说 ， 通 常 有 如 下 图 数 关 夭 : 
cc, Res @s ò; Ao d. 0,) 
cy 一 er 人 (ae ës w) | 《2-19) 
e =M, Bs un ous EI 
以 了 上 式 子 中 ， 分 号 GO Dionnknaä ez, AE AARAA o 
M zujati Mach asal h Eik SE LES gi A E Recul t ERA Beynoids#, 
RP PLAKTREZKE, =r EBE LKAA AERA. Ae, tT 给 定 的 飞 
AER, A 
M=M (u, h) \ (2-20) 
Re= Relus, A) | 
d. A, ERAMA, JAMRA e. 
在 某 些 资料 《例如 [7 ]) PRE, MD, PORMA SERAH R A PLE AP 
分 显 ， 因 而 列 出 的 运动 方 稳 与 大 多 数 书籍 《 包 斤 本 书 在 内 ) 有 所 差别 ， 达 是 必须 注意 的 。 定 
Dk Stonn D. Da Du, Ya, ZP om. "Dm WD, 


H. Z, 5i TEM: 
-Q -Qe 
yje e) 
FA Z, 
Q =Q0,cos ð - Z sinf 1 
Y=Y, (2-21) 
Z=Qsin f + Zol 
ZS a ARKIT ADO Di AA EEEn APARRA A8, R D EM, 
thx) . RJ EM, Ghau AMH HEM, Gpe) o 
SADT BARRER. AT KEILAA KIR 5 


ev IO >» 











ua, 1 vue 





2 1 
| 
1 | 
GR SÉNG -pu Si | (2-22) 
, i ; 
VM. cm, 2 an Ssh, 


Hop AHARI, BAIRE Cah; me mas ma AEE ARA. MACRI A AR 
o H EREHE -E. LAGERE- ER, A 
m= m (ÄM, Ps Ôr Âo Ors Oyo a} 
My SHu Ma Po Ås Örs war wer a, P) | (2 23) 
Ma =TM M> ds Ô, ws OI 
ZEPA AERAN KITA f +” Enn WIER, AMGA A TA) 
TAHE 人 ,和 冯 。 胡 达 式 中 的 参考 长 度 统 - Ei, 
式 《2-19) 和 (2-23) 式 所 代表 的 英 系 式 通 党 由 风 润 试验 或 理论 计算 得 到 ， 以 气动 曲线 
或 由 格 的 形式 给 出 , 灶 于 除 冲 数 和 je 数 以 外 的 其 它 套 数 , 在 一 定 范 团 肉 村 以 表示 成 线性 丽 数 
的 形式 ， 例 如 





Om, Əm 5 dm, Jm: a 
Ja að, zt w; Cla gå 


BEAG m a FU Ae A MA SE A AAE M B M GEER TE B Bo 

这 里 说 明 - - TRAMRAHEL OLE- o ò ARM fo RIE SARE 
W HURT AYE SÉSaiäNRImzA, SMIE. AAR mf 的 绝 对 值 
MPH. ARMARE Dirk he, CA, rr Dipen, , A 带 也 引起 小 最 
Gët Den, Ò AKR M EBRA, 3, AE. A i AEE EE H 
的 是 产后 俩 航 力矩 ， 但 也 附带 产生 小 量 的 流转 力矩 。 对 于 正常 式 飞行 器 〈 即 方向 能 在 重心 之 


OUT Ee a] 





图 2-7 BA M AER 


a 20 ba 


Fei, EMS, Zë, 而 对 于 鸭 式 飞行 器 《 即 方向 能 在 重心 之 和 前》 > E S, UE 
Din, 656- 站 升降 舵 俩 转角。 当 它 的 后 稼 向 下 仿 时 ，3 ,为 正 。 升 隆 舵 偏转 的 主要 昌 的 是 产生 
ii. ITERA, CAD, AIR, THO, EDER m. 


5-3 重力 6 
重力 的 数 个 为 
LG mn (2-24) 
ite KI RDR: AEH mE. eipmimiiänemat, i Aigh 的 负 向 。 
ut. SHEEJA me 
Gaa 0 


| EE -| -mg t 2—25) 
{7 ,1 D 


wo EJH ET EE ek A 
HARER TAR, JARE H mAN g SER. SAER MAE g Eë 
NEIE 


$ 4 飞行 器 质心 移动 的 动力 学 方程 


飞行 底 质 心 攀 动 的 动力 学 方程 的 一 般 形 式 是 式 〔1-37)》。 在 其 中 取 趟 司 的 活动 从 本 系 ， 
GARE CECR 


4-1 在 航 迹 坐标 乘 中 的 质心 动力 学 方程 
现在 皮 航 迹 坐 标 系 作 为 方程 (1-37) 中 的 活动 华 标 系 。 在 此 坐标 系 中 飞行 速度 的 分 量 为 


Urh v 
[> -| 0 (2-26) 
Vaa 0 
lz: NM. PEREAS papa, Pn, zb, Or d. Steg, 的 
e, D, e, ES bhea: 
prva 0 0 d. md 
| reya L a| d. J- 0 -| gh, cs ) 12-27) 
TEE 0 ĝ Å 


把 式 (2-26) 和 (2-27) RAH- RAR (1-37) ， 且 考 不 作用 力 F=P+R+G 的 分 
EREA (2-14) 、 (2-17) 和 (2-25) ， 得 出 





d. inf sap H 
A d | 0 | $. cns ) IW ])- 
0 0 
Prosy'p — Q d - 
GH ze, d Y Jsf -ma (2-28) 
0 Z 0 


把 BY、B {无 风 情况 ) 和 B33 的 表达 式 (2- 人 )、{2~10) 和 (2-41) RA bÆ MF 
后， 就 得 到 在 五 风情 况 于 在 航 迹 坐标 系 中 建立 的 飞行 跟 质 心动 力学 方程 如 下 : 


Zeie gelos fi -Q - masing 


du - 
DI y eola + gesin f siny, + sinf a+ KI zi 
di (2-29) 
+ Yrosys - Zsiny, ~- mvosd 


dp, ， - 
- murot d -= P| — coat ie Ain Bogsize + vinla t ppJsiay,]+ 


+ MN ane, 十 Zoos , 
至 于 在 有 风情 癌 下 ， 几 于 变换 矩阵 日 ?极为 复杂 [更 式 (2-9) 门 ， 上 训 | 得 到 的 运动 胡可 也 糙 俩 
der, Atum. 
4 2 在 机 体 坐 标 系 中 的 质心 动力 学 方程 
AE (EH Schi Ru va BUER A S. Di O E EL k 8 
的 ) 旋转 角速度 o, Ltd ep e mn 但 按照 习惯 简写 


e 
Wao., Oya Wo 


D 一般 公式 (1-37) 中 ， 取 机 体 坐 标 系 作为 活动 坐标 系 ， 则 写 出 


好 1 C x bS Ux D 
d . _ 
m dt Da: JE] ox Dyr = 
CG) 


Vae z Uzt 
DEM -Q 0 
-| Pinge L d Y J- d nl (2-30} 
0 PA 0 
玉 开 后 ， 得 到 在 机 体 坐 标 系 中 建立 的 飞行 器 质心 动力 学 方程 


el Ze, Tobias T Uye )= Peosp p ~ Qrasacos H -- 


+Y sina ~- Zeosasin f - muaing 


du 
ol gi Losse sin |= Psing, + sinceros f + {2-319 


+Ycensa + Zsindsio f — mgeostiosy 











wl 2 H OUy T yV: j= - Qin f + Zens + 
+ mycosgysiny 
泡 论 有 风 址 无 风 ， 这 组 方程 都 是 适用 的 。 
在 某 些 特殊 情况 下 ， 可 忆 作 如 下 简化 : REAMER, M e AR 
Us; T Veose fl =y 
Up 一 一 vainaro f e= — D 
Uzi =vsinp =vf 
Jä. AEREE TEIR eg 
UR ronst 


v 72 ew 





于 是 只 要 取 式 【2-31) 的 第 二 ， 第 三 个 方程 


du ` 
mof 一 了 六 FOr Ox fej )- Painf, Ke Y, 一 m gentes y | 











wl SÉ Z Get OI )= Z, Fimocotsiny | 
SM, 、Z, 是 气动 力 民 在 缚 y; 2 上 的 分 量 。 整 理 后 得 到 
32 w,- Gen, H , -2 Lo . 
i8 z (2-32) 
D £ 
Fr =y E t mu 一 -costsiny 


这 组 方程 与 后 面 的 转动 动力 学 方程 (2-41) 及 运动 学 方程 (2-46) 可 以 联 立 求解 ， 得到、 
户 ，?、 旭 的 变化 ， 但 不 能 求 出 ?和 高 度 ys。 这 样 的 为 程 只 适用 于 攻 嫩 2 不 大 ， 且 飞行 时 油 不 
长 的 情况 ， 俩 如 长 机 的 横滨 机 动 飞行 。 

4-5 在 地 面 坐标 系 中 的 质心 动力 学 方程 


对 地 速度 v ZE mn ke E, SECOND EE SEET IO DÉI 
于 是 ， 方 程 (1-37) 成 为 


, Uad Pegy = Q 0 
m i Uyd ]- vd Pain", J d Y J- 一 d (2-33) 
Ved 0 之 0 


展开 上 时， 变换 短 阵 B: 卡 式 (2-13) EA, Tal ERSTEN NEE 
的 质心 动力 学 方程 
m Zi = Pleosgppcosprosd +sinpplsinpsiny - cospsindcosy)] 
一 Qrosp scos + Y (sinp isiny, — cosp sinb eop} + 
k Zirin seot, + cop ssingsiny,) 


! (2-84) 
+ Feos cosy, — Zcos siny, -mg | 


m "uni " PL- cosp pin peost T sing p eospsioy + sinpsingeosy)] + 
十 Qsinp cost + F (cosp siny, t sinp asing rop) + 
+ Z {cosp og, — siny sinf siny) 


EARR. EBRR RILA (2~11) ]， 所 以 得 到 的 方程 也 非常 复杂 ， 这 里 
PRE To 





4-4 在 气流 坐标 系 中 的 质心 动力 学 方程 
ERREG Po SREE v 在 气流 坐标 对 S$S, 中 的 分 电 为 


Usg v 
| " -| > | (2-35) 
Uzg 0 





Trëtt A S AET IEE e, Sr Trëtt SN, Di bep DE eo, rei a. 相对 上 了。 
MELE GEE o,ao oi Err Sin, (rd, Bn hz OH (Dei) , Gei rb, 的 
mn, Clz-1n. Io 


Da, a 0 0 Pa 
iele hE 
人 0 d 0 
ind t pa ' 
= | roastleosy P, + sin, d | (2-36) 


— eos rinya d. -F rosy, Å 


krit AA {1-37) "GI 


opi sinp, pa s U 
i . 
LU Ñ 


— Éin, €. ops, 


Peot y — LU r D 
el Én ` | Y | Mi mG | { 2- 37) 
0 £ 0 


把 些 代 腿 正 ， 得 到 无 所 条 件 上 在 气度 坐标 对 由 建立 的 飞行 器 质心 动力 学 方程 


ni "3 = [eosl e ! y pho Q- mgsint 
S 6 | 
mol in "mei: Ee )=Prin(e typl H 


(2-38) 
+ 区 一 mo D rar s 
no( -cof DP, ie — Stan, 0 j= 一 sf Ju prsia Ai bk 
LZ + mgeot anya 
基 些 文献 《例如 [4 ]) 采用 这 样 形 式 的 方 种。 但 这 组 方 穆 还 不 能 直接 解 算 ， 必 须 疼 于 把 
第 二 、 二 方程 作为 代数 方程 解 册 dydt 和 ed9yudt。 而 解 出 的 结果 与 在 航 迹 坐标 系 中 建 赤 的 
AE (2-29) 完全 一样 。 所 以 ， 这 样 向 波 竺 不 如 直 域 在 租 迹 坐标 系 中 建立 方 但 . 


$5 飞行 器 矶 质心 转动 的 动力 学 方程 


第 一 膏 84 得 到 的 方程 (1-43) 鞠 是 最 - - 裔 情况 上 作为 绝对 刚体 的 飞行 器 绕 质心 针 动 的 动 
DER. E DEM, A. M Rau, EM, SM ， 苞 但 表示 气动 
JESH DAA 


j G 
Woh -Ty- Tmo EE EE CN twos )- ， 
一 了 (o ae, )=M, +M., 
* 24 =" 





Zut. -mw tm m2) -1,, ( dëi "Dia: ) z 
-fya ( gei "dei )=M, +M (2-39) 
1 eo ET EEGEN 
í 
faa ( dos -wy JEM. +M, 


F 
HAJO., wy, w KITARE A8 o CREIL IK paS: ah TS EE) 在 $, 申 的 
Tms Las 7 DU an Igas Tante KITAR NA, DOE AN M RARER: Ma, My, MIA 
Maps Mopo M22 9llETt a Jm augi DERS Pg., ATE 5 ARL. IR 
ERT FERo 

Hak EE KTR AAR EMO, y Am ER p50 TERRAE 
程 简 化 上 成 为 


i, dox 十 (7 -一 fps 十 了 (oswa 一 Ze = DPM, 


to de, ` , 
fy Kë H(I,- da) EE Lean "A (EE (2—40) 
i da i a na 

© p TUT Osy lay loz owp T EM, 


为 了 德 于 计算 ， 要 把 卡 式 化 成 解 出 的 形式 


Zon, 1 ~ 
Zeene DIEN EM fe Ty owt 
SIM AM 


LU Inf Tt, enen, / 


fü 1 = ! 
= My th Me ste he Toy et , (2-41) 


二 (7 


Zo, 1 
We 一 Ié KEE TEE dry — 62)] | 


轴 对 称 式 飞行 器 具有 了 机 个 对 称 平 面 ， 不 仅 7 padan D MAs D, Är. ;所 以 
转动 动力 学 方程 下 为 简单 


Mere 
EK L (2-42) 


7 ais ww SM. 


其 中 ?二 /二 1 。 





当 六 虎 故 动机 转子 的 陀螺 歼 应 时 ， 方 程 12-40) 要 稍 恕 修正。 发 动机 转子 轴 社 御 近 人 羽 平 行 
了 于 机 体 纵 轴 x;， 因 而 转子 的 动量 给 内 有 分 量 态 , ,二 站 ww;， 这里/ ,是 转子 的 等 效 转动 惯 贡 ， 
oc RR: of, Hri 而 及 可 以 认为 ;不 变化 ， 故 d 右 x/dt 王 0。 这 时 ， 根据 第 -- 
章 导 出 的 一 般 公 式 《1-45) 可 以 断定 ， 方 程 (2-40) 的 第 一 式 不 变 ， Sa LS 加 -一 项 
"Luet BAD lu Orya 


8 6 返 动 学 方程 及 几何 关系 式 


除了 前 两 节 建 立 的 动力 党 方程 以 外 ， 各 运动 参数 之 间 还 存在 着 一 些 运动 学 关系 和 几何 关 
Fo 这 旦 也 是 为 求解 运动 方程 所 必 丰 可 少 的 。 下 面 就 来 建立 这 些 关 系 趟 。 


6-1 速度 分 量 的 坐标 变换 
把 淋 讼 分 量 由 航 迹 坐标 系 变 换 到 地 面 坐 标 系 ， 可 利 几 伸 阵 B 


所 以 


在 这 里 考虑 了 


等 等 。 


Dygd Ié 
| p$ | i | 
Uxyd IW 
den ` H 
de =vrosd eosy, 
YP ie | (2-43) 
| 
{ze , 
An ~ ee sing, 
dx 
dE "Hen Yad 


把 吉庆 分 草 由 视 体 坐标 系 变 换 到 地 面 坐标 系 ， 要 利用 第 阵 B4 


HU sed Vai 
-Bda 


Urd Uat 


H 
Er =y, 0 costeosy + vyr( ~ eosysinteosy + sinysing } + 
i 
wos (sinysindecos — cosysiny } 

d ' 、 
Eat e sin Gin e eg peost — De sinycosy | (2—44) 

den . a 。 ， | 
Tdi = Us eosygsin¥ + Uy Ccosysingsiný + sin yeosy} + ' 


en, einysintsin p + cos ees 








6-2 角速度 之 间 的 关系 
K 行 器 绕 质 心 转动 的 角 丸 度 矢 晶 @@， 在 机 体 华 标 系 中 人 三 个 分 量 ws wys wo RH 
面 ， 根 据 机 体 坐 标 么 的 形成 过 程 【 见 图 2-1) ， 中 可 分 能 成 : 党 轴 ys H A i WA Së 
xpe Hofiy LT aA E za dcop Aiva Mi sinpo E 此 有 
ca Bsiry $ 
(eel le] 
EP Écosp 0 
EFEX, #7 
o, P +sintg 
wy rees Ë +sinyå | {2—45) 
名。 = — us fatn pi + cosy 
HA, oy os 本 由 微分 方程 〔2-41) Dm, ke bn, d, Aer DA Lé EDR 
数 方 程 ) 解 出 ， 结 果 为 


i 
Wi zz, T Lët am T wsiny) 
d 1 | 
= rost Loi wg -wsiny) ， (2-46) 
' | 
dë, L . | 
dë = sin wy | 


6-35 几何 关系 式 
EEATT. CATAE oo, H. 8. ba PIURE E Triho 
当知 道 其 中 的 五 个 时 ， 可 示 通 过 几何 关系 求 出 其 余 三 个 。 
A REHA EZANLA: TARARE ERENER, MARHE ,包含 
Di, BiaAy, diy BIE Ay., Are 以 及 B44 包 售 ?,，B# 包 含 和 六。 
根据 坐标 变换 的 传递 公式 《1-17) ， 有 
B'B¿=B;B4 {2—47} 
Krbnbr B}, BABIA HA (2-6). (2-10), (2-2) 和 (2-4) ëm. W kN 
HARRUA 3 行 第 1 列 元 素 相 等 ， 得 到 
sin ff =[sinyaingcos({p — p) + eos yain(y — Pa) jost — sin peosdaing (2-48) 
Am TARRAA 2 行 第 1 wk. a 
sing = {Lrospsingeos Cé Y) — siuusin(g- p.) eos — cos poosthsind } Zero (2—49) 
同样 地 ， 有 


ek KE (2- 50) 
全 两 边 的 第 2 行 第 3 ARR, Aa 
siny a Lei nfeosgain A — cos ycostsincesin P + sin yeosteos A) Zoe ë (2-51) 


5E ng, P Ya KÉ Pto 利用 公式 {2-48) 、 (2-49), (2-51) 7 可 以 先后 RHA, 


人 KÉ? 


D 27 D 


TA (250) 两 过 的 第 2 行 第 ] PKR. Sal 


sinl = singeoaco b ~ ee Men Ba nie H — ai ue vun dain g (2-52) 
全 等 式 两 边 第 2 行 第 2 列 的 正 素 相等 ， 得 到 

cosy, = (singeing + cos peusteos ci Yeoh (2-53) 
TRA 3 iTA LACR, Al 

sinf y, — ps (sinyinaeos f ~ cospsin H iepel (Zodi 


zm, dro, On, AJAR (2-52), (2-53), (2-54) ， 可 以 先后 来 Uya 
多: 。 但 是 ， 当 ?。 值 很 小 时 ， 利 用 余弦 值 公 式 (253) AFERKA. 这 时 最 好 用 
公式 (D REE. 

DE TEST P RHA A LE RE D ERE e fE IIo 

6-4 ”其 它 美 系 式 

1. AALER APRE ReRe RA 


EERI A 
Ux; L (OSC? rosji 
| Vyr Lei D |- d -siuu vs fi ! 
D 


Uze siap 
册 此 得 到 
Uy! 
tat 一 o, | 
{2-35) 
o Hai 
sin f = R | 
vs JIE e) 
2. 出 地 面 坐 标 系 中 速度 分 量 求 4 和 好 ， 
D ed V roit voa, 
HF | Upd | | Ü L "| siu fl 
Vad 0 = coh sinp, 
DEIER 
所 以 {2-36) 
ind = Eui, | 


3. 6 和 后 的 计算 

在 傅 定 气动 力 第 时 ， 带 村 4 和 | 记 的 值 。 虽 然 可 以 用 数值 微分 法 计算 导数 4、 记 ， 但 这 样 计 
算 的 结果 赐 差 较 大 ， 跳 动 厉害 ， 而 且 逢 始 值 不 易 计 算 。 所 以 希望 有 直接 计算 的 公式 。 

机 体系 Dx,ay,z: 和 对 于 气流 系 口 xsyszs 的 转动 角速度 可 以 表示 成 驯 - 0, 【注意 此 处 名 
实际 上 上 是 @, 的 简写 形式 ) ， 由 本 章 2-3 中 的 论 泛 大 参 看 图 2-3， 可 见 .二 述 相对 角 速度 巾 洛 轴 
ya 鲁 户 和 语 轴 z; 的 组成。 将 上 壕 美 系 投影 到 气流 坐标 对 中 ， 就 是 

. 28 ， 


co DEE å 

DIE d IER 

TS mp di O, II Ê 1 DEER ! 
Ww,- È 0 Cg 


把 BY 的 公式 (2-5) fo, DARAAN (2-36) 代 人 于 式 ， 得 到 
| Oeae H — oying- f+ Lo, — ART 


ein ët ovcgeg - f 


~ O1ettpadein f -+ o psingsin f Lin, -deosp 


Psin + Pa 
-| d gn Éegeiz, + Hsinys | 
一 $ asf sin, + Bees, 
合 上 式 在 右 陋 边 的 第 2 和 第 3 元 素 征 此 相等 ， 经 过 装 理 ， 就 得 到 zi 和 有 的 计算 公式 


= ainat O yens — P ee Hem ie, — iny, 
| (2-57) 


t= o," -iy (o sosasin é -cousingsing — 
-~ #seoslsia y, + ees zz, ) 
通常 ， 对 和 的 计算 精度 要 求 不 获 ， 击 且 & 和 彤 不 天 ， 所 以 可 以 得 到 较 简 单 的 公式 
=g -Pe cosy ~ Bsinya 

Gier 一 | (2-57 ) 

公式 (2-57) 或 (2-57) 适合 于 与 是 迹 采 中 的 质心 运动 方程 相配 合 。 

震 质 心 运动 方程 是 在 机 体系 中 的 ， 则 采 用 以下 的 公式 (2-58) R (2-58) 。 取 方程 
(2-55) ， 对 时 间 求 导 








i 二 和 
一 GË ZE" dÉ a 
cu d Te 
Üzr — Yz: 
TE ELA 
考虑 到 
ETIKETT 
Sim a 
> Uer fer TF Har fan ` 
= costa E y = 
Un l 
` (2-58) 
$= 1 UDr: — Un. Webs RRE ELCH | 
eusp HN 
当 & 和 局 不 大 有 时 væve ， 所 以 有 近似 公式 
GEET SCHÉCK | ETA 


EINEN = Uribe MVE, 





§ 7 关于 运动 方程 的 一 般 讨 葵 


7-1 方程 的 体系 


以 上 在 “平面 天 地 ”和 “ 硼 直 大气” 的 条 件 下 建立 了 Rimm ch JP. 

飞行 虞 质心 的 动力 学 方程 《为 方程 ) 可 以 分 六 列 写 在 四 种 不 同 的 坐标 FS, MALER, 
机 体系 、 地 醒 系 、 气 请 系 。 蕊 经 指出， 在 气流 坐标 系 由 列 写 的 方程 【2-38) 经 过 周折 乒 份 好 
结 为 航 迹 悉 中 的 方程 。 丰 地面 坐 标 芝 中 建立 的 方程 (2-34》 特别 复 杂 ， 而 朋 不 使 十 线 化 ， 所 
以 通 常 不 采用 。 实 践 经 验 尖 明 ， 行 之 有 效 的 质心 运动 方程 是 建立 在 航 迹 毕 标 对 出 出 碟 (2--29》 
和 机 体 尘 标 系 中 的 式 (2-31) 。 

至 了 改行 器 绕 质 心 转动 的 动力 学 方程 【 力 儿 方程 )》 ， 则 具有 有 唯一 的 形 忒 ， 屿 是 企 机 体 玲 
标 息 中 建立 的 方程 (2-41) 。 

国 此 ， 飞 行 器 运动 才 裕 有 两 讲 不 问 的 体系 。 第 一 套 体 系 是 : 力 方程 在 航 迹 坐 标 系 中 写 
Wo HARTEILLE A PS 出 。 它 称 AE-L” A AT- GR 
是 ; 力 方 程 和 力矩 方程 都 在 机 体 坐 标 系 中 建立 。 它 几 做“ 机体- ILAR A a UETA 
两 套 体 系 中 所 需要 的 运动 学 关系 式 和 几何 关系 式 也 各 不 相同 。 表 2-1 列 出 了 这 两 从 方程 体系 
的 性 质 及 组 成 该 体系 的 各 态 穆 的 号 码 。 

这 斋 套 方程 体系 各 右 其 优点 ， 应 根据 实际 杀 题 的 具体 条 件 来 选取 。 


Si 方程 体系 的 特性 








keka | DN , 力 年 方程 的 DIE 
Sätz j BA et 系 
E ALE | KERRAK HERRE {2-24)， {2-41),(2- -43) ， (2 -48), 
{IIT} i | (2 -48) +(2749), {2- 51), {2-57} 
roaa W - 
DG | ALE AE j DEES E (2- 31), (2-41). (2-41 ， (2 一 18), 
CEE) | 8 (2-55), (2-58) 





7-2 迁 动 方程 的 方块 图 

为 陈 象 地 说 明 存 求解 运动 方程 的 过 程 申 各 运动 参数 的 来 龙 去 及 ， 副 了 两 个 方块 峡 : 
图 2-8 (对 -HT 体系 ) 和 图 2-9( 对 子 TT 体 柔 )。 在 这 些 方程 中， 操纵 面 俩 转角 5。 R), 
9。( 方 向 能 ) ，5。 (升降 能 ) 和 油门 开 度 8 是 作为 已 知 参数 。 ROAA. TRANTE 
量 以 外 .所 有 方块 左 过 的 变脸 帮 在 方块 右边 出 现 ， 这 说 姑 ， 这 些 方程 联 立 起 求 是 封闭 的 ， 是 
可 以 求解 的 。 

可 见 ， 飞 行 器 的 适 动 是 出 一 些 非 线性 - 阶 微分 方程 以 及 代数 方程 所 描述 的 。 为 了 得 双飞 
行 器 的 运动 现 律 ， 必 须 川 数值 法 求解 这 些微 分 方 各 ， 并 辅 之 以 一 些 代数 运算 。 求 解 的 方法 将 
EEA HE. 

a 30 e 
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图 2-8 GT 体 柔 的 和 运动 方程 方 其 图 
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区 2-9 TI 侍 么 的 竺 动 方程 方块 图 








7-3 特殊 情况 ~ 一 铅 垂 平面 内 的 对 称 飞 行 

前 而 刻 出 的 是 最 … 般 情况 下 的 飞行 峰 运 动 方程 ， 邮 飞行 咱 作 空间 运动 的 太 程 。 在 : 些 特 
殊 情 况 下 ， 返 动 方程 大 夫 简 化 。 霹 在 我 们 求 矿 讨 -种 最 重要 、 也 最 常 光 的 特 琶 情况 ， 妈 飞行 
器 在 销 重 平面 内 对 称 飞 行 的 情 现 。 在 这 种 情况 上 ， 如 果 选 段 的 地 面 译 标 系 的 Opx oyo "EIST ER 
KIPER EMAA. M KITAA A A B A e e By, 2、 Pss Yar Usps Pa Zp Oxs 
Oys Öxi Ou, Z, Ma, METR; Më, PIIRE, Tiaia R a 坐标 系 三 
EERE- À, H2- mbh i SH asp, 





留 2-10 EREE A HA Et 


Kiri EEE MARE KITI Ea A TE AR e a A Ee R D e 
DAANEN. 
HT 体系 的 运动 方程 为 


一 一 Posla 十 Pal- Q- mosiud 


H 
ali Psinfa + pp) H am grosd 


L, M.-Pe, 
(2-59) 
= dxo = veol 
dyo = ysin 
dë 
dt Dr 
KZ ER 
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TTIE SI ah ir 


cf ` . 
m TE OT rosty'p ~ ru + Ysing ~ ngsing 
H 
dn , . - 
- e = -MsU + Pringy I Qsing +Y coa - myeosy 
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loz l; f 
fi - Pe, | 


dt 
cn —— i 
SÉ SCT EE 
y 2-60) 
- P Tig siatt u, pco 
d 8 
dt "ir 
tuit 一 ii 
h Ue 


D EH Dk ,十 v : 


sing = "2 (vya =dYpfdt) 


7-4 几 点 补充 说 明 

1. 气动 操纵 而 仿 角 和 油门 开 度 

本 章 7-2 中 指出 ， 当 税 纵 面 伍 角 ax，6，、6。 和 油门 开 诬 6 为 已 知 时 ， 整 个 运动 方 税 组 
Sum, dm, II, bb, RIDE DL MEET 
是 步 见 的 。 往往 给 定 了 对 飞行 器 和 运动 的 某 些 要 求 ， 根 据 达 些 要 求 反 过 来 确定 操纵 古人 人 角 和 油 
门 开 谋 ， 问 时 确定 其 它 和 运动 变 避 。 对 改行 器 运动 的 要 求 可 能 是 ， 例 如 要 求 飞 行 器 作 水 平 让 线 
飞行 ， 或 者 以 给 定 的 过 载 发 给 定 的 流转 角 作 雹 侧 滑 的 等 速水 平 盘 旋 运 动 ， 或 者 按 某 种 叶 引 现 
律 跟踪 -个 挨 已 知 规律 运动 的 上 日 林 ， 等 等 。 这 些 对 运动 的 要 求 称 为 运动 学 约束 ， 也 可 以 写成 
方程 式 (参看 [ 8 1) 。 

2. 飞行 更 作为 刚体 和 作为 质点 

完 玫 的 飞行 占 运 动 太 徐庶 包括 质心 运动 方 穆 和 绕 质 心 的 转动 方程 ， 亦 即 把 飞行 器 当 作风 
体 看 待 。 以 上 的 万 程 组 就 是 这 样 建 立 的 。 

但 愿 在 某 些 情况 下 ， 当 仪 要求 计算 飞行 器 的 航 述 ， 而 不 美 心 它 的 稳定 性 分 析 ， 盾 电 知 这 
各 具有 相当 好 的 稳定 性 时 ， 可 以 忽略 和 行 器 的 转动 技 性 而 把 转动 动力 党 方程 退化 成 力 征 术 衡 
方程 。 这 就 是 把 飞行 器 看 作 是 可 控制 的 质点 。 这 样 能 使 航 秋 的 计算 大 大 简化 。 

3。 飞行 器 的 质 基 变 化 

在 计算 飞行 器 航 迹 叶 还 应 考虑 到 ， 发 行 器 的 质量 着 不 是 不 变 的 ， 这 是 有 由于: 为 了 使 发 动 
机 产生 必要 的 推力 ， 必 须 清 耗 燃料 ,因而 飞行 器 的 质量 逐渐 Rho ENARE Ru AE phi 
耗 的 燃料 质量 ) 取决 于 油门 开 度 ， 或 者 说 取决 丁 扒 力 、 速 度 和 高 亩 。 因 此 ， 在 整个 运动 力 程 
组 中 还 应 补充 一 个 质量 方程 式 

。34 ， 


= p=-/(0,) (2-61) 





或 KÉEN vs h) (Zei) 


EROEM tE h ER 

但 是 在 以 下 两 种 情况 下 不 必 委 出 质量 方程 式 。 第 ~…, 燃 料 请 耗 率 & 不 大 , 面 比 蕉 虑 的 时 间 
不 长; A, HATAR SEN Ki (期 发 动机 工作 结束 以 后 的 飞行 )。 

7-5 四 元 数 在 飞行 器 运 动 方程 中 的 应 用 


目前 国外 在 计算 和 飞行 模拟 试验 中 已 广泛 地 利用 四 元 数 《quaeternion) 代 赫 Fuler 角 ,这 
放 而 的 问题 详 见 附录 《四 元 数 芷 飞行 器 运动 方程 中 的 应 用 》 。 
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第 三 章 “ “旋转 球形 大 地 ”情况 下 的 
Kýr aDJ E 


Sch, ÆRA E-A RA ERE E) EARR Za kinns 
I o PETERE RIR RA AES EA R 0 e e NG San. 

kreuk A, EEDE RAJE Esh A o P AA A AR EM EE a 
为 了 芷 省 篇 幅 ， 儿 是 与 第 -: 便 相同 的 邦 分 ， 本 章 就 不 重复 叙 进 了 。 

在 坡 后 一 节 中 进 述 赵 高 可 飞行 器 的 一 个 特殊 类 型 一 一 弹道 式 火 箭 或 航天 器 运载 火箭 的 运 
ECH 

扒 导 运动 方程 时 ， 仍 认为 大 气 柏 对 于 地 球 而 车 是 前 站 的 ， 因 而 改行 属相 对 于 地 球 的 速度 
CHERE y) 与 机 对 于 夫 气 的 速度 【〔 对 空 速 度 u ) 是 - - 笋 的 。 当 有 必要 考 虚 风 的 影响 时 ， 
可 邦 照 第 五 章 的 方法 推导 有 风 时 的 运动 方程 。 


81 坐标 系 的 定义 
当 把 大 地 看 作 基 旋转 的 球体 时 ， 侠 究 飞 行距 的 运动 ， 宕 要 以 下 儿 种 坐标 承 . 
]-1 地 心 坐 标 系 
地 心 坐 标 系 O.x sizs《 见 图 3-1) HETER. RARR E HE 转 。 原 AO, {EARE 


心 。 轴 六 在 赤道 平面 内 ， 道 过 Creenwich 子 午 线 《 即 经 度 7-=0 的 子午 线 ) 。 轴 x。 也 在 赤道 平 
面 内 ， 遂 过 东经 90 (7 二 90" ) 的 子午 线 。 轩 z. 生 是 于 赤道 平面 ( 即 沿 地 球 轴 )， 指 向 南极 。 


1-2 泽 动 地 球 坐 标 系 


浮动 地 域 坐 标 系 Qxayszs( 见 图 3-1) BRAO KITI O. $ ye 沿 Hiit, di 
向 上 方 。 至 面 x4z4 为 当地 水 平平 面 ， 轴 xs 指 向 东方 ， 轴 za 指向 南方 。 根 据 A E e 
标 肝 可以 是 作 “ 东 土 南 坐 标 系 ”。 由 于 大 地 是 妹 体 ， 当 飞行 器 运动 时 ， 这 个 坐标 系 的 原点 和 
指向 都 息 变 化 的 。 

EE: 此 坐标 系 无 论 就 位 置 或 指向 而 冯 ， 都 不 与 地 球 国 连 。 只 不 过 在 定 艾 总、# y 
划 。、 日 ，?y。 等 角度 时 ， 此 坐标 柔 的 作用 相当 主 第 二 章 中 的 飞行 器 牵连 地 面 坐 标 系 的 作用 ， 记 
坊 两 者 采用 了 同样 的 符号 。 


1-3 DL, Sat ges 


MIRER SOA us, SERRO E ER atz EE ue 
ëmgi. Driträizkiskbëzidp "ep" A EMAA “eg 
ÆR” 


D 36 H 





8 2 各 坐标 系 之 阐 的 关系 及 将 换 短 阵 
2-1 浮动 地 球 坐 标 系 与 地 心 坐 标 系 的 关系 
这 两 个 坐标 系 之 间 的 关系 用 下 列 三 个 参数 来 俏 定 〈 药 阅 3-1) 。 
Pr: ia, RTR OONA 


Be EGIARI MIR E I9 P fA ro(= 7 R? 
6,371,000m) 15 KITEE A ZM 
r= rn (3-1) dE 
Hi mc EA: MGreenwich T F Sp $] 
UE GENEE fhe ELR za 【向 南 ) 


- 180°< /e180:。 

地 理 纬度 gg; Mën, nt, 
HERAT ELR one DD, 

这 三 个 参数 r 、4 、p 也 决定 飞行 品质 
PDC SO. vn, 中 的 位 置 。 
RERA 

Xe Z regain | 





(3-2) 


Ye = reospeosÀ 

Ze = -rainm 图 3-1 (ERER e 
(Re, xez PARER R Oxyz kR o 以 RTR AE 

一 次 移动 来 实现 《( 见 图 3-2) o FEO: xeyr Ahe (EER -A RAD iy za 


Fe 





由 3-2 地 心 坐标 驼 与 泽 动 地 妹 尘 标 驼 的 关 和 对 





Séit, AG go WE Axiyri: Bairy PAE r ， WEO xayda. 
ei. GARA Gi-Auun AEE, JMA AEG ts aR ETE a) hb ER Atr AI ERE M STI 


B!=B,(-p)B,(-4) 6-3) 
A (1-15). CAD SRA RA LR 4 
ve A - sin 用 0 
Bi 一 | nog gin viedfcn- d. — siuy (3-4) 
sind: (ni Site os A Cosi: 


2-2 ”其它 坐标 条 之 间 的 关系 


IcCahkka (nik, RUE. zk, SC. ESRAR) Zopp, Diät 
d. WMA. ëm, MR y, BRIAR AO, GR RE 的 沪 转 pu, Wafi 
P. fjale L. ME EE HAE EA EARE m ADART E l KiTA 
HEIE n e REHE E = 

Eih ERPE, MB, Bi, Bi B3, BL. Bii A CRRA E 
Hiao AR AIEE MKA o 

E Ep., 8 ove pa’ 0. mR DEARI EE Sien AiE. (0 bur 5 
WFA eDi KILA lr GI mE GEA) D. kl, D KITA aA He EnH 下 
HEER EL Ho 

Bhp IU Sabine ZHN Cem Re NSIi3. ofh, Ob 
的 大 部 分 均 与 图 25 相 同 。 





B-53 各 誉 标 系 相 直 关系 示意 构 





8 3 飞行 器 质心 的 运动 方程 
5-1 在 任意 上 坐标 系 中 列 写 的 质心 动力 学 方程 


现在 考 瞩 地球 绕 自身 边线 的 放 转 (自转 ) ， 但 不 考虑 她 球 绕 本 出 的 旋转 【公转 ) + BEA 
为 地 感 中 心 的 吉 速 度 为 雳 。 这 时 ,地 心 坐 标 系 〇 xsgeze 是 一 个 绕 定 辅 作 等 速 旋转 的 坐标 系 ， 
其 从 速度 外 6 的 方向 为 指向 北极 Khan? BD 10 ?rad/so 

届 a .和 iv ,为 飞行 器 质心 相对 于 惯性 岩 著 条 的 绝对 加 速度 和 绝对 速度 ，F 为 作用 在 飞行 器 
KEJA JR AFI, mI KERRE, W 可 以 直接 写 田 质点 动力 学 基本 方程 
d, 
dr 

在 实际 问题 中 大 们 关心 的 是 改行 器 丰 对 于 地 球 的 运动 , 即 飞 行 器 在 地 心 坐标 系 C-<egr2r 
中 鸭 运动 。 根 据 肝 速度 合成 定理 ， 有 


a 二 有 a， Lë, 8, {3-6) 


ma, =F 或 am =F (3-5) 


Hop R= H A B A E F: 
a, Etma, MAAE (与 地 球 因 连 的 ) Abt RAEE EE EFTA 
IREE RIA aE A LER mn) 中 的 相对 导数 


EH on 
Bony, E KITA LD teo p ARARE E o 
a EREDE, B RAA 
a. än texr+ ox oxr)} 
此 小 ai, 是 参考 系 原点 的 加 速 弃 ， 在 这 里 就 是 地 心 的 加 速度 ， 前 已 指出 ， 它 等 十 葛 : SS 
mp, Zem, 四 旦 参考 系 的 角速度 ， 在 这 里 就 是 地 球 自转 角速度 四 ni 
r 旦 质点 位 置 的 矢 径 ， 丰 这 里 就 是 矢 和 从 O.O( 见 图 3-2) 。 因 此 
a= Wpx(mpxr) (3-8) 
A E FE REER KEFE 9 AE TH P ME E o 
a 是 Coriolis 加 速度 ， 由 于 全 = gn, EFT 
a. =20pXYr i (3-9) 
把 式 (30) RAR (35), H 
ma, =F-m(a, +a.) 
或 ma, — EIB, Lä, 
其 由 S,=-ma,, SZ. -ma。 (3-10) 
分 别 为 牵连 慨 性 力 和 Coriolis 惯 性 力 。 
考虑 到 式 {3-7) ， 写 出 


wl Är ) = F+S.+S。 (3-11) 
C 
mae EGE E a E LE OENE 它 相 对 于 动 奉 标 系 口 .xctrze 且 在 相对 角速度 
e, Ro, ;都 可 以 在 Oxyz 中 分 解 


+ 39 时 


o. sonal loune t rK 1 





(3-12) 
ËTT TT E sk Í 
hi, j, KERRO xy M ERER i 
Bim (1-28), Jr 
ï D f HI 
Lo ) =( A 1: "Ti (313) 
- ÖY, y o drs den. däs 8 
Wi ( 3) "nat: "Ach (3714) 
R av, EAER O xur HAH R 
ELR (ON RAR (3-11) ， 得 到 
5 
ell "leen, Ir: S, ES. (3-15) 


ASA (1-29) PAE (3~15) 写成 矩阵 形式 


Dr > Gry Sesr Urz Fol Sert Sex 
BE Ury 1 Ory | | Ury J- Ep Soy tSa, | (3-16) 
Urz Dra Ure E.A, Fe 
每 次 说 明 ， Gm Tia, y RREI KRE Sak ZC ushpüäR, Yr E Kibri 
对 地 坏 〔 即 对 地 心 坐标 系 ) 的 相对 速度 ; 局 ,是 所 选取 的 参考 坐标 系 Oxyz 对 地 球 的 相对 角 玉 
RW. FF 代表 作用 力 ; $。 和 $,。 分 曾 代 表 由 地 球 族 转 引起 的 牵连 惯性 力 和 Corioli; 惯 性 旋 。 
以 后 为 了 书 寄 方便 起 见 ,就 用 符号 Y RE KAO IRIRE HIR ENE v 《省略 下 标 /)。 





3-2 在 航 迹 坐标 系 中 的 质心 运动 方程 
TERRAE ERRO yn EAER RoB BARE Ev ONERA 相对 速 座 ) 企 GE ke 
Zrbm ARD (w o oh, 方程 (3-16) 成 为 


d D ex +s U Fort Serr t ders 
EN 0 Bi eren | 0 J- Fynt Senrt Seya (3-17) 
0 CR 0 LEE ACEN 





下 面 分 别 赋 帘 各 个 部 分 。 
航 迹 系 相 对 十 地 心 系 的 转动 角 巡 度 @，-.。 可 表示 上 成 两 部 分 之 和 
Dio = Did E (3-18) 

Ona HALER TEP AR RRD AE, BA A EE e te 
等 十 

Od = sinbi, + Pecol ja + Ôk; (3-19) 
Oa HARRAH TiO Aa EE, HE -l pat HeRR, E 
等 于 


DNA 一 - àk, -piu (3-20) 
IA (3-19) 和 (3-20) RAR (3-18) : HEIER R AA. fF 
d. ainf 0 


DIbnek — 
| PE BEZ EZA 0 DI 0 | 
[| å -1 IW 


D A0 D 








Peeing + Àlsinfsingp + cosh sin, geg) — Øros cosp, 
peoh + À leosfsing — insta, cs) + psiaheop s | (3-21) 
6 一 heosy cosg 一 siný. 
为 了 消去 上 述 裘 迷 式 中 的 变量 和 和 9 (经度 和 纬度 变化 率 )， 要 利用 如 下 的 运动 学 关系 : 
rp = -v,a 5 veol sing, 
reos À =u, a7 vosh cosp s 
持 此 有 
p= -veosbsiny:/r 1 (3-22) 
Å= ueegsflonsg , /reop J 
作用 力 F 的 分 量 仍 然 与 第 二 章 中 一 样 可 以 写成 
Po Peosg - — D 
[= Lee zo, Je Y EI -mg | (3-23) 
Fah 0 Z 0 
牵连 惯性 力 $ ,为 
S, — mopx (Opxr) (3-24) 
首先 写 志 wp 和 Fr ERGA ENEE 
Opere 0 
H: 
DDprrec " Gp 
rye Xe reosg nd 
| m. D y. -| eng sai | 
P ae Ze rsing 
TA (3-24) FARER 
Seere 0 ses 0 siy POOog sint 
| Serye LA 0 | 0 | reosq cas A 
Sarac 一 Öp -0p rsing 
Sin 六 
= ma Bees [2 {3-25) 
0 
洪 后 转换 划 航 秋 坐 标 杀 : 
KEE sin A 
[se |- seeeie: E | 
Serzh 0 
singoeos® — costsing ssing 
— mo Ärensg | vos oos sindsiny sinp | {3—26} 
coa% sing 
Coriolis WIENS. 3 
$.=-2m@pxY (3-27) 


《注意 这 里 Y 代表 相亲 速度 ) o EE Eh APRES 





Sars 0 、、,， 
EJ 1 
Serzh — wp D 
0 
= Zanen | conp eos {3- 28) 
cosdsingp — sin fl sing cosy 
Ikik ajA M EJE gay% ik Hit 


2 
BEEN rD ) (3-29) 


Bin g 59.8204m/3; AHIRA MH Bi D me, 
rp=6371000m ， 为 地 球 的 平 瞎 半径 。 
JEA (3-21), (322), (3-23), (3-26), (3-28), (3-29) EA Ñ (3-17) H; RF 
EU i BAE aE Ae br A E a R E AF ALG a DS A fE 
Ku Kë rela + yrpyemsph - Q- mg (2Y SEA { D } + H 0 } -i 
- {ono reonip sinf eos -sos sinpasinp) } 


Ié 
mU TE = STEE ke ge, Jain Heintz, Laint g + op)cosy ,d+ YF cosy, 一 


Soe, wl Dë ) eo V+{ 4 mo + | 
siny, mY Ty cos \+ y D rg , (3- 30) 
十 { Hamme peasy cosp LA { mo greosple seus p + ein atu 9 sug) H | 


dý, 
(dl 


- 0 | 
一 mucas = {PL ~ os(at gp)sinfoos ys + sing +ppjsiny:]+ Ysiny ,+ ' 


. 1 | 
+ Zeosvws} + Leier "gent ‘gr)+ { 2muo plcos Osiny' 一 ， 


-sinfsinparosp)} + { mo Sreged ein Ges, } / 


为 了 说 胃 各 个 内 素 对 运动 方程 的 贡献 ， 把 上 面 每 个 方程 的 有 迅 分 为 四 部 分 ， 用 世上 标 
出 。 第 一 部 分 是 直接 作用 在 飞行 器 上 的 力 ， 第 一 章 (平面 大 好 情况 ) 的 方程 (2~29) 购 有 边 
正好 就 是 这 妆 分 ;第 二 部 分 表示 地 球 曲 率 玖 影响 ;第 三 部 分 是 由 地 球 旋转 而 引起 的 Corioli; 惯 
Hi, mn Ed 地 球 旋转 而 引起 的 牵连 惯性 力 { 芭 离 心力 } 。 这 由 部 分 的 顺序 E 
是 按照 它们 的 重要 性 的 顺序 排列 的 。 在 本 章 呈 ， 将 分 析 这 四 部 分 力 的 作用 的 量 级 。 

因为 飞行 器 的 位 吊 要 用 距离 " 、 经 座 4 和 纬度 "这 一 个 参数 确定 ， 而 且 在 动力 学 方 称 
(3-30) 中 包含 变 基 4 和 9， 所 以 还 需要 补 光 -- ech ES 求 描 让 这 三 个 参数 的 变化 。 为 
此 ， 把 滞 度 Y 分 解 到 泽 动 地 球 坐 标 系 ， 得 


dF “Uyd ~ ysin 
re = =y, a5 Ons sing, (3-31) 


TA 
reog Gi -= Upd Uo cosp 


H 42 M 


im 


5-5 在 浮动 地 球 坐 标 系 中 的 质心 运动 方程 
现在 选取 浮 芭 地球 坐标 系 Dx vaaze 作 为 获 考 坐标 系 。 于 是 方程 〈3-16) FR 
J f Lya JU 一 es DÉI DU, d Fedt Serxd T KEE 
SEO Uyd L Did + yd | Uyd 1 Fyd t end te, kd | (3-32) 
` ad CO dæma ze Uad Fat KE T KÉ 


FERAS BRAA o 
EPAR AHI) Tib AAIE A AAEE EL (3-20) 式 ] 








Oae ™ AK, -pia (3-33) 
KEER RA 
的 有 一 wx 0 — 0 — ġ 
| TEE ! 一 dl 0 站 0 J- Asing (3-34) 
(D d—e: zd — A D 一 Avon 
"o | (3-35) 
Los Än, A 
doerxd ËTT: 
Wir | de yd D Ux dtgE tr (3-36) 
Lee ed 一 Tt 
作用 为 FF 的 分 景 为 
Fg Deosgp yp 一 Q Ô 
| Fa =Bs | Wë 十 Ch | Y Jf -natan | (3-37) 
Pa 0 E 0 


牵连 惯性 $3 .= 一 加 px(@pxr)。 由 于 


Wp-xd 0 0 
| Op- yd el 0 -| sing ) (3-38) 


WD sd -0p — ro 
red D 
HD | fyd |-| r {3-39} 
Fad 3 
KEE d D d sy 0 
Cd | Seyd |--03| sing | sing | r 
Serzd — Cosi -- cosp 0 
d 
=mobr | egal ` | (3-40) 
singes 


Coriolisi tE Dë. = -2m px* 的 分 量 为 


KEE ð aay Usd 
| Ad Lal sing | — | 
A -eop Uad 


D 43 D 


Uy neg EU, däin 
~ -2mp = Uy gtu {3-41) 
"e dl DÉI 
把 式 (3-36), (3-37), (8-40), (8-41) 代 人 式 (3-32) ip, ATEA EaR Ae Sr 
ri KITA Oba UE p T: 


du, ` ， : 
mai d =f PLeosg peos Ve kaing elatnsinie < rospointeosyp)] — Qeop cost +. 
+Y (sints saat, -cosp ineo + SZ (inf, rou, + 


Hi 
Foorg singniny l +{ r { gEV, EE varena) JA 





+ { -2moplangt a + cospu k ti) 


dv 
nm" A7 AJ Pleosg psing Haing penste os ) - Qasin |! Hau, Hues, — 
dla 


ao ' (3-42) 
~ Zoo 一 DD E ) p[i eta tot) + ， 


+ { Zentren, } + {mopreo g } 
Tv, 
le = { PĪ ~ cosg apin Gees + sing pleospsiny + at gin fege ) | 十 


十 Qing aeos h + F (cosp sinn, + sing asini cosy s J+ Ziro acosya T 
-aing -indinp.)} +f (euni 一 ul 
+ {2mopingu,g } + Imeärvgsgsirg } 
Sr (3-30) 那样 ， 在 方程 组 {3-42) MA TRAR., HI 上 分 成 四 都 分。 第 
一 部 分 中 直接 作用 力 ， 第 一 部 分 表示 地 球 曲率 的 影响 ; 第 三 部 分 是 由 地 球 旋 转 引 起 的 
Coriolis 惯 性 力 ; 第 四 部 分 是 由 地 球 旋 转 引 起 的 牵连 惯性 力 。 


因为 最 后 必须 用 球 坐 标 r，4，FF 博 定 飞行 占 相 对 于 地 球 的 位 置 。 所 以 需要 补充 下 人 列 送 动 
学 方程: 


dr | 
di "Zeg 
de 





roj Eed (3-43) 
dÀ 
reosp qg Urd 


EOTF EI REENA AAA A e 





zd i 3-4 
eg. — — Uri { 4) 


v= Nuitvtdt vi 





3-4 飞行 距离 的 计算 

按 以 1 方程 计算 的 结果 是 球 侍 标 r(7) AG), G); 它们 能 单 值 地 箭 定 飞行 器 相对 于 
地 球 的 位 置 。 当 计 工 飞行 雌 的 航程 (或 射程 】》 财 ， 还 要 求知 道 立 行 器 离 起 飞 点 《 眶 发 射 点 》 
的 距离 ， 这 跑 离 是 拨 在 地 球 表 面 上 的 铅 荆 投影 点 米 计 异 的 ， 
而 地 球 甫 面 上 商 点 之 间 的 拭 离 是 指 通过 这 两 点 的 大 大 的 骂 
长 。 

设 : 起 习 点 为 名 (45，oo)， 改 行 器 在 地 而 于 的 投影 点 
APC. ei, RPOAK (A34) o 

为 此 ， 分 别 通 过 外 点 和 请 点 作 子 午 线 QN MPN, æ 
六 是 厂 极 。 在 球面 三 角形 名 亦 疡 中， 


- ba T 
ANSA- Ze, 进 p=- "Wan A 9 一 -有 一 人 








图 3-4 球面 三 角 祖 中 NP 


PARURA. ARIEL na 
Coen = cospeosg + sin psingeos fY 
所 以 cost = sinp asing + rosg pcospuos( A = Aa) (3-43) 
w KITEN POAR) s 则 为 
Stan (3-46) 
HKE sih (PaL FAR) 时 ， 公 式 (3-45) 将 造成 相当 太 的 合 天 误 状 。 注 意 
HO ijn, 一 mo 和 /一 4 者 是 很 小 的 角 ， 式 〈3-45) 可 以 变换 成 


1- 





L n? esing osin® + cose yeosg i1 EE 一 HOM 


1 .. ， 
sl — ga) 一 cosp ocosp H(A — ża)? 





=] Ly a)? eos? > Coto Jog A-Aa? 
由 此 得 到 距离 的 近似 公式 
S rpt (p - po)? teo ep CA 一 Zei (3-47) 
其 中 引用 了 平均 纬度 


1 , 
Geier (PrE) (3-48) 


8 4 KFE EK 
仍 假 设 飞 行 器 具有 纵 同 对 称 平面 ， 因 而 7 ,二 7 了 := 二 0。 转 动 动力 堂 方程 仍然 是 


dor = let, EN EN, -Fold td, 一 TI. Jee 十 ! 





dr 


+F- Tlg- Iloro] 
à 45 à 








do, = EE -JÈM >a DM, -ig EA 一 i, — leen, t 
HESE EEN (3-49) 
a | 
[EM (fy DEE | 
z bh 


dou, 
dt 


ZE, A Hëlen, wys uE KITE UE 对 角速度 名 在 机 体 Air 系 中 的 分 最 、 EO 
TEEN ECKE CH 
SEET 
O=O Oudt Oe 
Ate, ERLER RAJ ES ARE. MEFR A HARE op 
A. HAK (3-20). A 
eu O- Oder Gase: i) enk, 


把 上 式 分 解 到 机 体 坐 标 系 中 
0 


Ox Cr Ki 
EHEJ i Jaf 
0 


E pda sl WÉI Aen 
展开 后 ， 得 到 


pader: Tuet CÅ- ouni iodsing — malen gees) Ggs Bes d 
Ordet 5y (A pE — ceosypeostsingt (eospsingsin g + 


+ sin prosperos] + BC -cosvsindeosp + sinysiny} (3-50) 





(dat Ty] (让 一 opl isinpeostgsing + È- sinysintsin y + 
| cos gege eos ]+ (sin asi ufeneg + cos sing) 
第 =. 章 得 出 的 运动 学 方程 (2-45) 仍然 有 效 ， 但 其 中 的 os o wx* 现 在 应 理 解 为 机 体 
系 相 革 于 溯 动 地 球 系 的 转动 舍 速 度 分 量 ， 因 为 角度 、#、? 是 机 体系 与 洋 ann SS: pt f 
HE o 于 起 ， 运动 学 方程 是 


dr . 

di TER redex d -tahlo rade yr CSF 一 dr zrsinp) 

dä 1 
"in oop Gidrat a y COSY T O ed a rainy) ， (3-51 
dù ， 

dT LD taga z COSY 





LEEK. IHI KI RR DNA ZG R E REA (GAD Rhos wr os PER 
(3-50) Ro: ye RERA (3-51) Ry, vo do 这 是 严格 的 处 理 方 法 。 
但 实际 上 ， 太 究 飞 行 器 的 转动 动 万 学 特性 时 ， 四 p 和 是 -都 远 远 小 于 飞行 器 转动 的 台 
(BU e, 1 ， 内 而 可 以 到 
O; ag DN 


Fie aaa LGS (3-50) 和 (8-51) 可 以 简化 成 
NM 46 h 


dt "Gs" (ell yeop — riny) 


dy 1 l 
di = coah Co pensy = @ siny) 


(3-52) 


dô 

di 
Aa. TWFA EA ai ER a OURA EER A R e au SE, RH jem (3-49) 和 
(3-62) ， 与 第 一 章 一 样 。 

ERB ve, ont. 区: 以 后 ， 可 利用 第 一 章 6-3 中 的 关系 式 【2-48) (2-49), (2-51) 
计算 #8、a 和 y,。 在 必要 时 ， 仍 可 利用 该 童 6-4 中 的 关系 式 (2-57) 计算 A Ma. 在 “旋转 球 
展 大 地 ”情况 『， 严 格 地 说 ， 所 用 关系 式 中 的 wss w o: Æ X ML iR E A THRE y 
ordres Ordregrs Or doaro 但 如 前 所 瀑 ， 可 以 近似 地 取 为 绝对 角 速 庶 分 量 ， 

若 质 心 运 动 方程 是 在 浮动 地 球 坐 标 系 中 给 出 的 ， 则 式 〈2-57) 中 的 6 和 细 , 述 要 利用 由 式 
(3-44) TH TARHI: 


二 wyi y T w atr Y 


Vyd Upd DI dÉ FI | Heda ln 





d: ie oe 
ÉIER 
d, "eps fg, lege ider 
§ 5 Than iE R 
5-1 两 套 方程 体系 


本 章 在 把 大 地 当 作 旋转 的 球体 的 条 件 下 建立 了 飞行 器 的 适 动 方 称 。 

象 第 一 章 那 样 ， 木 章 也 采用 了 两 套 不 同 的 方程 体系 第 一 套 称 为 “ 航 杰 -机 Eë f 
Sg nur, Com, kB adr 系 中 建立 ， 力 摔 方 RER, 第- 
套 了 你 为 “地 球 -机 体 体系 ”， 缩 写 为 DT。 在 这 个 体系 中 ， 力 方程 在 AREARE 
力 生 方程 在 机 体 华 标 系 中 询 写 。 表 3-1 总 结 了 这 出 赛 方 程 体系 的 主要 转 征 。 


表 3-1 方程 体系 的 特性 


力 方 程 的 JEB | 3 a 








i 3-30), (3-31}. (3-49). (3-50). (3-51), 











E R | , ! C2-48), (2-49), Gah, (2-57) 
(HT) i HEERA HERR R REH EE 
` g ,9-42), (3-43), (3-44), (3-49), (3-50) H 
上 地球- 机体 (3-51), (2-48), (2-49). (2-51), (2-87) 
{oT} FaR  ， DR EE {3-50 
OSTEN EE EE KEE MK 





5-2 运动 方程 的 方块 图 
图 3-5 和 现 3-6 是 这 两 套 体 系 的 运动 方程 万 块 图 。 


Ör, du i y da 
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AFE Ca 207 
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atA E 


{3-31} 


dan ipta d. 





o ah KÉ sel — Gs 
DA? 转动 动力 学 





ey p, to 
Hai 方 释 C Ia 
Hr 



















Bis bhy is 
`: ZIER 
sf 
WU 
rh (3 ~ 207 ba 前 
P AE 
t3 = 52) 
it 
oau, : —[ DE 


运动 学 :关系 式 


(3: A1? 





几何 美 系 式 
{2 = 48} (2 -49) 
EE 


DH VE) 


KI 


pr 
Bagi ZER F} 
PETA 《2 一 57》 $ 


图 3-5 LIT 体系 的 孝 程 方 执 图 


18 。 


Za Én an Én 
ITA? 
DÉI 

aa 质心 动力 学 
ixy py Gr 
Va Vyd se 方程 《3 -42) 
Tap 





UN 
Vat V Vas 


hP 


Fa Ha, Rat 


Ze siyi z BR 
Ae, p 转动 动力 党 


方程 (3 一 49) 


















运动 学 关系 站 





PA { 3-50) 可 及 近 航 
p ediz 地 用 式 
i = 52) 


i—i wl 


(NIE DEM ` 
CH 运动 学 关系 式 
e 
d t 3—51} . 
, 
几何 关系 式 
Pth, nt { 3—48), © 2—49), as P, Ze 
【 23—51} 
Pry | EES 
DBnps (3-53) č.p 
FaF h, at 2-3) 
Beil nF 


ae 1 和 棱 系 的 方程 方块 加 


5-5 特殊 情况 之 一 


在 芳 虑 大 地 的 曲率 和 旋转 的 条 件 下 ， 现 在 来 听 窜 -种 特殊 情况 ， 即 在 杰 道 束 面 内 的 对 称 
飞行 《{ 风 图 3-7) 。 在 这 种 情况 下 ， 纬 度 w 二 0， 飞行 器 的 所 有 横 侧 运 动 徐 e a, 
Pes Br pa Osr Opa A An Z, Ma, ij 一 -均等 于 霉 。 因 而 适 动 方 种 大火 简化 ， 
不 难得 出 如 下 结果 


Lé P Y 
Deg A 
Va "Li 8 N " 
Vë A 
Y li > 
hyj 
d 
SA 
ai 
T 


Xa 


D 
"CM 
us 
Dr 
图 ?3-7 在 赤道 平面 内 的 对 称 飞 行 
TT 体 系 的 运动 方程 汶 


` 
di = Pista bpp) -QO my ( p ) sinf + me Braing t 


dÉ 2 
muyy = Psin(a >e) Y -mgo ( js 十 1 "mue + | 


— 2mvw p Moreo 


d 
Sg gien (3-534) 


A 
r d4 We ueosh 


» K >» 


dor, 





f di =M. 
dt 
di Si TAT UD 
a=ġ- 0 
pT 体系 的 运动 方程 为 
bi di | fos (pt KK 一 Oeos — Ysin 一 m Dd T Mme pya 
v r m 
Tur = poin(gptd) -Qeind + Yent — mg ln Set 
— Imap t Mohr 
dro 
BEA y 
r dt xd ER E — DE: E A (3-55) 
由 DE ^ dÉ 了 ah X k ; 
ii = Gate ` ` d’ R ' 
x ` ` - i K j 
` — DN kk H futa cob) | 
v=y uia T Upd D ~ 了 | 
ie m, EE 
: ， EE 
t sgeig n Esg eil 7 
th bR f "vt teh WW n , 
= 4h-niat “一 
a=i-ø 、、 ak o 7 


5-4 特殊 情况 之 二 

在 某 些 情 况 下 ， 地 球 族 转 的 效果 远 小 
于 地 球星 率 的 落果 《 详 见 本 章 5-5 中 的 讨 
IE ， 这 村 可 以 采 击 善 虚 地 球 曲 秦 但 忽略 
地 球 旋 转 的 运动 方程 。 同 时 ， 还 假定 飞行 
器 在 铅 重 平面 内 作对 称 飞 行 。 中 于 忽 覆 地 
球 旋 转 的 效果 ， 这 样 的 对 称 飞行 可 以 在 任 
何 大 圆 乎 面 内 〈 不 一 定 在 赤道 平面 内 ) 实 
现 . 这 就 是 常常 遇 到 的 桂 殊 情 况 之 一 。 

Ah Kiria o Hb Ei aJ LRI E 
标 Cre a) RE { 见 图 3-8) 。 飞 行 器 的 
运动 方程 可 以 由 式 (3-54) 简化 而 来 ， 只 
SS Rb ën bn, IRESE Ise D er 








图 3-8 在 任意 大 见于 面 内 的 对称 飞行 


e Al 








du 2 
morn = Prosla tpp- Q-Go (2 ) sing 
i r 
dü KR | 
MU gt — Priula tpp l F -mg ( 2 ) rost t ` muiecosff | 
do 
d “二 gens 自 , (3-36) 
lwa m | 
Ti E, | 
dd 
jt D 你 
他 -下 
o-ë- Å 


5-5 关于 地 球 曲率 和 旋转 的 影响 的 估计 


goen, RR (把 大 地 当 作 施 转 的 贺 球 体 》 的 质心 运动 方 各 (3-30) 成 (3-12) ， 
与 第 二 总 《把 大 地 当 作 平 面 ) 的 相应 方程 (2-29) 或 《2-34) 相 上 比较， 在 方程 的 右边 增加 了 
所 部 分 ， 第 一 部 分 是 由 地 球 曲 率 造 成 的 ; 第 二 部 分 是 由 地 球 旋转 引起 的 Corioli' MHE J: 
第 二 部 分 是 由 地 球 旋转 引起 的 牵连 惯性 力 ( 离 心力 )。 现 在 米 分 析 这 三 部 分 的 影响 上 大小， 

为 了 明 最 起见， 考虑 -- 个 特殊 情况 一 一 在 赤道 平 量 内 的 水 再 飞行 。 这 时 9 一 0 ， 方 各 
(3-54) 的 第 一 和 第 二 方程 成 为 


A bi eu 0 


0 = Psinte rep Y- mg, TEE i (3-57) 
-| Bm nt Mmeo Br | 
可 见 ,这 时 三 邵 分 附加 为 仅 出 现在 法 身 力 ( 即 第 二 个 ) 方 程 中 ,地 球 曲 率 的 影响 表现 为 m2* /7 ， 
ES TT GEI, E 离 rle: 田地 球 旋 转 引起 的 Coriolis 惯性 D'äZmzzun Dë 
KSC, mies: CL, 由 地 球 旋 转 引 起 的 牵连 惯性 力 mob 与 距离 + WEI, MER 
ER. 
BE KATERE D300km; KIPIR +( 指 相对 十 地 球 的 速度 ) 分别 等 于 3800m/s，1000w/x (大 
PAI M =3) ,3865m/s，7730ra/s (卫星 国 周 轨道 速度 ， 即 第 一 学 审 速 底 ) ， 取 
ons 7,2921 x 10 F rad/s 
rp=6.371 x 10°m 
930 二 9,8202m7s2 
SI 
r'=rpt+h=6.671 x 10°m 
9g=g0lrpir)’ =8,96m/s 


ki -Sia Mii DHRU] mgle {以 百分数 表示 】 和 加 表 3-2 所 示 。 


表 3-? 三 部 分 附加 力 的 大 小 《示例 ) 











1 
EX E 地 球 曲 率 影 响 | Coriglis 惯 福 力 影响 FERE y EN 
i (== e ' 2mvo pimg | morig 
(ufs) ro? (X) | 
| WK : (%9) 
300 ， D.15 0.49 : 0.40 
1000 ， 1.67 1.63 n. AN 
3865 25 i 6.29 i 0.40 
7730 } 100 ， 1 6 0.40 








nE. IFEA KI iR wenig, EA Coriolis WEM: Dag 
EED IE: FERIE UE h, TAAR. 

EIR, TEA TERRAN g Mie DAR NEE. EEGEN 
(Pint JIDD MHK ENIGER” Ep biH A S mk RE nein m 
速度 的 人 台 成 结果 。 当 使 用 这 样 的 重力 加 速度 数据 时 AERE EE A a To Im P 
程 帆 不 应 再 叔 虚 这 一 项 。 

AAJA LS Sri E T A H ENR SETOMAA ER. HEERE MET E KITE 
离 和 飞行 时 间 ， 而 是 & 行 速度 。 这 个 结 诊 也 许 是 与 人 们 的 直观 印象 相 予 扶 的 . 


$ 6 EE E 动 方程 


6-1 概述 


显然 本 豪 前 儿 节 所 古 记 的 运动 方程 ， 原 则 王 适 用 十 任何 类 型 的 超 剖 速 飞行 辣 ， 但 由 于 强 
道 式 火箭 的 调动 段 具 有 其 些 特 点 ， 相 应 地 要 定 
dëst p Sachs, Si eso 
运动 方程 。 

弹道 碟 火 科 的 主动 自 是 水生 发 动机 开 作 时 
的 飞行 贞 。 它 有 以 下 证 要 特点 ; 

lL. 以 铅 生 状态 发 射 ， 并 且 在 发 射 后 有 -- 
Dirk kE, Hi TER T E AE 
JARRE; 2. TE È SFE MHRA TF AS 
AMAVA ELLE) 都 很 小 ; 
2. 主动 段 的 水 平 距 凯 不 大 ( 见 R3-9). Æ 
节 所 建立 的 运动 万 程 也 适用 于 航天 器 的 运载 火 
葡 的 飞行 ， 因 为 它 旭 具备 土 述 等 点 。 

在 火 葡 发 动机 熄火 以 后 的 飞行 腿 ， 称 为 接 图 3-9 弹道 式 火 箭 的 主动 自 
sO ERRIA A BRIR. 


„7 


一 卫星 
d 





6-2 坐标 系 和 变换 矩阵 

L MEERI GG br Sne 

AER h M EIRO sx Yoz EE IHR, Dap Eer, Ré, Jg Sin, fr 
发 则 点 铅 重 向 F。 轴 xs 在 射击 平面 内 ， 水 平 向 前 。 射 击 末 面 即 是 由 大 射 点 、 日 标点 和 地 心 
质 决 定 的 平 而 〔 见 图 3-10)》 o hr KE, EATA GFE. KEB OD 的 加 迹 - 一 -暗道 
kt DAS E Aan z, Miao 





图 3-10 HEBRAEK eA 


ZE kb ke BO, 原点 口 在 火箭 质心 。 畦 < 人 zy Et T ie, 1 
图 340) o HA RLE PRA J bio EE E A i Apr AAE E a h i A pr R AE E y Ha ii 
Hh fo HEOx yie BAR AOX, Yg. 

DS bb Ce TALER RO. X yeze (ERESI) WER H RIR 
ORIRE EMBE DEE EIERE 

地 心 华 和 标 系 口 。x-yeze 首 先 绕 轴 ze 转 过 - Aa CR -pA Aa a 转让 ga 
(BIS - x“ 转 过 pu) ， 插 且 作 平移 ， 就 得 到 通过 发 射 点 的 东 (x) .上 (yy) 南 ( 2 ) 瓦 标 系 。 吕 后 
DS Ain, at Ia Wb EG Ee 


Ze 


A 


W Zr 
3-11 MAERO gx yyg2g 与 地 心 举 标 系 品 cxogeze 的 疯 系 


IC, IEEE $r E St kk Spe ée ue 


a gjd a 


”cospo sing 
Bs=B,COB -POOB -A mings (3-58) 
“ Oso OS 
这 显 仪 号 出 了 以 后 有 用 的 第 3 Un vs, 

注意 ， 固 定 地 而 坐标 系 和 地 心 坐 标 系 各 是 后 地 球 固 连 的 ， 它 们 之 间 拜 次 有 相对 运动 。 

2. 机 体 坐 标 系 及 其 与 元 面 举 标 采 的 关系 

HE GEHIO ERRO z MEDEE TAM METES AER 
e ECKE GE EE CES ARAE e e CEET -12 
Hel EF É xy LRE e AERA AmA o, 向 上 为 应。 机 体 的 
Ah, 与 精 纺 之 间 的 角度 称 为 滨 转 角 ”， 这 里 的 轨 凡 在 平面 xwgye 内 ， 垂 让 寺 轴 + IW SE 
?为 正 。 

三 个 涯 态 角 之 所 以 这 样 定义 ， 是 因为 : 第 
一 ， 上 粹 前 定义 的 偏 航 角 水 适 几 丁 改 射 后 的 铅 王 
Kirit Plán: AAi RR TI, 
H. Camus a, Æ 
按 以 前 的 定义 ， 则 偏 航 角 为 90 ,显然 不 合理 ; 
第 二 ， 现 在 这 样 的 定义 与 弹道 式 火 箭 中 陀螺 
ER EEGEN 
献 [ 11]，f123) 。 

DR 3-12 q m., iri M dh B ERR 
Dxrsesr 册 发 ,首先 绕 轴 zs 转 过 角 中 到 
达 Ox yz'， 然后 绕 轴 yy 转 过 角力 、 到 达 
Oxyz’, Rai t tAr. 就 与 机 体检 
br SO, Eë, SES ei ArbztO-1e 





RE EE EE moi 机 全 从 标 夭 与 地 面 尝 标 对 的 闫 和 
IZ voa cosh 1 cosy sing Po -any 
-cosp simf — : va eost + siu% coy 


B;=B.(r)}B,(4)B.( =| ! ain4 sin peost I +einy sing sing : {3-59} 
| siny sinf 1 : -sinp rost + 
| H H 
' + rosy sing cos i -+ eo3y sing sing 3 


KU EN RE A PE EREA CN RE aT AME 
rost sing 一 功 
| 


CP c0» D 


1 


-sinf coi y (3-60) 
psint t peos? ` — jeos psing ] 
3. 气 沪 玲 标 承 及 其 与 其 它 坐 标 系 的 疯 系 
气流 华 标 条 Dxoy。zs 的 定义 也 与 第 生意 糊 同 。 但 它 与 地 面 符 标 系 之 间 的 s 衣 度 也 应 重新 定 
Xe RERE v《 指 相对 于 地 球 的 速 庶 ) 与 xoys 平 面 之 间 的 角 亩 称 汶 的 下 偏 转 md, Wi 
REHE Së v ETH x oye ERTER- fhaoi 的 角度 称 ADER, DERN: EE 


D 55 + 


IE. a, Eh Hau" al, lët pr Dm gäe mp, CW Ei2181 。 

i.e, rz gr täkniä, e yo me, ft, t iR i mp, Il bet 
Bizge kpo Denke 

eash >in e, 
| SI cosd rs (3-61) 


od rs — in 1 
D 


WÉIEE g 
ROx yz LARK, PMRAM 
ou E E E E du a2- 
中 式 (2-6) ; HE Eja, PiE 


Eo ie 可 APRS ETICA ERAR EE 







Ze 


la -ÊP 
"7 o 1 0 | (3-62) 
#0 i 


1， 各 坐标 水 AAR ag 
以 上 从 个 伍 标 系 之 问 的 关系 可 以 
H-A RR E A3- o 
a 儿 何 关系 
为 了 找 出 nO ep, Pa E 
ZARLA. 利用 如 下 组 阵 等 成 H 
Di-D'rB: 


SES Oreyee 


Re (el 


ELORE A Üo Ae kee 


图 3-14 弹道 式 火 第 所 用 的 华 标 又 关系 简 图 


FEA (8-62), (3-61) 和 (3-60) 代入 上 式 ， 芽 展开 后 ， 得 到 九 个 等 式 ， 取 出 共 中 两 售 等 式 
-a = -eoseid + sindeosd — py ) 
P > yaingeosð + dusel — ege fett H + psingsinð - y: 
Hi- JC SI) 
-a7 - sinfi- 0) (d-t) 
LI a=- E (3-63) 
有 第 二 式 ， 和 有 
P = psin(d - 0) + peost- 0)- yp: 
= Yaind prosa — Ps 
— Got, 
D 4, (3-64) 
6-3 质心 动力 学 方程 


现 选取 固定 地 田 坐标 条 〇 sxsyszs 作 为 质心 动力 学 方程 的 央 考 化 标 系 。 它 相对 村 地 球 是 
Rm. Wi, Lë (8-16) 成 为 


Vag Frot Sasso heen 

1 a ar 

A i | KÉ | haet Saeys t de nn {3-65} 
Vas Faot Saers toss 





再 次 指出 ， 这 里 vE Ki a D AHER PARSER o 
疯 先 考虑 作用 力 F 的 分 量 
Q 


Fas Py: - Grs 
r ] Pye | el Y OI Gye | - (3-66) 
Eao Pa Z ds 
HPP, Pyr Pa ARRDILE AAHAS AS nPE, ， 游 动 小 发 动 机) 产生 的 力 的 分 
Si -0, Y, ZATEA RRADHE Jaos Gen, He BIER. 
H FEH mE gA a AER o. EJLA R Ro 不 难得 到 它 在 地 面 坐 标 邓 中 
EW 





X 
gs5 -g 于- 
I 
Ys trp. ` 
Gëff? r (3-67) 
z 
Ger D 





Drzüäeatrairi? 
共 中 ra 为 地 球 平均 华 径 


CN sit (yo tro)?+ af (3-68) 
为 飞行 器 质心 至 地 心 的 距离 。 
EEE D 
S.=-mõopx(@pxrj= -milost Op- eär] (3-69) 


D 57 D 





SITE Q cosy p SH 
m | CAD an ! =B? 0 | sinP eg I 70) 


DTN -wp — eos p cosp 
r.g Xg 
BIER 
Pry Co 
Fik p.l = wpa 
其 中 a Xocogosing + (ys Trplsind — zaeosiacosit 13-71) 
ER TRR r ÆR E Lii. JEA LAFAR (3-69) 中 ， 得 到 
S perg - Fang apiiii + Xy 
| KE "mon | - Fainga — Yg En | 13-72) 
Saras Fotos acos u + Sa 


Caiolis WEDS. = -2mp xvj iA 


Sarrg COs SiN aay Usg 
| KEE 1 -ameo| singa | | Uyg | 
S 


trzy 一 rosi orostt Jan 


= eoa E peos HU pg T ENPU 
= Pmp | Goal ups ln, o COSE oëin Dt ey (3-73) 
Sin ul, g ~ eP psiniivyg 


wR. 把 式 (3-60), (3-61), (3-62), (3-63), (3-64), (3-66), (3-67), (3-68), 
(3-72), (3-73) RAR (3-65》 中 ， 持 年 若 虚 到 气动 力 Y 、2Z 和 外 可 以 表示 成 


Y =cjagsS 
Z=cf fgs = -cepys | EE 
GËT EN 
Ahe, Lë Ya Fan P yA Eh EARL 二 是 得 到 质心 动力 学 方程 如 下 : 
d ze Ab - 2 
m Reie := Purcos- Psing— Ce heel -ceogS t- 8) sia — mäer Gi 十 
LZmont— e2aiemsiiitwn —ëlogian, g) + Mohl- rocosmosintg Lei 
Dou 十 
m der -= P. atnf t D gef cgs sind +eggS ip- Hoost -mgr Sege +! (3-73) 
+Žmopl paos HUV gg t eosPasinuY og) A MORC- rosino Ziäs Frp) 
dn, , z 
m= -Peat rhs be Ob, ER D Wl mëi er ai + 


Zmeonntsieg att, copsintvya) 十 MOD oco pocosH + zr) 
其 由 的 谈 最 rf。 用 式 (3-71) 描述 。 
6-4 转动 动力 学 方程 


ARAA AREARE Amta S (以 后 用 7 RR) o m WER 
=Í, a= 0o IEAk, AR Kirn E REGLE A mie E EERE e, DEL 
转动 动力 学 方程 应 从 第 一 章 4-2 中 式 《1-42) HIE. HIRA (1-40) o ZEL Sa, F 


D 38 D 


动 动 力学 方程 成 为 


duu, JI, 





Ta g vos, 
dey dL , 一 
dr T WÉI ety T Wi oro: H, 3 ?6) 
: li 


Hp EF UAEM GE e E DEM o ERM DAN ep3 

WEDA EE DEM o, AAEM a WRA NEM p: 

ISät IM, ARE EM op, RADA- M ERRI DENM ppo 

为 了 分 析 的 便利 ， 把 - 之 .ov 称 为 深 转 内 阻尼 力 焊 ， 写 成 Mo 把 - Ta, 称 Ai 
ARLE JIE SEM 2o E- Teo, EA IR AR E DE, GA, LAS] 1114 
Lol, Apap aH RE Latr NE A RRE R H TAERA TE TE EE A SAE 
Prät OS ët nt E E Dei, A mg EERI Corielt M D SG, HALAH TA M 的 计算 公 
式 。 BAE JHS, MERA A E pE e A EE A ARL E A Ho 

于 是 方程 (3-276) 可 以 改写 如 下 : 


7 = (WIH M? ost Map | 

, | 
j dor =Q -i Down Mot Mm M -ptM ON (3-77) 
T EEN -os 下 星人 中 十 上 < 人 


6-5 运动 学 方程 
质心 的 运动 学 瞩 程 可 以 直接 列 出 如 下 : 


dxs 
dt "eg 
dys 
mn "Dan {3-78) 
dze 
di = Vau 
v= vs +u, tuže 
CEET EEN | (3-79) 


Lef = VaVe 
H DEDZA AH GH re, MERIA F E ARH AA 度 四 ,-* 由 
mëi 【插图 3-13) Wå, e Eër paR phs pR. MOLDA A hy, 的 #cosy 和 
"EN, i- psinys Ao Eia RhA 


a 59 D 


CO rt í 1 0 了 -på 
EEM jee HE 


求解 卡尺 ， 得 到 





d? S= Warea CP patse 
d | ， 
WW 有 人 【人 


A Sorg at | 
转动 动力 学 方程 (3-77》 ho ho ES) E KERE) 绝 An, ieai PIR 
(3-80) 中 的 四 ，* 必 飞 行 器 相对 于 地 球 的 角速度 。 它 们 之 问 的 关系 是 
"Go" D: 


所 以 在 机 体 坐 标 系 中 的 分 董 的 关系 为 


SET EE Gi WEIER 
SÉ D o, Lei on. (3-8 
Dr CD = CD 
共 中 的 wp.zs 等 由 式 (3-70) 给 出 。 
当 研 究 转 动 动力 学 特性 时 ， 有 有 时 可 应 急 略 地 球 旋转 效应 ， 内 而 可 以 取 


{ 一 如 x 
| DIE EEN "| y | (3-82) 


EELER DJs 


6-6 关于 操纵 力 和 操纵 力矩 的 补充 说 明 


ah DS kën hp... pe LR DRAAM, Hun, Mapte HIR MAET EIR J 
GEREENT ed ao EISE 
KIRA. BRIER A E HEA PRAT o 

HARARE i aiH 1R. AE 2 A A At, ie RRETH: REI nn." 
m EE A RETEA = Pela i Pe J EEI 3 E " 

He az, Saäomtbr, wë sai RI 
转角 为 
ô, -0 一 六 | 


5, =- (8, Hia) 


” | (3-83) z 
, 1 
= 2 (6-6)) | 
zS, 84, A FREA Eo ô Mó; AEA Mte 
为 正 。 
它们 产生 的 力 可 以 表示 为 图 3 35 WARNAE OH) 


+ ô D 








ye 72A: 
P=- Ali + 6s) | 
其 中 4 表示 艇 偏转 单位 角度 时 产生 的 法 向 力 。 
ët ei ENKE 
M p=- LEER \ 
Mopo Alel, hòa) | 
Ms Ae, (ës 4,) 
其 中 7 是 堆 舵 面 中 心 到 重心 的 纵 问 距离 ; 
ec 是 从 能 面 中 心 到 纵 辅 的 距离 。 


6-7 小 结 


本 省 中 葵 谍 到 典型 弹道 式 火 箭 的 主动 段 的 固有 特点 ， 


{3-84} 


(3-85) 


AAR RRA ET T 


BIES A F” BAA (371), (8-753), (8-77), (3-78), (8-79), (3-80), (3-81), 
(3-84), (3-85) o MAAE ZABE MA F Ma ARAS ATR EA Ak 前 的 主动 段 的 
BaH, HRH: E ATE AS Sth t- AEE KITAR m EAE ah A RIR 


不 相同 的 。 


在 实际 进行 暗道 优 火 箭 的 主动 段 的 计算 和 分 析 了 时 ， 还 可 以 极 据 不 同 阶 段 的 具体 策 件 .把 
方程 作 不 疗程 度 的 简化 。 关 于 这 个 问 巅 ， 在 交 献 [10] 中 有 间 细 的 论述 。 


= 6l e 


第 四 章 改行 器 的 小 扰动 线 化 方程 


d :和 第 三 人 章 建立 的 运动 方 称 烙 为 公 而 地 描述 了 了 飞行 器 的 运动 。 只 要 其 备 是 肯 完 会 的 你 
气动 力 狼 据 、 发 动机 数据 以 及 控制 规律 ， 求 解 邢 些 方程 就 能 得 到 飞行 器 的 延 动 过 程 。 包 括 飞 
行 问 质 心 咯 移动 和 绕 质 心 的 转动 。 内 而 那些 方 释 可 以 用 末 不 仅 研 究 飞 行 航 水 问 是 ti HI 
稳定 尾 和 柑 级 性 问题 。 但 足 那 些 万 各 都 吓 变 系数 、 非 线 必 的 方程， 绝 大 多 数 情 况 下 行 不 测 解 
折 解 ， 共 能 出 数值 法 求解 ， 即 在 县 体 的 数值 条 件 小 得到 具体 的 数值 结果 。 这 样 ， 误 不 容易 三 
纳 出 带 有 普 源 意 之 的 - - 般 性 规律 。 

为 了 便于 研究 飞行 器 的 稳定 性 和 操纵 性 ， 通 常 要求 把 运动 方 称 化 成 常 和 水 数 线性 方 称 组 
《包括 微分 方程 和 找 烙 方程 ) ， 这 样 就 可 能 用 解析 法 求解 或 进行 解析 噬 究 ， 帮 及 居中 是 纳 出 
一 些 普 沾 规律 ,提出 一 些 飞 行 号 质 指标 ， 作 为 飞行 器 设计 的 指南 。 

KAEI EURES A EEE, FR Rie A RR A AOE AE aE A 





论 。 


Si 革 本 概念 

1-1 基准 运动 和 有 扰 运 动 

KIr ERMAT WERA TE tbe HA E TEGE a EOS EE A 
IK aha ch, 基 淮 运动 做 数 用 MER * FERI 例如 Pa ago go 

FRAIES RERE Aë Sit 2 RAAE a ARES a E e dra e R ti 
CHAVE) HEERO W o piE MIRS BN AR m EARR M A ESTER oE A 
Së, PH mi i nc, a: 8. 

DiR ER RRT AEE a E A Dim, BHEE a a aE ah n ax 
有 一 定 程 度 交 大 为 因素 。 但 在 天 多 数 情 况 下 ， 达 样 的 划分 是 自然 的 ， 公 这 FERI 

E EE EE EES EE EE ECKE EG EE P a i EE 
示 。 例 如 

Jusen ves AJozsog-uae d= - Has {i 1) 

丁点 要 建立 的 方程 起 在 一 定 的 条 件 下 得 到 的 忆 微 小 偏 量 为 对 象 的 线性 化 方程 ， 所 以 称 为 
小 皖 动 线性 化 态 程 。 与 此 相 尘 地 ， 把 第 二 、 第 -: 章 建立 的 未 径 线 化 处 理 的 原始 方程 称 为 全 时 
方 称 。 

1-2 SS 

TERREA PEPHER ERKKA E FA E E ATE KAT e e ii A FE 


乎 以 线性 化 。 双 于 考虑 地 球 明 率 和 旋转 时 飞行 器 的 和 运动 方程 的 线性 化 ， 在 必要 时 读者 可 以 参 
照 本 剖 的 原理 和 方 活 来 完成 ; 在 非 平 静 大 气 中 飞行 器 的 运动 方程 的 线性 化 问题 将 在 第 撕 章 中 


e 2 > 


齐 述 。 内 此 ， 我 们 已 经 降 含 地 采用 了 下 列 盆 说。 

1. AARE EPRA, Een, 

2. nët tram, ARENA KTA kirëte hon, vr Dm 
气 摔 动 所 产生 的 力 和 力矩 作为 输入 量 加 到 运动 方程 中 来 (参看 党 五 章 84) 。 

3. 试 为 飞行 器 是 刚体 ， 不 考虑 结 移 弹性 变形 、 操 上 然 面 特 动 的 惯性 影响 以 及 转动 邮件 的 
陀螺 效应 。 

4. AHELA., Ria ac A B AERE Oxy, 对 称 的 . 

此 外 ， 还 作 如 下 三 点 柱 庞 假设 。 

5. 在 基准 (未 受 扰 ) 运动 中 飞行 器 的 纵 对 称 平 击 处 于 铅 重 状态 mA KIE ERE 
ALA Eer, Im 0, af, genat, LIRE ZER II ES ES DI Mr EE k ein 
BEAPHAR TEE 

6. 假说 基准 运动 是 定常 的 ， 即 其 淮 和 运动 邦 数 不 随 村 间 而 变化 。 实 际 王 ， 只 要 基准 运动 
参数 的 变化 不 太 画 烈 ， 则 在 所 和 研究 的 扰动 运动 的 短 时 间 内 就 可 以 近似 地 认为 基准 运动 是 定常 
的 。 这 样 处 理 法 称 为 系数 你 结 法 。 

7. 假设 在 扰动 运动 中 运动 套数 的 仿 基 《例如 4v、A&a、48) RS, Dëss 
的 一 醒 及 更 高 阶 的 最 是 可 忽略 的 。 这 就 是 通常 所 道 的 小 近 动 很 设 。 

一 般 说 水 ， 在 小 扰动 假设 下 就 能 把 飞行 器 运动 方程 吉 行 线性 化 。 再 如 上 上 上 述 其 它 假 谨 ， 
就 可 使 线性 化 运动 方程 分 离 为 纵向 和 横 侧 两 给 彼此 犯 立 的 方程 ， 并 且 是 常 系数 的 线性 方程 
( 详 见 83) 。 介 是， 如果 雌 兹 运动 不 是 定常 的 、 厕 是非 定 常 运动 ， 则 线性 化 方 桥 吝 是 变 系 数 
的 方程 。 

对 飞行 颖 运动 方程 的 线性 化 尽管 作 了 这 么 多 限制 ， 但 实践 表明 ， 线 性 化 运动 方程 的 适 骨 
范围 仍然 是 相当 广阔 的 。 至 于 这 些 条 件 完 部 可以 放宽 到 什么 程度 ， 例 如 ， 运 动 套数 偏 野 可 以 
大 到 什么 程 论 ， 这 个 问题 很 难 作出 明确 的 答复 ， 这 只 能 由 实践 来 检验 。 在 理论 研究 阶段 ， 可 
以 通过 把 允 最 方程 的 计算 结果 与 线性 化 方程 的 计算 结果 进行 比较 ， 来 制定 线性 化 方 穆 的 适用 
性 。 

当然 ， 也 有 一 些 飞行 情况 是 不 能 用 线性 化 万 程 来 解 完 的 ， 例 如 飞机 的 失速 后 尾 旋 《 这 时 
扰动 是 相当 天 的 ) ， 用 回 体 火箭 发 动机 进行 加 速 的 飞行 阶段 (这 时 基 淮 运动 套数 变化 大 剧 
A) 。 这 些 飞行 情况 还 是 最 好 用 人 孚 其 方程 进行 硼 究 ， 不 要 损 强 条 用 线性 化 的 方法 。 





1-3 方程 线性 化 的 一 般 方法 
飞行 老 的 任何 一 个 运动 方程 〈 微 分 的 或 代数 的 ) 可 以 表示 成 如 下 的 一 般 形 式 ， 
Ze, y, o w)=0 (4-2) 
其 中 的 x， yp ovo weba gik etib a. 
RERI 2ahsëu URT e Ea -5 EAN 
zz pat Áx: Mrs äech Au to w=wyt Aw 
于 十 方程 (4~2) 写成 
fxnt Axo yst An: o wat Aw)=0 
ERER (ego Wx， …， we 处 展开 成 Taylor 级 数 ， 厦 且 忽 略 二 阶 基 更 贾 阶 小 县 ， 则 
写成 








_ f 
Í xy Hgt, Wajo (a ) Ax (下 Annne ) Aw=0 WEEN 
Efast MaA FE (4-2), Bu 
Fixar Han tto wa) =0 (4-4) 
E (43) RER (4-4) ， 得 到 


这 是 任何 一 个 运动 方程 (4-2) MAREEA Rain SEIT DU BE E, 


$2 力 和 力矩 的 个 量 和 导数 

2-1 力 和 力矩 的 偏 量 的 分 析 

iah A BRRR. BAA DADER T BAA DHA hl E HERE 
AAEE, REH ER EARE. 

设 此 个 力 或 力 夺 4 是 若干 个 变量 的 丙 数 

A= Ala, bs pw) 

其 中 的 a bn, 与 方程 (4-2) 中 的 xo ys ee’ ww 可 能 有 一 部 分 是 重合 的 。 SEI, A 
的 全 时 应 为 








CR! | 2A L ana d fin 
A A ER Aa SI? Ah + 3g Ag { 1—6) 
或 写成 
JA A7 dat Ah Abt A dn ERC 
其 中 的 俩 导数 基 在 基准 点 《xs e’ wa) HDD, HN 
可 一 2A dl ven 
过 =( ES 1 " Zo Lon: bgo s Ga) 


对 于 作用 在 改行 器 上 的 力 和 态 征 来 赔 ， 所 到 的 基 崔 点 〔 蜜 基点) MAE KRA. H 
米 扰动 运动 状态 。 和 前 已 指出 ， 把 所 劳 虑 的 基准 飞行 状态 限定 为 钻 匡 平面 内 的 对 称 定常 占线 飞 
行 .要 据 宏 气 动力 学 理论 和 实验 数据 的 分 析 ， 井 结合 有 关 的 物理 概念 ， 可 以 知道 ， 在 这 样 的 
基准 点 的 邻 域内 ， 作 内 的 力 和 力矩 到 类 于 吐 些 内 数 。 
发 动机 推力 已 取 并 了 和 速度 V, mE AARIA Eâ, A 
D Di, hs Gul 
Bu CIS 2) 
AP=P" Avt P Akt P’? d, 
但 是 由 于 发 动机 系统 《特别 是 洞 纶 奔 气 发 动机 )》 有 相当 大 的 惯性 ， 从 油门 村 位移 到 惟 力 变化 
之 间 和 有 相当 大 的 时 间 延 退 ， 所 以 用 忆 2> 8， 表示 古 不 通 必 实 的 。 最 好 用 可 榨 制 的 推力 偶 量 
Ar, RRE TE 
AP=P dat AA AP, (4-7) 
气动 阻力 包 通常 取决 于 、 A OB, MEAE. AEAEE d W, our REN 


sa d >” 





"tt, ën (LEI ite) OQI) a150. AE 
AQS" dv t Q Ah+Q Ae tOiz Aë, (4-8) 
其 中 的 Q 项 有 时 也 可 以 忽略 。 
升力 了 主要 取决 十 v0、 和，&， 其 次 也 随 5: 而 变 ， 轩 严格 地 说 ， 还 受 迎 角 变 化 它 a SO 
FERE o s DI SC, BLUE 
AH =F" AvtYtAA+Y daty’ Aó tY Aa tY Aan, (4-9) 
在 大 多 数 情况 下 可 以 忽略 Y“ 和 了"”*， 在 粗略 分 析 电 全 时 还 可 以 忽略 了 *。 
ioo, RI, EHME AE tik, 和 但 在 扰动 运动 范围 内 这 个 e 
小 ， 不 必 荐 上 谍 ， 因 此 可 以 认为 
Atmg)=0 (4-10) 
mu. wm, A, a, å, oe, dë, WIE 
AA, e Mr Aa At ARAME Aa +MEAGtME Aost Mi AA, (4-11) 
HPA KHR Aim LM Aio 
气动 便 力 Z ARTo, A, MEAS, REARMARE o., o, ARAARA A, 
对 此 本 来 应 该 写 
AZ= Z Av+tZ* dát AB+ Z° Aost Z dot Z Aò, 
BITERI, war 0p, Ós, FED EZIEH TRS REAGE Ao, AREE GER 
Ee AKARA ARED o ARA 5H 
AZ=Z AB +Z dos +Z™ jw, + Z A8, (4-12) 
E S oii Fe TEAR ZMZ éi 
BR DAM mt hurt, h A FERRETA, or, oy, SARE 
Hafi. HEER ERKERE RAE v AAR Su 
Me= MÉ AL MI Aost M Ao, Mi dôs +M Aò, (4-13) 
AM,= MS APBtMe Aen, +M Aen, tMi Ai + Mes do, (4-14) 


2-2 力 和 力矩 的 导数 

现在 求 讨论 在 上 一 小 节 中 周 现 的 为 和 力 短 对 运动 究 数 的 偏 导 数 ， 例 如 Pr" Ver, MEH, 
JI, EETA Däreg. 

第 -组 推力 的 导数 

发 动机 推力 对 速度 和 高 度 的 导数 P* 和 2?*， 可 利 楼 模 据 发 动机 特性 曲线 (rz, A) 用 数值 
微分 法 在 基 淮 点 Co, ha 人 外 求 出 ， 不 要 再 化 成 其 它 派生 形式 。 

第 一 组 ”气动 力 和 力矩 对 速度 ”的 导数 

WEG: Dogisn, "ab DOS 


o=o, Laun 
D ET pu, MEMA Marh 数 影 响 c，{ 至 于 通过 Reynolds 数 对 cx 的 


影响 则 可 以 忽略 ) 。 于 是 


e 65 e 





f 1 
On SC . EN 2 PLAS Gei tiatie hb 


TA 0 io, WOM :Do 一 1yas。 用 cy 表示 3cy (244, Ek ëlo oan Mia, ME 


ek 29 = OU - 了 CHAf LC 
RI 2 


(4-1) 
同样 地 得 到 

Y 、 

YY "一 Sé = Pyg (-; cy Mut ess ) {4-156} 

KEE oMa = Patyo ba (-ż = mM Ms TMg ) {A lae U 


EE TEE EECH 

第 三 组 TADAN E km 

IO Dm, Ami, REER EA TEE o "Hä SE a 的 变化 影响 好 数 进 
TRWA D Aceso WIH 


‘x oM 1 ,9 1 











AFM =vjas kk 
oM ww Ais 
ða 和 ER (a 
A. 29 
3 ， ^ 
局 e = 2- DH cH y 了 十 deeg (4-1€60) 
同样 地 得 到 
EY 1 H St 
Y'= E pavå zl. ee TT Ceä A 了 下 六 
2M. 1 、 oi oi 
Mi= ay = 3 -pav Sba = em H Vë. a Myg Pa ) 《41-16c ) 
这 里 的 梯度 ae* 和 2 可 按 标 蕉 天气 的 公式 导出 。 已 经 知道 ， Wika HAR 
` ( D 4.2558 
1,2250 (EEGEN Ae 11000 
p= p11 [kgy ms | Lmi 
A- 11000 
0.36391exp ( - ajre ) 11000 <4<20000 
=f 20.047 M288.15- 0.00654 `  _  A<11000 d 
295.07 ”000<A<20000 
肝 此 可 以 求 出 
po — 4.2559 
E en h psi1000 
-5 一 = | [l/m] Cm] 
一 1 1 Lü be: 20000 


其 中 


0.0065 




















Əh -| 2(288.15 一 0.00656) Flya] 
a 0 
第 四 组 "ai: aime CSS 
这 些 导 数 可 以 徊 接 写 出 如 下 : 
Q7 -3 Garg 
Qi: = Cd 3 Garë? 
re= IY copd pucks 
ki ZS Sek va 
KETTER 
HS ECH = må? > Pavi Shy 
Z= SCH -8 2 progS 
Zr = > Pauko 
Mi= Er 一 ms > PS 
Mis 3 = mi 3 DT Si 
Mi EC = mir- È pots 
az Sat = mf R et Si 
vie än 一 min L DautSi 
Manz Séi mg” 2 pyu S! 


CH Sged o 
第 五 组 ”纵向 气动 力 和 力矩 对 角速度 的 导数 〈 纵 向 动 导数 ) 


首先 看 了 


Ae UIO00 i 
[m] 
11000420000 


id-17) 


(4-184) 


(4-1865) 


{1-18c) 


{4-180) 


a 


4-l8e) 


Ba) 


4-18 f} 





一 


4-189) 


í 4-184} 


(4-18; } 


{4-187} 


(4-184) 


{4-181} 


(4-i8m) 


(4-187) 
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引入 无 系 纲 变化 率 





T= at 
Da 
且 引 者 导数 符号 
EIS 
"7 "a 
(NS 
dk 1 , 
Ya 一 Ee Sch? Dain ANDA 
WE, A 
oz oY a 1 
NS Be, 一 eg CAS 
Mi= D mit pavas bà 
o2_ IMa ə 1 
EE EE 
其 中 wi ， mis ÊT., 
Ce 一 Zo o. = fg 


第 六 组 ” 米 侧 气动 万 和 用 御 对 角速度 的 导数 〔 横 便 动 导数 ) 





o Osl ` Oyl. 
BT 2o SEH 
并 共 引 入 导数 符号 
ov Om Dn — dm, ` 
Me "ae, mx T Ae, 
me 二 om, my = ge 
wr. Oc wy _ _9cs 
"e "Ae, ` z — J0, 
于 是 横 侧 动 导 数 可 表示 为 
ZL"*= Z = 021- puvaS! 
A Tp- 22 =e% L PUS 
H 


C4-10.:) 


(4-195) 


{4-19} 


(41-19.7) 


{4-204} 


(4-20b) 





Me = DT E PDU Si’ {4-20c) 
Mo = Be me i Pava Si? {14- 20d) 
Mi = mo l eiis Sl” (4-20e) 
Wie = m” 了 Data Si? (4—20f} 


以 上 仅 列 出 了 气动 力 和 力 第 的 导数 的 夫 达 式 ， 而 革 在 表达 式 中 叉 引 入 了 气动 轧 采 数 和 气 
动力 矩 际 数 的 导数 (例如 cs、mmx2*) ,这些 导 数 通常 称 为 气动 导数 。 本 书 不 打算 进一步 研讨 
这 些 气 动 导数 。 气 动 导数 的 研究 足 飞行 器 空气 动力 学 的 重要 课题 。 通 常 包含 以 下 卡 要 内 容 

1. 各 气动 导数 的 物理 本 质 ， 产 生机 理 ; 

2. 各 气动 导数 的 工程 估算 方法 和 实验 广汉 ; 

3. 各 气动 导数 与 改行 痪 几何 外 形 、 部 件 参 数 的 美 系 ， 以 及 随 邮 数 的 变化 - 


8 3 和 运动 方程 的 线性 化 
3-1 动力 学 方程 的 线性 化 
质心 称 动 的 动力 学 方程 { 虐 称 力 上 方程 】》 用 川 航 迹 坐标 系 中 的 方程 组 (2-29) 的 类 项 个 万 
程 和 机 体 坐 标 承 中 的 方程 组 (2-31) 的 第 三 个 方程 。 
首先 进行 式 (2-29) ER ten 
m o = Peos{ a + Ppop — OU 一 ngsing 
的 线性 化 。 按 一 般 公 式 Bag 5) WE, "ëmt H 常数， 有 
Ka Gef 一 d'Soen? géet pter Lee Änne, - 
~P eos HEET AQ- vg Gene, EE) 
由 了 寺 基 准 状 态 的 2x 一 0， 和 且 把 4 忆 和 4 的 表 太 式 (1-7) 和 (4-8) RA. "aa 


m ddr =(P" Aut P* Ab: A Peosay- pp) Pasin ay typ ju- 





- (Q" Av+ Q Akt Q dat QÎ dé.) -mgood d” {4-21} 
HEAR A EACL- RAA A EREE, A EA EECHELEN ETA 
az, Zoch: 样 及 把 4PP 和 4 了 的 表达 式 代 入 ， 最 后 得 到 


ta" 147 =(P" Aut P’ Aht APedainl rgt Pp) t Pycorl ryt Pp Auch 
+Y” Avt PedAh+ Y AatY daty” Ae, Y’ gS) + 
+ ap Geif AË {4-22} 
把 式 (2-31) 的 第 三 个 方程 线性 化 ， 
= 69 ， 





f : 
人 Zi “十 Aach, ak om AJ Du: 一 AJextieea yw d Ei 
-= — A sin Dy — Late xd fe Ka A Zos- Z anin fi A É = 
一 mnte Jinny ap Ausl asya A Y 
HF oss wyg batya an0, 于 是 
d gege KREEG or Vy )= -Qy ABT AAL mm Gross Ay 
H (2-55) afg 
Av =u A P Vera ug (Ug TU) 
Dyrt Urep p " Ugi 
HERH A Zm ELE (4-12) ， 得 到 
MU AAE 一 Uag A Oy "iis doe 
= -Qg AU ZEAR +Z Aen, tZ” do, t Z’ Að, + 
-mocopBe dd" (4-23) 
REDRAW DEAE ORADE 卢 程 ) 的 原 姬 形式 是 式 (2740) o JEH iE 
包括 气动 力 夫 和 推力 征 , 颁 进 性 动机 安装 是 对 称 的 , Git H Aa ng MACH), 
于 是 
SH, - Ma D SM, =M, b = Gë =M, 一 Ze, 
Schu. AM. Ma tcahhä, fe, KE ahoi, 
"A4 (2-40) 的 第 一 个 万 程 线性 化 














了， "Ac Te KAya Fg Ao) 
“Fay (Aowa t Osud tr — dder) AM, 
H Fos séien un, ERIE AM AIREA (4-13) ， 最 后 得 到 . 
d Ax d Awy Uy EI 
Weer ckt femme =MEAP M Ae, M Awt 
+M AS +M dò, (4-24) 
同样 地 进行 式 (2-410) 的 第 二 和 第 三 个 方程 的 线性 化 ,结果 是 
dA do, o DI 
I, SEPE La E ME A e Me den, tM Awt 
Mis AS, +M? AS, (4-25) 
p, ŻA Ad MAR MS Ao ME A Me Ao, +M? AA 
- (P Avt P*Ak+ dP.)e, (4-26) 


3-2 运动 学 方程 和 几何 关系 式 的 线性 化 
首先 把 运动 学 方程 (2-43) 的 第 一 式 线性 化 
* 70 = 








d 
d xn. = ege fl geosp es — agin fa A Begeff, a — OmeosO sing, A Ya 
HOF #50, BEI 
idx . 
Zä SEH DE EE (4-27) 


MIE RR AA RRR EE. ATHERE dhidi RTS duo dta 结果 


p 


ddh Sinb a A+ uscod y d 14-28) 
ddz 
` o ?=- serge Ha dp. (4-29) 
转动 运动 学 方程 (2-46) 第 Amie 
d 
- Ar = dos- Ated(e, weosps — Ossinpy) - 


-ted#( A ycosp ~ Owi Sin Ha Au: A sinra™ ww 05 a du) 


HF Ory T Oyy TD YaTT 所 以 得 到 








aay = jo,- igs Ao, (1-30) 
同样 地 ， 第 三 和 第 三 式 的 线性 化 结果 是 
TK A 
dt costy 4 SG C1-31) 
dA ` 8 
"ur "d: (142) 


ILME (2-48) 的 线性 化 
vos ls A P = [vosy A Ysindgeoslp — ps Ja tsin yacosha do — Ps) 
— sinyasingysin(p gé, lat A9 A, sintya d ysinl y — yr)x™ 
+eosyxcos p- Pdl AP- Ap, lege — Lin Masina ý — p, yt 
+ us apiui d — Ja ein fl AO- rosya peosdwsind s+ 





1 siu yasinta d'Zain y ~ sin Lagos wcos d y A D 
B T Gaza: PaVa 0 M 
AB = lsinda AYT AP- A oe — wossindr d Y 
: (Ap Aë, ro tsina- Hei AY 
由 于 zz - 8 二 vw， 最 后 得 到 


AË Seol Ay- AY) tinag dY (4-33) 

类 位 是 ， 式 《2-49) 各 式 (2-51) 的 线性 化 结果 是 
da= åf- A8 (4—34) 
ys=t80#d P + Za AY (4-35) 


5-3 纵向 方程 和 横向 方程 的 划分 


sit, 已 经 得 到 了 所 需要 的 翁 部 线性 化 方程 ， 即 不 《4 -21) AA G35). ARE - 
Lann, 这些 方程 叮 以 划分 成 互相 独 站 的 两 组 。 

六 各 (4-21)、 (4-22), (4-26), (4-27), (d-28), (1-32), (4-30) M um nl Le 
at Wio, Al, Ah, Aen, Aa, Aoa A8, AP, Aao AE IE H ER AA hii 2 
IERI Pie 

diie (4-23), (4-24), (4-25), (1-29), (d-30), (4-31), (4733), (1-85) IUEL 
BEREDTE E AE, Ady, Aos, Ap AP, A Pe, Azp Ai, An, du, 这 些 
JE Ee Be A R E aa a A Ea 

J, EEGA EARE E KEA fe Ko bE anes A Fen NAs 
ROHT MAR ARARA MA EH mifi AE A ARE - 
Bri A Ea R EE A TEHE. og A MH A AEE A gR A LRR E EIERS 
SERE A FEM ARRI — Ae Ar E ITa tR To REME. 

MERS A A HITER NETSA aia h A A ME ai A FE o 


Su ”纵向 小 扰动 方程 组 
4-1 纵向 小 扰动 方程 组 的 自然 形式 


『 征 于 分析， 现在 把 纵向 小 提 动 方程 组 的 各 方程 进行 必要 的 加 工整 开 。 把 彼 下 正和 相 影 
y OLARA) 的 方程 收集 到 一 起 ; 葵 利 用 关系 式 40= dé- da 和 4 =d 8dt， 滑 去 
ERAMA 把 49。 和 ,1P。 看 作 输 入 变 最 ， 放 在 方程 的 右 中 ;把 4v，4Ac、.1#、4# 
看 作 和 输出 变 央 ， 上 放 在 方程 的 左边 ， 且 拨 厅 序 排列 。 于 古 得 到 焕 合 方程 组 如 下 : 


D 
d Sidd [- Dcostant TP) tA] J Au +C Pasina t Ppl +Q- 


- a one y] Aa t mgeosby AT 一 Poslat ga) rQ'JAR 
= Q AA. 3 stoe tpp AP, 


. d 
-[P"sinl agt pp) +Y"]Av- (ue: 十 mpelze 二 [了 Dags ag rE) 一 





d -| e 
- mgsint e] lAa 一 [er 一 mes) p H mgsinbs | 一 (1-36) 


"Das as tpp TY AR Y AS, tintas i Pp) AP, 

lfa d Ifa - d` KÉ d 

-CM pre,) Au (7. +rMs)dart (T. Ze MI 
pe 
dAn 
dt 





Jar- 


-siny dn Ugros g Aa- vecos a AA+ -0 


Krie AC., Aos, Axpo SETER A JT EREK GUAR UN E A AE a A H 
ev 72 >» 





A0= Ad- do 
d A 
Jon, zs 
di {4-37) 
xD 


di = rosy AV- Gesin fl d 


HPR KRR ht EIN, EAHA NH (4-36)。 当 给 定 40- 和 4P。 
由 ， 利 用 此 方程 组 可 以 得 到 运动 变量 4un，4der、dzg、4 的 动态 响应 ; 由 此 方程 组 可 以 扒 - 
导 思 飞行 器 的 传递 画 数 ; SS 46,=0，4P。=0， 则 此 方程 组 撒 述 飞行 器 在 操 杆 状态 下 的 
纵向 自由 挑动 运动 ， 从 中 可 以 研究 飞行 器 运动 的 动态 稳定 性 。 

线性 化 的 方程 组 《4-36) 仍 保 持 其 原来 的 【线性 化 以 前 的 ) WEEL. Tagih 
上 方程; 第 二 个 是 法 向 万 方程 ! 第 一 个 是 依 仰 力矩 方程 ; 第 四 个 是 高 座 变 化 方程 。 因 此 ， 称 它 
为 自然 形式 的 方程 组 。 

为 了 今后 研究 和 计算 的 方便 ， 引 及 如 下 的 纵向 动力 学 导数 : 

x” =| D'Senstäs tpp) - Him 
x” =[ - Pasin ay tpp) Q7] m 
zech zs D'ieost ay tE) Qim 
xte -Q jm 

和 2 一 cos[ aat Pp)/ m 

vz P sinf agt Pp) tY]? mog 
y [Y+ Pacos ag + pp) li mvs 
y= H €i mug 

y =Y”) moy 

yt =[P*sinlaa tpp) +Y *] mus , (14-38) 
vis Kei mog 

y” 一 sinf agt Pp) i mis 

us=(Mr- Diesiil- 

ns= Ms, 

EKOS 

u SME d, 

ut (AE Pep) 

pte Mejl, 

HR= -epii: ， 

这 里 必须 指出 用 点 。 第 .…， 这 些 导数 仍 是 有 量 岗 的 ， 没 进行 无 县 岗 化 ; 第 二 ， 这 些 导数 
仅 伐 者 推力 和 气动 力 ， 汽 代表 重力 和 愉 性 力 ; 第 三 ， 在 这 些 导数 的 表达 式 中 出 现 的 气动 力 和 
mm kB (Pmr, Ai, MI) MRAR (4-15) Æ (4-20) 计算 ; 第 四 ， 每 个 导 
数 都 有 骨 显 的 物理 意义 。 

利用 以 上 定义 的 动力 学 导数 ， 并 且 引 用 微分 牌子 p=df dt 则 厅 合 的 纵向 小 扰动 万 程 
(4-36) 可 以 写成 


D 73 D 











ip- stadu- ixt tychy) Aat goba dt- xt AAT t A t AP 
= Ae Juni em -- 4 auf, | Aa- D (t-i H sinda jdd 
Za Da 


— DO 
-yt A= AA + yr AP, dek 


SIE DEET EE EENEG EE TEE 
= sia y dv 一 Varot) 3 d E Varo y A i+ Pd h=0 


4-2 ADRA RNEER 
H ERTE R EA d IT E B a T ERR ES ie E, 
ax =Åx+ Bu ii- t0) 
其 由 vleit Gas 
AEn s nE RRI RE Ok PRIRA HRE) ， 
B 是 ”xx 二 维 控制 定 阵 ; 
u Aarb, 
这 种 标 玲 短 阵 形 式 的 状态 方程 庆 许 多 卉 成 有 向 成 熟 的 分 析 方 法 和 计算 方法 ， 侧 且 特别 便于 
利用 于 入 机 求解 。 所 以 希望 把 飞行 器 小 扰动 方程 组 的 自然 形式 转化 成 标 竣 仿 阵 形 式 。 
Du, Tram AtIti 291, METR SE Io, Ja, Jee, A8, dk: FE 
经 过 简单 的 代数 变换 (由 第 二 个 方程 得 到 pde 的 表达 式 ， 代 大 第 二 个 方 释 中 ) ， 得 到 式 
(4 一 11)， 见 第 75 页 。 





4-3 几 种 简化 情况 


HRA RRHH (4-39) 或 (4-41) Æ “RRRA” SeMS AEH eA Rg FAA 
Si, "HE HRAT kE KIT É RREA EMG RA ERE YEA E 
TERG. EJEA ZE ESEA p EER ECE. a E k E f e o A 
此 ， 侣 理 的 情况 只 有 了 止 种 :如 果 考 虑 ARKA. MERE = 0 如 果 Bass0， 则 应 不 者 
E 4 站 的 影响 。 当 可， 在 以 上 前 担 下， 还 可 以 作 进 一 步 的 其 它 简 化 。 以 下 喜来 研究 这 些 简 化 
情况 。 

向 化 情况 工 基准 和 运动 为 水 平 站 线 飞 行 ， 即 gz =0。 于 是 方程 的 自然 形成 为 

(p=- vj)Advu- (x tgjdargAt- x’ Ake r Aia A Pe ` 
- y” du- [yt py Aa- -lpt yt Ah=y EEGEN 
-p3 dv- (utp uI dat (p ~ Hop) ET UE hni dd tusAP. 
vadu- vgdtt pAk=Ù0 

方程 的 矩阵 形式 为 式 〈【4-43) ， 见 第 76 页 。 

这 个 五 阶 方 程 组 可 以 用 过 研 党 飞行 器 的 水 王 骤 航 飞 行 的 畦 性 ， 畦 别 是 飞行 高 度 稳 定 化 问 
题 。 

简化 情况 ”如 果 在 挑动 运动 中 飞行 贾 度 的 变化 4 对 力 和 力 扰 的 影响 可 以 忽 酷 不 计 , W 

耦合 方程 组 的 自然 谋 式 为 式 〈【4-44)》 ， 见 第 77 页 。 
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(p- en Än (att gosta) Aat gesty dte xt dé, ef A, 
` Ié 
uy" da Dk Liipituë- inte än ru" (pt Cé sinda IA? 
SNE ry dP 
-al DEET CHE EECHELEN tus dP ; 
此 方程 组 的 华 阵 形式 为 式 (4-45) ， 见 第 73 页 。 
多 数 情况 下 还 可 以 忽 赂 Y AY” MA Meu ia, 小 必 方程 组 进一步 简化 成 为 





(4-44) 








-AV | 
a de | 
di Ao, | 
dð 2 
: xë x+ gosh y D 一 eos y Ar 
E SN 一 8 十 gsinf yi Ug 1 — gaint y Uy Au i 
tu HZ- Mt O- Sain Éasipel Ho Tis - pagsinds/va! Aw: 
` 0 0 1 0 "Ai 
` x? yP 1 
| a p | Aë 
-47 -y det 
+ A - å f ， | {4~46) 
B THR A An, 
、 0 0 
如 果 除 了 上 以 上 条 件 吕 外， 还 加 上 旋 淮 运动 为 水 平 飞 行 ， 即 #0 W fk M y e A 
f Av a" xtg 0 pai An 
d Au | =w e 1 0 | Au 
di Aos! Tiir 3-uy portas 0j do 
A8 1 0 0 1 D. Aù 
z db wf 
_ dz pË AÓ., 
y y b 
十 åz A p pP {4-47 ) 
Ha THY TI HaT HGY A Pe 
0 0 


在 研究 飞行 器 的 纵向 运动 特 福 时 ， 式 (4~46) 或 (4-47) 是 最 常用 的 方程 。 

简化 情况 下 ”在 研究 飞行 器 短 周 期 社 态 特性 时 ， 所 闭 虑 的 时 间 握 ， 速 庶 还 洛 有 明显 的 变 
化 ， 因 此 不 仅 Ah mh A. 4v 的 影响 也 可 以 忽略 。 此 和 外， 通常 还 假设 基 淮 运动 是 水 下 的 ， 即 
ba 一 0。 还 有 ， 在 短 周 期 摸 态 中 考虑 对 推力 4 已。 的 响应 是 汽 有 实际 意义 的 ， 故 设 4 PP 二 0。 
IERT: AIE 不 是 硬 合 变量 了 ， 所 昼 短 周期 运动 方程 写成 


a l- KÉ _ Öss 
Aa -— ZS Ae ` , H 
d 工 十 3 +y" ， lty" 
di = ; a |+ a, 146， 
go) (p A peep AD | o) Ia ER. 
` "tu ` Ë "TTT (im 


id-48) 


(Sp-r) 


一 ` -— e Wa 





a+] aft I 
SET 7 ` BE: 
gonep 一 D 











0 0 ` 
oSI ay ËA ey 
| ar D Bi "es bf Ze 
2d F afan Dean 
oF af +T DOEN 8 
T qf 7 zf 
aY ze 
0 0 
effi D SE 
(*a TEB 8 RA " DEM 
pf+I IT. 
HATER — phi afi 
# geog +A aA 
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DES IESSEN 





EE KE ENEE EE H 
rd An 1 Jo ur - 
d | | e ETEL Dis H II 下 > o DEM {1 19) 
Jo, Auf" "te dw: DCK 
向 化 情 识 下 在 飞行 名 的 长 周期 柜 态 由， 情 仰 迹 动 进行 得 很 缓 愤 ， 因 而 与 和 各 在 美的 
奶奶 第 可 以 忽略 m HL AA TE ci Agape. (Ga, E Jag. ELA 
"FE, IS IHN Adm das WA (1-39) 简化 成 为 
【六 一 | Gees a AD az AÂ, | PAte | 
47 Av- yt Ag- ip- gsind avyd =y’? Aå, tPA be 
=u} EENG da -ns do 1+ut AP. ] 
Jun "dé CDEN III 得 到 的 
Ar= (pl Avni Atu A Pe) HE 
EAMH: Jio A FA “FEESTE mv, 





, . HEEN 1 Hay” 
KE — KE 一 D = a - 
b ut Bi ur 
d o dan 
EIER Ha xX EEGEN B: Y 
È HS 8 DE 
De P.a 
nb He Papp HUO .5 
a? A uz H y u í] o0) 
T 4a R a a e Ge 
d ki p 一 Groty Jeppi ia USE 
= Hä Jl, 
dt AU uf Gelinfls/taJ- Af yt Bu AP, 


4-4 飞行 器 纵向 运动 的 特征 方程 和 传递 函数 
为 了 便于 过 阅 ， 作 达 里 把 纵向 运动 的 师 种 重要 的 常见 的 情况 下 的 灶 征 方程 和 传递 别 数 刻 


二 出 来 。 
首先 考 虚 由 方程 (1-46) 所 描述 的 情况 ， 即 不 若 帖 高度 变化 的 影响 而 号 pm 
HR. ATAR, SA TFE 
(4-52) 


Ye= goly Gao T Gini g Ua 


着 凡 把 rs: 简 生成 5， 把 上 标 @。 和 6: 简 写成 w 和 6。 
Art E 211 
ALG) -ast taast ass teste {4-53} 
(4-5379) 


Alsas taas? aast tas tag Û 


H 

为 特征 方程 式 。 

其 中 a, =l 
ga Emx? Hy HPT 
aa 5a [Lata go] (x ge Cy Oo ) 


BNTGE ET 


or=geH Hyp) (x HN EEN + (0 -9.0) "1 
FIs uT- eur? 
ou sel -YH is 0, Ann la, Tenue tie tu, Aur" 
以 升降 艇 偏 角 35, 为 输入 民 时 传递 画 数 为 





NS,ts) bass + host tps)h 
i Ie l _ a basit bis” Ch 3 -5 
D HO AG AC) 人 
共 中 KEEN 
h =| e” DTN a GARE -ixit ge Jy’ 
bi=[- (ye gi u+ gop lxt Egel nt) + 
+H x Hy rd 
hor vi x TI tt ~ EE u? 
Ng (s) tasi terasi tes He 
{ Io -一 z 一 - “3 20 J D. EE 
2) Ziel Ats) Ts) (4-55) 
其 中 ca 二 一 
C= äre Lian, x yt oi 
SEENEN EE pyt ga tx ut 
SEET KEENT 
NZ (sl) 二 十 ;十 
(3 Wg {(s)=—— u f ES Ea -5 
2) Ziel Als) Ala) i! dch 
其 中 BCEE 


el 一 [一 加 Re 十 Re)xe Mt xn y? ly- gso x pe 
eg SD YPU? (YT gsn xe tO- ixt geu +x uy? + 
+E lyt- Ysw) +t {x gy" lu? 
14) Eirsiëog 
Wg (5s)=W§ (s)-W3 (s) 
Watls)= sW (s) 


(1-57) 
W5 Va Tata Alf (s) + wacosd al (Cs 
A JPAERRAM FERAE b (LS, vi, aRt, ui, ei, 
RARR AE (4-49) 所 近似 描述 的 短 再 期 运动 的 情况 。 
frz A FAE 
Als}=aast toas ra = 人 0 (4-58) 
其 中 aa 一] 
ai 一 2 (He 十 AD) 
Qo= — HH" yA" 
以 升降 艇 偏 角 5, 作为 输入 时 ， 和 传递 丁 数 为 


rb c=- y? 
SEENEN 





KÉ ENZ — DEN + Ié 
" Pao iey EIS (4-60) 


Hp fiT RË- uyt 
foSy u? — nuch 
(3) 其 余 传递 两 数 


W3 = i W:s) | (1-61) 


W$ CY WE C) Wg (s) 


$ 5 模 侧 向 小 扰动 方程 组 
5-1 横 侧 向 小 扰动 方程 组 的 自然 形式 


DISK ah 《未 扰动 运动 ) 是 对 称 运 动 ， 其 中 所 和 的 横 侧 变量 中。、aw。*、 

Du, Pas Per un, Zopar yx 者 等 于 0 ， 因 此 有 有 
DEET TEE ET ET 
， 楼 侧 变量 的 偏 景 就 等 丁 谈 变 量 本 身 。 应 后 不 必 使 用 前 置 符号 4 。 

"enne seen, Ae (4-24), (4-25), (4-30) 是 彼 此 相关 的 ， 形 成 
烽 合 方程 级 。 把 副 寓 俑 转角 d. 和 方向 舵 翁 传 角 ó, 看 作 输 入 变量 , 放 在 力 程 的 省 边 ,把 侧 滑 角 
D, "iaa, (GER ov REA p 当 作 输出 变量 ， 放 什 方 程 的 左边 。 十 是 ， 模 
便 小 扰动 运动 的 看 合 方 程 组 是 

(nos -+ Qu ~ ZP) 8 ~ (mosas+T2Zoz)as- (nos+2Zeo)av- 


一 mm Gens Y = Zivë, 


-MEB HI Me) os (Ia A Mjo, = M6 +M, 


(4-62) 
f: d Oy d wy ER Aa 5 
一 MEB 一 ( sy J tM, jos 04 di -My Je, sM tM 全， 
一 xt lg, + EE 
共 余 的 方程 (4-29)、 (4-31), (4-33), (14-35) RIERS ts, SU mt, 
di — costy y ， 
= B- sinday 
P= P- osla 
(4-63) 
dzo = — Uncosd a H, 
costa. 





Y:=tgła P + cosb 


- ai 





CLN EROP G H EEA (4-62). RUINES TAA W AR (ò. =0. ò, =0) 
TEIRRA ETE. mH oSA DAT KIRAR EEN, o, iw, "AER A mO., y 
RA Zu or PPIE < 

经 性 化 的 方程 级 (1-62) DH if ER XL Sr Mm ing, 第 二 个 是 湛 转 力 
Sir, Stin im bp, 第 四 个 是 适 动 学 英 受 方程 。 

Tice, (IA mr ig, 

z= {ZÊ Qi} nmtis 


2 mv (FO Os 0,) 

















r AM Wi (1-81) 
Hen Halda 1 (r= Wrs in Öz ÅL) 
Hga sA, d : 
Ia, EE ER EI 
GEM \ {1-63} 
Popty i, Í: 
ARMERT ERMI ATO 模 侧 小 扰动 耦合 方程 组 {4-62) 可 以 写成 
" ierch $ Au, — (geosdy vap ze SO, , 
ELTERE ! (1-66) 


: ， o o — ð | Å 
e, Cesp tip od (ënne ein, TAO, 
— ws (ëss, t+ pr:0 


5-2 用 修正 导数 法 化 成 标准 形式 


由 于 惯性 积 1;, 的 存在 ， 横 侧 运动 方程 组 (14-62) 或 (4-66) 不 是 一 阶 线 性 微分 亡 程 组 
HIRERE: Amap A A RERE AT ERER E s, 的 作用 ， 可 以 采取 了 珂 个 办 
法 。 第 - :是 “修正 导数 法 ”; A E "PEER" o MEIRE- AI 

BEIA AE T FER EARI (4-24) 和 (4-25) 


do, dwy 











fa a de de Me 
E KI 
dmx dgy _ 
i grt fear TM 
出 此 解 出 
dos lt rts yl, o Af, ayid y) My 
dt o aI,- T tx -Tayif, 
> (1-68) 
doy TMy bf ey, L My T ai DM a | 
dt IER 8 TE ii, 
DEAM mu. EFR (4-13) 和 Ca , 并且 | 太 修正 的 动力 学 导数 
Mt OMY 5 
Hx Fa Tai, | 
Mit IOME i 11-69) 
Bg U Í, -jala 


(rp sto? GPE d. 


TARRE (14-66) wl DL 
(p z0) fP (oe Lazio, (1+, (god Un) BEE? 

SRSP Tip Hr Oa- Ng = 

- ná P 一 To, | CECR rä, tds, - {1-70} 

-s Tiga Oy t Dir 


vr Ay REAR e dEr A 





| Ui zf agta” li zf Or onfair D 0 D 
d o uË Die Diet 心 Or D D daf | 
155 ， Wen 了 了 
i x g ， - Üx d 
di Gy uf Hy Bd 0 Oy Jup” H” ò, 
Io L - ee 0 Iw .0 0 


企 正 导数 法 大 不 改 恋 运动 灾 遇 的 含 允 ， 而 且 计 算 也 不 太 复 杂 ， 央 此 是 相当 有 起 的 万 活 。 
ZEIL, Zo A Zr 可 以 忽 腾 不 计 ， 划 z?* 、z2 可 忽 腾 。 因 此 方程 (4-71) DS 


数 竺 阵 的 第 … 行 可 简化 为 
(2f Më 了 Cous fa Un) t 1-72) 


5-3 用 惯性 主轴 法 化 成 标准 形式 

"H rinRImik an Ah fr SS be, erzu Fo AC, WIER 
ESET 
BEHT RAA HF HER KAREO ARRE AnaL h E E ) 
或 测量 仪表 的 安装 基 淮 线 选 定 了 盯 机 仁 华 标 起 口 xjy;z1; 大 且 已 求 出 了 转动 刁 量 /i TT, 
Sigk SI... Ae, SEKR EI Le Ou, H FEMA” Jann e De 

ip a, 之 间 的 鲁 座 为 e。， 其 正 向 如 图 4-1 所 六 。 根 据 惯性 积 的 换算 公式 ( 见 后 向 的 


公 (4-81)} ) ， 有 


Jann 一 sinepcorepf(7y 一 LE Long ep — inep) u =0 


Ep ` -stn2es(7，- /Jeos2epf e "7 H 


Ar Ae 





图 4-1 惯性 主轴 的 确定 


由 此 可 以 确定 轴 x, 的 方向 
， 83 “。 





27 。， 
EE E see 
注意 这 里 和 角 ae, 是 有 正 负 意义 的 。 
人 然后， 按照 后 面 将 时 由 的 公式 (4-84) I EH ERRONEA opa Topi 按 公 式 (1- 851 
至 (4-90》 计 算出 对 主轴 的 气动 导数 对 ss。、 季 ?其 (M4).，(Ws*)s] 等 。 
最 后 ， 得 到 持 阵 形式 的 横 侧 小 扰动 方程 组 


CE 73) 





fej zË apy Ee l+ sin Zeus fl Us 8 >ò 
d Oep -lufo ez We 0 lo 
dt yp ， de? Ke: Ka D | Gab 
Yp "` ! 0 1 WP D | Yp 
A O dé 
noh Ha fo 
十 sp xp | {.1—7 1) 
nie uy Ia 
"ap "ah i SOn 
0 0 
按 此 方程 组 进行 动态 分 析 ， 得 出 的 是 万、@ pp、@yp、Y5 的 变化 规律 。 如 时 涉及 到 安装 
在 飞行 蜗 上 的 角 或 角速度 传感器 ， 则 必须 换算 到 安装 共 淮 线 所 对 应 的 @:、w;、Y。 


总 之 ， 膝 训 惯性 主 销 法 需要 进行 相当 复杂 的 换算 【〔 企 性 和 让、 气动 导数 和 运动 变 昌 ) ,不 
如 修正 时 数 法 那么 简便 。 
5-4 飞行 器 横 例 运动 的 特征 方程 和 传递 函数 


mist (71) 所 描 迹 的 横 侧 矢 动 方程 ， 在 简化 和 落 件 【4-72) (县 忽略 Z"* 和 
Zer) 下 列 出 其 畦 征 方程 和 传递 画 数 。 





畦 征 多 项 式 是 
中 (sy 一 Gd 十 Gas2 二 uiS 十 Go (4-75) 
D Atieizo,zitoszi koëikoustoazzü {4-75} 
HEIE 方程式 。 
其 中 ea 一 1 
ua -2f AE uy" 


ga SINT tuy H uR Hy" A get n 


人 一 
- afp tous ik NW vosthe (tp ns - uË) 
Go = -A costa [ighe unt — Sue") kann Ainsi 
D Fmptrsna, 作为 输入 时 ， 传递 丙 数 和 如下; 


Nes) besi tbs? + hystho 
Atd ACs) 





(1) Wi (Cs) 一 t-76) 
其 中 Pa = 


+ RA ew 


ED 

by (Enn pe) UA- staht i H+ Geo on) sr! 
ti -Haye ~ ( Gegen vs) ga Ju 8 

ha 5 {one nall- uiga- POTIN + 


2 (rosa Za lf? 十 DDT 








Nae*(s) east keaxt Heste 
H xis} = NA 1 " .j 一 了 
(2) aa) ALO At) { 7) 
-Erp Ca u?” 
c= pz t (Gëtt EL E 
o (uf wa) st EETA TETA 
Ca = {Yoo ia gäe AEN — fan) 
(Cs) eas? tes’ teste 
A 和 一 一 ER E (Est 本 
(3) 1 iz Cen . Are {1 78) 
Hpi Ce vu 
e, ier over Lef HO) dn 
pr= (Eust 一 Aen len ona bf — fuuss wie 
ey = (geste alt in RS) 
Ni (ei fast +fis+f 
tx 一 一 一 
{4) E (oi Als) Ats) {4 79) 
其 中 IEN 
EHNEN DATS + 
LI A — ëng Nët 
Ze (ungt — ARA At Af gba- Ho” u tgða)z + 
+ (HB sn? + Gett lëlle +t ab ow yd 十 
+ (ahut -aaaf talus ghe- oan het) rr 
(5) 其 它 传递 画 数 
根据 非 艳 合 方 程 -6 可 以 得 到 其 它 的 传递 画 数 如 下 : 
1 1 
Hären -Wr {g 
s, (5) costy A SA 
- t . 
It Lei: Ws (Cs) 一 SCH KE 一 inay (ei. 
(4-80) 


Wans) = — Oacus fl a 3 H Wéi 


- cost 
Hs(s)= 18d sl rei rg 


Bu s,s) 


当 以 副 届 仿 角 56, 作为 输入 上 时， 传递 画 数 同 上 MEE in, vin, UOR O, vis, 


= ND a 


So 惯性 和 气动 数据 的 模 算 

6-1 两 个 机 体 坐 标 系 

丰 实 际 虐 作 中 常常 有 必要 把 惯性 数据 和 气动 数据 昌 一 个 机 体 怀 标 承 民 息 到 另 一 个 机 体 鞋 
标 系 。 例 如 几 给 定 的 机 体 坐 标 对 到 正 嵌 性 轴 系 ， 由 络 定 的 机 体 轴 系 到 稳定 性 辆 系 【 于 和) 
HEKE FEADAR Ao 

h PALIER RR z Bema S KITAA RPD. gun", AEEA fi HLE 
Hp RI EE EE A y e io SE Ei E A B A ALe E R On Ai 
Ox, Ha ze fiz, z, Zë, xyi F Xy "0 HAEE Sr euk 4-2 所 示 。 

















Wa" Ping 
WS, HEA Oxiz, 到 Oxyz: 的 坐标 变换 矩阵 为 


cose sine D 
ul 一 Sine vose o] {4-8]) 
0 0 1 


HUHA (4-81) ， 可 以 直接 得 到 角 训 度 和 力矩 的 换算 公式 。 角 溃 度 


Oya Ort 十 15ine 


TEE ENEE Eh EE \ (1-82) 
ma TO] j 

HE 
Ai “MM, coset M isine 
Ai. = M, sine + WM, eose ! 14-83) 
M.A l 


6-2 HEERA 


A F Ca Ea ERR ag k aA (1-40), afate 
ME EEA RAS RH MEREK Ban RAE RN AL A 
(1-18) MER, BINENE | 可 以 按照 一 般 公式 《1-22) EAR. TE NB ZIEN H 
Rat, Tt, et, ERR 





Ts 一 了 ys 0 f I foiyl 0 
- f siya Í pe IW Le Ten Taya DW B: 
0 0 Í s 0 0 上 
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把 坐标 变换 公式 (4-81) RA RF H> 得 到 惯性 数据 的 换算 公式 
laa- costed zi + sini ef yl 一 2sinecosel yiyi 
了 ,一 rin2ef si 十 res2ef nr Zsinecosed x1y1 | 
EE 


Í saya 5 — sinecose(d yy — To) eoste ~ sin e) siy 


{4-84) 


6-3 气动 导数 的 换算 
dr DIE op e, OIEOëD DE AM, AM, AM ARTIR DRRR 121, 
(1-13) #1 (4-14) ] 
AM, Ais Mat 0 {x 
| AM, | An 0 | 
AM, o 0 Mee, 
此 关系 式 具 有 式 {1-18》 的 形式 及 相应 的 物理 意义 ， 所 以 可 以 按 式 (1-22) din 
(zen, (Mich, 0 CMa) (Me) H: 
Luc (Av, 0 上 aleco， OM'en, D ER 
0 0 MP), 0 0 MI)’ 
若 险 趟 《4-81) ， 展开 后 得 到 
CA Gei, (Me) oe (Mes MO) sineroe 十 (Marsin e 
(M2r), CM) cos se +MY” Mo) sinecoe — (Mgt Än sine 
(MEt) = (MZ eoste + (Mp Mizi,singeesg — (MEH) sin r EE 
(Minis MoI) co ee (Cp rM Zei sinerosg + CAM Zei, sinte 
(Meis ME a 


H ARE D E PER ERITA HH 

-AQ 0 0 -QO 

ae g e 
AZ: z9 ZW 0 8 


Een wy O: È Osn Oyn o, 的 省 写 。 同 样 地 ， 利 用 公式 (122)， 展 开 语 得 
zl 





(Q?) = (Qt ) eose — (UN VT sing 
(Y2), =Q? isine + (Y %7 )icose 
EST )icose + (Z5) sing 
(Zr), =- (219 sinet+ (Zr) 1cose 
通常 有 用 的 只 起 屿 后 归 个 式 子 。 
诗 干 由 和 角速度 引起 的 气流 坐标 对 中 的 气动 力 ， 虽然 也 能 号 成 
- 0 0 - QO 
RE 
Z Z'z 


ZY 0 EA 


{a—B86} 


JI O.V, Zio o, PAWA TTAR mA. MAREAAR (1-22), 
须 另外 扒 避 和 如下。 对 于 岗 个 视 体 华 标 条 (1 和 2 ) ， 可 以 所 出 


Q Wi OD - Q”: Ù y 
| ` ! | ` 8 r ` ! [e 
2 ZL"* z Ut; > 


Edel 0 
和 — d o - Q” Ox 
MV l - | Ü 0 re] dëi 
£ Lx ZY nn “1 De 


比较 以 上 上 师 碟 ， "äi. ërem a 


(Q7 =R), 

("ch K’), 
(Zant I) coe- {Zt ising 
(Z°) = `- (Zei, sing ! CZP} one 


(1-87) 


ZE, RUMMA ss, A. o, ERE IPLE ARER EEEN, mia Ime 
Su: MÊNE) 也 不 能 按 公 式 (1-22) 进行 换算 。 推 导 如 上， 
AN. M3 Mir 0 Ags a o 
fas] Ju Mr o ! 加 DN 
AAT 0 o je 
AM, Aix Mis 0 


sil ai] Ju Aix 0 | H 
o "AM." 0 -0 MÌ Ss, Ni D 
Estro, Gambia, Raa frets a e 


(M3); TATA TER CMO) sing | ` 
(MSI), = (MF) cose + (MY), sing 


(Mei, ~ (AM Zeie sins + (CMI) rose | (1-88) 

(Mr) = - (MŽ) sinet{ MSY) wose 

(Mira ME, | 
由 便 清 角 如 8 引起 的 气动 力促 也 具 丰 上 述 性 质 ， 因 此 有 

(HE), 一 Tag8)ieose 二 (sinoe 

(ME), — - (ME) sine + (CMA) cose | (1-89) 
最 后 ,由 å, SG dl amb, MAE 

(Zi) = (Z0) 

(Z6), = (Z9; | (1-90) 


总 再 亏 之 ， 气 动 导 数 的 焕 算 公式 是 各 不 相同 的 ， 必 须 分 仔细 地 性 虚 。 


88 ， 


§ 7 按 移 趾 性 坐标 系 建立 的 方程 式 的 特点 


在 美 亲 飞 行 力学 文献 由 最 流行 的 做 涉 是 采用 稳定 性 浴 标 系 { 套 看 文献 [1]、[2]、[31) 
Oxwywzw《 见 图 4~3) 。 宅 的 原点 口 仍 在 & 行 器 质心 ， 轴 xw 与 未 扰动 运动 的 速度 矢量 vs 在 飞 
行 器 对 称 平 面 投影 方向 一 狼 : 轴 ww 在 对 称 平面 内 ， 导 南 于 轴 xs， 指 向 上 方 ; 轴 zs 垂直 于 对 
称 在 面 ， 指 向 右 方 。 在 扰动 运动 中 它 基 与 飞行 吕 固 连 的 ， 所 以 臣 实 质 上 稳定 性 化 株 杀 吓 一 种 
特区 的 机 体 坐 标 杀 。 星 然 ， 系 定性 坐标 系 So Luise AN, 的 轴 xw 和 x 之 辣 《因而 办 
vu, E) 相差 一 个 按 栅 体 坐 标 条 衡量 的 未 扰动 迎 角 wa,x。 


Ne ka 





图 4-3 稳定 性 坐标 系 的 定义 


HS DEA, DEI, =0, 因 而 横 侧 小 扰动 方程 略 右 简 华 [大 看 廊 程 
(4-62) 和 (4-66) ]。 但 必须 把 惯性 矩 和 亩 性 积 以 及 各 个 气动 导数 由 原来 的 机 体 华 标 系 换 
算 到 稳定 性 坐标 系 ， 为 此 而 策 用 的 公式 是 式 〈4-84) ER (4-90)。 共 中 ， 如 果 “1” 代 天 原 
来 机 体 毕 标 系 (1!) ，“2” 代 表 稳 定性 坐标 系 《w)， 则 应 把 e 改换 成 一 2,%w。 而 且 这 样 得 到 的 
方程 组 也 在 是 标准 形式 。 如 果 要 化 成 标准 【矩阵 ) 形式 ， 还 必须 再 采用 修正 导数 法。 

最 后 指出 ， 无 诊 择 样 选取 和 山体 坐标 系 〈 即 普通 机 体系 或 怕 性 主轴 系 或 稳定 性 轴 肝 ) ， 对 
本 章 所 建立 的 锥 间 小 扰动 方程 组 都 这 无 影响 。 困 为 切 向 几 和 法 向 力 方 程 是 沿 航 还 坐标 系 的 ， 
与 机 体 坐 标 县 的 选择 无 关 ， 而 情人 匈 力 第 态 程 是 绕 机 体系 轴 — 的 ,但 在 不 同 的 机 体 坐 标 系 中 轴 
z 是 一 致 的 ， 所 以 纵向 小 扰动 方程 不 受 机 体 坐 标 和 对 的 选 基 的 影响 。 当 然 ， 如 果 纵 向 运动 的 力 
上 程 是 挡 机 体 坐 标 柔 建立 的 【做 看 文献 [3]、[23]) ， 则 机 体 坐 标 柔 的 选择 就 会 对 纵向 小 处 动 
方程 的 形式 存 所 影响 了 。 


~ K9 >» 





第 五 章 ”在 非 平静 大 气 中 飞行 器 的 运动 方程 


SI 关于 风 对 飞行 器 运动 的 影响 


在 第 二 进 81 中 己 指出 ， 习 行 器 的 对 地 速度 《飞行 器 质心 相对 于 大 地 的 速度， 简称 地 速 ) 
va HIRE (KRR DIFESA RE) u -BAR wAn FRA EE 
二 十 次 {5- 1) 

BEEP ET EN t RIESA RE 
ZE TH GETA ENZ. ， 介 在 建立 运动 方程 时 假 
定 了 大 人 乞 是 稻 止 的 ， 亦 即 风速 W= 0。 

详 际 上 大气 往往 是 不 平静 的 ， 即 有 网。 风 对 飞 
行 圳 的 运动 起 作用 ， 首 因而 产生 多 种 多 样 的 影响 。 
例如 ， 克 影响 飞行 区 的 续航 性 能 、 起 飞 着 陆 性 能 、 
飞行 器 的 平衡 和 操纵 ; 风 还 影响 飞行 器 受到 的 裁 地 站?、 罕 省 4 与 内地 ww 的 关系 
荷 ， 引 起 飞行 中 的 题 壬 ， 影 响 乘 洛 的 名 适 性 和 结构 的 首 劳 寿命 ， 风 的 作用 影响 KIL3e ti DEN 
行 跟 逊 、 射 击 等 使 命 ; 至 于 风 影 响 导弹 的 命中 准确 度 更 是 明显 的 。 所 以 对 二 飞行 器 的 设计 和 
使 用 来 谱 ， 风 的 影响 是 十 分 重 旨 的、 不 可 名 视 的 问题 。 本 章 的 日 的 就 是 建立 存 非 平静 大 和 气 中 
飞行 名 的 运动 方程 ， 利 用 这 些 方 程 站 有 可 能 研究 风 对 飞行 的 各 方面 的 影响 。 这 里 仪 服 于 平面 
大 地 的 情况 ， 当 必要 时 读者 不 难 将 这 些 方 程 推广 到 球形 大 地 的 情况 。 

一 个 很 重要 的 概念 是 运动 的 惯性 力 决定 于 地 球 Y ， 而 作用 在 飞行 器 上 的 染 气 动力 及 和 
发 动机 挫 力 P { 火 稍 发 动机 除外 )》 则 央 定 于 实 速 u 。 用 直子 表示 ， SR, 

dY 

Tan 

E HLE Kiik, Mach 数 、 Reynolds MARED A Pe Cys Cza Mas Mys MAUJ uo EBE 
速 & 决定 的 。 

再 回 题 一 下 ， 在 第 一 章 定义 的 坐标 系 中 ， 地 面 坐 标 系 .$s 和 机 体 坐 标 系 S ,是 与 风 的 存在 
EKA RRS ,是 以 地 速 v 为 基础 的 ;而 气流 坐标 系 $。 和 中 机 体 坐 标 系 S$ alku 
为 基础 的 。 在 从 个 角形 中 ， 鼎 航 角 革 、 米 人 委 角 站 和 湾 转 角 ?与 速度 无 关 ; pa at fym 
航 迹 倾角 4 决定 于 地 速 ，; 迎 角 “和 侧 滑 角 1 决 定 于 空 速 上 ， See, 仪 在 无 风 时 寺 有 
EL Rh EEEE S, 与 气流 坐标 系 S$, 之 间 的 肉 一 角度 ( 见 图 2-5) 。 

有 有 : 风 对 飞行 的 影响 方面 ， 有 两 类 不 同性 质 的 问题 。 第 一 类 问题 是 有 规律 的 (虽然 是 想 定 
的 )》 风 的 作用 下 飞行 器 的 运动 。 这 属于 确定 性 的 问题 ， 第 二 类 问题 是 连续 天 气 案 流 对 飞行 器 
运动 的 影响 。 这 属于 随机 性 的 问题 ， 即 系统 对 随 宙 输 式 的 响应 问题 。 

另外 ， 在 运动 方程 的 形式 太 面 。 有 两 种 不 同形 式 的 运动 方程 。 第 一 种 是 全 RA frisch 
方 企 。 沁 对 反 速 的 大 小 漫 有 和 任何 限制 。 泡 鹿 适 用 于 确定 性 的 风 ， 世 适用 于 随机 进程 的 克 。 方 

。 DD e 


的 





zc RireitPtooh mec 


程 中 可 以 不 仅 反 上 映 风 束 本 身 , 而 昌 反 册 风 速 棋 度 。 这 种 非 线 恬 的 点 程 只 能 用 数值 广 法 解 出 ;第 ， 
二 种 大 人 坊 量 形式 《或 小 扰动 形式 ) 的 方程 。 它 急 以 无 风 时 的 对 称 定常 直线 飞行 作为 基 淮 状 
态 : 而 把 出 丈 引 起 的 附加 力 和 力 捉 作为 外 加 下 扰 硕 。 当 然 这 只 适用 二 风速 不 大 的 情况 。 这 种 
上 各 对 于 用 频谱 法 解决 飞行 器 对 大 气 率 流 的 响 店 问题 特别 适宜 ， 虽 然 它 同时 也 适用 于 研究 飞 
APARTIR ETER, (nët ei) KIIRE o 

E 82 JET A a AM ETAK ETARE ai A Ri 0 Rue E EA i EA 
Du. $404 T EAR EIJERRA OR KI E A i SURE -Ahipu BITE 
Wäzd KETARA A AIE r EET ng 2g. 


$2 在 非 平静 大 气 中 全 量 形 式 的 运动 方程 

2-1 概 述 

作 关 于 风 的 影响 的 第 El (REI) ho ARE w ah A 
IC 0 E 90, JAER iA PNE wa, wys waa o EHRE 
Ip, AURRA, ARAR ie RIE EBER. ER REBT A RNE 
HAIR IREA R” ERE ERG RREK. wE EER M pE EDMEE 
情况 。 在 对 飞行 试验 结果 进行 复 楼 计算 时 ， 风 的 数据 应 以 试验 中 实际 测 其 到 的 人 气 策 件 为 俯 
Hi 

AE EIEE THAR IEA KEAS A ATE ERER EEH. AR 
在 第 二 章 那 样 ， 有 儿 个 不 同 的 方案 。 下 面 分 别 讲述 。 

在 所 有 方案 中 ， 第 二 剖 的 转动 动力 学 方程 (2-41) 和 转动 运动 学 方程 (2-16) 仍然 保持 不 
Zo DU. 仅 列 出 质心 称 动 的 动力 学 万 程 、 质 心 运动 学 方程 以 及 有 美的 英 系 式 ，。 而 且 为 了 
节 阁 篇 幅 ， 仅 写 出 年 阵 形式 。 在 实际 应 用时 ， 读 者 自己 可 以 把 这 些 方 伍 展 于 。 





2-2 TT 体系 的 运动 方程 
有 风 时 ， 机 体 坐 标 邓 中 的 质心 动力 学 方程 (2-31》 和 质心 运动 学 方程 (2-44) 仍然 可 
用 。 特 动 动力 学 方程 (2-{1》 和 和 转动 运动 学 方程 (2-16) 也 仍然 可 用 。 
根据 速 庶 关 系 式 (5-1) AAA u 在 机 体 坐 标 和 中 的 分 由 
Hyg Uye U y 
fee Jofo] Bafe] (5-2) 
Hzr Use Wed 


Zu, a, D [EER (2-55) 的 推导 ] 





H 
tga = -= ` 
St 
I 
sin Hl = Sé (5-3) 
1 
- | 
u= duž, ni, +u, / 


总 之 ， 束 个 体系 由 方程 〈2-31)、 (2-44), (2-41), (2-46), (5-2), (5-3) 给 成 。 


* OI 


可 多， 与 无 册 时 的 运动 方程 相 比 ， 仅 有 少量 的 改动 。 在 有 风情 况 让 委 甩 这 个 方程 体系 足 
最 为 向 便 的 。 所 以 它 是 利得 推 着 的 。 





2-5 “DT 体系 的 运动 方程 


行 风 时 为 了 把 气动 力 由 $s 宰 换 到 $s， 不 能 利用 式 (2-13) 给 出 的 短 阵 BY 
5 ,传递 . 所 以 地 而 毕 标 系 人 Ss 中 的 质心 动力 学 方程 应 写成 


dv, a'di Propp -Q 0 
deu, dien? J+ dd Y II d (5-4) 
du aldit 站 


0 


， 了 必须 通过 


质心 运动 学 方程 为 


dxp 

dt "ri | 

Zi: SE (5-5) 
dën j 

di Wad / 


转动 方程 (2-41) 和 (2-46) 仍 保 留 。 


HUI FAREED E 
Hys Uxa T Wyd 
Lal {5-6) 
H. Drd Wyd 
然后 利用 式 6-39 HBe, f; u, 
罗 个 方程 体系 由 式 (5-4), (5-5), (2-41), (2-48), (5-6), (5-3) 组 成 。 这 个 体 
Ait TT 体系 复杂 - 些 ， 但 其 些 情 况 下 还 是 可 以 和 采用 的 。 
2-4 HT 体系 的 运动 方程 


HAR EEEX, o (2-10) ARDT BiM 不 再 存在 ， 所 以 质心 动力 学 
方程 应 写成 好 下 形式 ; 


roidt os 人 -Q 0 
"| od idi Le il Punos |+ d Y j} d ns] {05-7) 
-veost dy, dt 0 pa 0 

Je (2-43)、 (2-41), (2-46) 仍然 保留 。 但 儿 和 何 英 系 式 〈248) E (2-19) 有 有 成 


Dei 


机 体 坐 标 季 中 的 速度 分 量 按 如 下 公式 计算 ， 


He U Wad wi 
[e L le 0 J- j) 5-8) 
Ha: 0 Waa 


D 92 D 


然后 接 式 《5-3) He, fou, 
整个 方程 体系 由 式 (5-7)、 (248), (2-41), (2-46), (5-8), (5-3) 弓 成 。 
Jiii (6-7) 的 展开 形式 极为 复杂 。 所 以 这 个 方程 体系 起 最 复杂 的 ， 建 议 不 要 合川 它 。 
Zb 有 风 时 飞行 器 的 对 称 运 动 方程 
上面 列 出 有 风 时 改行 并 在 钳 垂 平面 内 对 称 运 地 的 方程 组 ， 从 这 里 穷 易 了 解 各 个 霄 种 体系 
的 实质 和 特点 。 而 导 这 也 是 实践 中 常用 的 。 
HEA, DAES E ER 5-2 中 表 小 。 





IC 


l. TTE 














m( SE 一 asuy) = Peosps ` Qeose +Y sing - mysing 
du, - 
wl - 十 ode) = Pasing + Qaing +F ceog a Gens 
dor, Ka 
fe dr TAM.: 
xp . 
di "= Ved E D cos- uv, :sin 
dyp ， 
SEO S=Vyd = Un H vr cost (5-9) 
dë 
di "D: 





Wes = Ur ne giny 
uyr SUyr + wdsing — wyarosy ! 
— 2 pO 
u= W už Tuj, 
Hai 


lag = 
> Hy: 


us y uÉ nf, 


Iyd 
` H 
: ing -二 -一 一 一 一 


2. DT 体系 
m ,i = Paeh- pp) -— Qro- a) -Ysin -or) 


m oge ` Paint Pa) Jsi- g) tY eoat- a)- mg 





DEE 
H 
Jaen 
, D, d 
dr | 
dyo _ 
dE "ed 
D (5-10) 
i Siha : 
Hued T ui, d hud ~ (Uya 一 ue Ji 
Hye = — {Uxd — td laind (Uyd 一 i ya Jos 


u= N uł, Lui, 








v= y vžątvłą 


Uyd 
tg -- 
Uyd 


3. DT 体系 





m se = Pealt- U +p) — UHeostiä Hoi: | 
| 
— Land Hoi mosing d 
8 
mot, = Panl- Lei - siui? t-a) + | 
+ Yoos{(yd 8-ol- mgit | 
do, ` i 
i, di =5M, j 


fe 


den 
SB Geen fl 
di Ud vreo 


dup 


-- gy Dad using (5-11) 
dë ` 
di 二 fda 
Her 一 (st Wrq) eosd Lag ~ wya)sing 
uyi = ~ {Ura T Wasin + {va ~ wa deont 





u= W už, +uż, 


Hyt 
Hee 


tgi = 一 


§ 3 有 风 切 变 时 的 运动 方程 
3-1 概 述 


在 地 面 附近 ， 由 于 空气 粘性 的 作用 ， 凤 回 往 往 是 不 均 与 的 ， 形 成 如 图 5 3 所 村 的 忆 剖 
Wio RPR ROAI E. REER, MOREX KITAN E AiE a KITA EAR MA 
2. kA E w E E A g 





h 
w: wih) (5-12) 
只 要 给 定 了 上 上述 医 数 关系 ， 这 就 对 应 二 在 已 Sab E 
行 的 情况 ， 完 休 直 以 条 用 本 章 8 中 建立 的 有 风 有 时 的 运动 方程 W CA) 
组 ， 计算 和 分 析 飞 行 器 的 运动 。 
企 理 论 分 析 中 应 用 较 广 的 怒 一 种 方法 是 考虑 风 切 变 的 影 
响 建立 飞行 器 的 小 扰动 线 化 方程 。 
图 3-3 MUE 
3-2 考虑 风 切 变 时 的 纵向 运动 方程 
我 们 仅 研 究 风 切 变 对 纵向 运动 的 影 啊 。 
飞行 的 地 速 ¥ , Cu 与 风速 w 的 关系 仍然 是 式 (5-1)。. 
LURRA 
KEE) (5-13) 


Eh CEAR U 对 水 平面 的 倾斜 角 ( 见 图 54) 。 
PA (5-1) 投影 到 气流 坐标 条 中 ， 有 


Dee Ha Wya u penal u 
[oe jefea DER 0 Ll wen 
Vog KÉ LD 0 D 


应 由 第 一 章 中 质心 运动 的 一 般 方 程 (1-36) . Sum REIRES RHEE R. E 
.的 角速度 分 最 为 


es Ob e 





TEAM 


d u + wos u DÉI ui mens, 
lan J- 0 | | — wsinĝ u ] 
0 da 0 


注意 其 中 于 标 ss 表 示人 针对 称 给 阵 。 上 式 展开 后 得 到 











du + dw g : 
f di 人 ` i Toast o | Prl- Q- mysingdy 
| 
D _ dw Auf, PELE = | Psinfa t gp) tF - mgeosh u | 
f di " di ! 
l 0 | 0 
显然 ， 共 中 第 三 个 睫 程 吕 以 抛弃 。 
由 二 
dw dom dh 
di dh ut 
而 在 风速 为 水 淮 的 情况 下 ， 有 
dh 
di KE? sin u 
HHEN 
do 
7 dh 


表 修 在 给 定 高 诬 上 水 平 风 速 的 铅 垂 梯度 。 寺 是 上 妹 方 程 可 以 化 成 


v HD e 


Wa ~ 
~ CO ` 
rv E 
Za A 
4 yu sinte coss + 9 siau) Foetz tel- Q | 
Rat R (3-14) 
E - quint Os + geosd a) z Peinflat g +F | 
ke ` 
a J, ŠPE M,- Per 
d o | 
di "= , (5-15) 
| 
fi) 1 
ci H sin, , 
a=ġ- d. 


3-3 所 切 变 下 的 线 化 方程 
接 第 四 章 中 采用 芍 一 般 上 方法 ， 把 方程 (5-14) 和 (5-15) 线性 化 ， 共 结果 是 


dA f 1 - ] 

Sr = t IT Drot Gas pp- Q"]- 一 z DEET u+ 
| 
| 

一 L- Pasini äg Fpp- Q"]Aa-— huyo 2 suy + eos fug) Afa i- | 
| 
toi [Přeoslaatga)-Q']Ak- A QAS | 
1 m COSL Eag {p d m z | 
dAl: _f_1 Fai D - 1 ain Ä i 
di TA mig LP"sinf dyt pp) EFSI + va Dain Gun) Aur 
+- l [Paola rp) reda +( 2 einfa t 9 sta PERS 
Mia yos Ga z)! 4 uy P e% : uy eP; 
1 , ， 1 ， (5-16) 
+ mug E seint Ga HEI +Y Dart Aĉ, i 
I 
dAdo: 1 H u LI E LI a - | 
kd E= po (MP ech Anc MAr + y: Mida t i 
1 o 1 , 1 
+y M Aost y (Mi- Pep dht) -MI dós 

d 

At = Aw, 

dh WM sain ra A u | sal Kë 

da= A8: A0, 


显然 ， 方 程 中 未 考虑 发 动机 的 油门 控制 。 
引用 第 四 总 式 《4-38) 所 定义 的 级 向 动力 学 导数 符号 zx" ，x"，…，43= ， 但 其 中 的 za 现 
在 应 改 戌 ux。 十 是 可 以 把 方程 (5-16) 写成 莘 明 形式 ， 然 后 再 化 成 状态 方程 式 的 奈 淮 形式 


e DI e 


AC 





EE 


si 


Sa DAN 


0 
名 
(rage t "ënn In 一 
- 
Ka 
ža ggo — Ër Out H 
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Zn ggwn i =n 15008 一 

















la Až a3 
8 ， D , 
eF | 
j zf- | 
— X 
\ J 
9 ën prostr 一 Xa gurs 
[ 0 0 
+ 
Ru DM (ozmse ->* gms ik -人 一季 Logan Hafi sH — an 
Eai uis f Ts Wé Io (* TS 一 f) - 
1 E guts H afi = offe 
0 Ka? Aa KA gep 一 KI nr penei : — ach 
b ypo 
i i 
i 
| zep 
iAP 
i 
el P 
I 
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这 个 方程 组 有 如 下 上 点 值得 注意 : 

第 一 ， 与 第 四 章 的 方程 组 丰 比 ， 除 了 状态 变 划 中 的 Ao ordre Ae 以 处， 系统 
BHR Eraus zap hum, drop (MEAE) 时 ， 两 组 方程 就 重合 了 。 

第 一 ， 这 里 所 有 导数 的 岱 交 及 开 达 式 均 与 通常 情况 的 相同 。 蛙 别 EF, oi. MRR 
五 大 气 密 度 梯 度 p* 和 吾 速 梯度 u 有 关 ， 帮 不 涉及 风速 樟 度 dzv7d8。 当 庆生 可 以 忽略 时 ， 

变量 上 jj#9 吕 以 去 掉 。 
=o AIH Eye ENOTE D. ERR-I RS-A mR O EW: 
BH HEEE ek KIT, BU, 0s R A A KITAR AE ETET n o 


§ 4 IPRA A H fE ERINE h iE 
4-1 E 


-发 于 RTRW Stater 是 一 个 十 分 复杂 的 问题 。 它 包括 ; KAEI ARARE 
LD Soe, CC Chic Kim rr: RÉG Ah 
SEMET, DEE EN KITARA Y E ESE A E EN i Si 
E GA Ap me FEAA, kim LENA AA. Bic LS 
ITA RDE ah A E Ra 

ERIE KITARR AATRE M, A TE bR, Bu ir H e Caus 00 Ay 
(HRA ERA HHD o FRH; aam, MEERE. a, F KITEA HH. ef LAA 
ARIER pE FA” o MEES AARAA m EI, ATEA E. i 
最 后 这 个 假设 是 必 要 的 。 

ATHA KITAA TAWE BARA KHAA AR EARS A Beeb 
RJAR. HLERA WEES ER mwa AAEE ARANA ENE. 

现在 ， 仍然 用 代表 飞行 器 的 地 速 ， 用 w REAR TAERE RAEE. t u akM 
CREED” Rizo um Liebe, MAMEET xan 
Du, WERO” mM LRA t RAR., E ER ir aE ab RO xyz, 
UE i E E e A en AEF EE 

Wys Wys W5 f lxs Y» z, 1) c5-18) 


HAAPE SERR ER E KITA Wat, RHR EE EA R EALER 
IE E PA EBR fio 

用 于 飞行 器 在 y 方向 【 即 大 体 上 错 垂 方向 ) IRTEERARA RRN hRS., 所 
以 可 以 不 考虑 雍 访 速度 在 4 Ami Ee, B 

D Wys We= fxs zs 1) (5-19) 

J EALA a DHE EE A LS, 

HASERREA TERE, MARAT ETARA EAD GA 
Zi, mP RARS. wDI a eE EE a e KI ootte. dr 
表示 成 





= DO e 


ER x SC 














W(x Tt) = w0, 0,7) + <0, Ü, 1)x 二 一 

ayla z [J= wd 0,2) + aer EE äng | (5-20) 
， | 

lt D (0, WÉI tx T Zone = 


为 了 以 后 节 写 方 恒 起 兄 ， 条 用 如 下 的 简化 符号 : 
We =W 0,0, F); wy =w, 0,0,7); w :=w:(0,0,1) ` 
dru, Ow 





LI ee 一 ag 《0， 0, t), Wer T gd D, DN 
Ow Jw 
wys a (0,0,4), Wya =; {0,0,1) | (5-21) 
a (0,0,), wa = D mn ` 
x b 


4-2 素 济 速度 的 气动 等 价 作用 


大 气 寄 流 速度 的 作用 实际 上 就 是 改变 了 飞行 器 上 的 相对 气 判 ， 从 而 改变 了 飞行 器 所 受到 
的 气动 力 ， 耐 这 种 附加 的 气动 力 可 以 通过 飞行 器 运动 状态 量 的 等 价 改变 来 表示 。 人 分析 如 下 
( 见 的 5-5) 。 
x AARRE ww 的 作用 就 气动 效果 来 说 相当 于 减 小 前 进 速度 ， 这 个 效果 周 vy 吕 示 
pg 一 一 as (5-22) 
这 里 w. ige -210 中 规定 的 含意 。 
# 方向 的 这 流速 度 w, 的 效果 是 增加 迎 角 


Cg = yl Us {5-23} 
DG IGE E w. EH EAS H 
a= np, Ug (5-24) 
FERREA We e, ETRA REA A E F AR UC E 
man Elie IXE ~ Uy (5-25) 
HJE gw,y7az 相当 于 下 的 滨 转 角速度 
rg TIW y O5 Wy (5-26) 
梯 座 Bw ,+9x H T E UL fE BE 
Opg PDW, OASE W; (3-87) 
HOJE Aen, /ä: AF fa E f RI 
Wyp T T Oli DS -Wre (5-28) 


以 上 两 个 偏 航 角 建 度 的 气动 效果 厦 不 相同 ， 所 以 用 下 标 1 和 2 相 区 别 。 对 于 无 后 壤 机 村 
的 飞行 路 来 说 ， 梯 度 3uw :/3x SE RRn ahh, MEE äu, /9z SEL SE 
气动 力 : 对 于 后 控 机 杭 来 说 ， 两 种 梯度 都 会 影响 机 村 的 气动 力 。 因 此 ， 汪 两 种 等 价 角 速度 


- 100 e 


Wyr 


Mis 


DC 








dw afi 


图 5-5 大 气 素 廊 回 论 的 气动 效果 
up ERE SS NU zb ah, Mixt, Mixt, Mir, Mpt? 
HE 3w./3x TELE A EE EC CEET E OS 
梯度 Aen, /3z 的 效果 也 可 以 忽略 不 计 。 





最 后 ， 还 有 
å 二- REF — di me /vn) — ER 
7 gi E EIST 
BU 
KËSSEN {5-29} 


D 101 a 





4-5 纵向 运动 方程 


到 第 四 庶 中 已 经 线 化 了 的 级 向 方程 {4-21) 、 (4-22) 和 【4-26) ， 厦 把 和 案 洲 速度 的 影 
Wo, os e, Ġo BIRIA (5-22), (5-23), (5-25), (5-29) SEH kR. mie HAE 
KA iip iga i REN Eo 

方程 (4-21) 修改 成 


Pasint en pi) des [QCA -t dQ dh 





+ Gel Ae, r) ENEE mgosba AË (5-30) 


方程 (4-22) 成 为 


ZA 
cai =[ Pe Au rw A PARI APolsinCos i pap) t Parol ty ts! AJoa"t 


Loi än 


ENEE KEE rd Aat POL Fél titw): 


YOC Ae, Wyr) Vie BE mysinfy AH {3-31) 
方程 【4-26) WD 


cl Aa, 


d di 


= dv- w) t+ Mi Abr Maf Jaat se) 
| Met Ata dw) + Me Ao, ~ we)t MI A’ 
-[#"{ Av- w) +P AAt AP. Je: (5-32) 
gem (ls, AE o E KIPR: 和 是 地 速 的 倾角 : e aie E LEA 
ZAHER. RIAIR RRE. Aea 是 它 的 迎 角 坪 最: w/vs: 是 MEW 
FEJERE R, 所 以 Acata = dc 就 是 气动 迎 角 增 最 。 其 它 的 最 (dr 局 sh 
CA 一 WAX HGRA. 
现在 ， 除 了 利用 趟 《4-38) 所 定义 的 动 hRS, xt, e DL, 再 补充 定义 两 个 动力 
导数 


xt- im e, 
x Uz, (5-33) 
YSS Y] mug 


ru, B (5-30). (5-31) 、 (5-32) ， 再 扣 上 另外 二 个 方 保 ， 可 以 写成 如 下 简明 形 
式 : 


a 12 e 


os 





di 
3 LI - 
— el A D. Aerm, +x? 7 
dA? `, d un at po H ja oyi d di WÉI N 
de SY dvo duat ginda O ot l Gp T dos 
, w 
— ui jA vis Ié. (uf An, trau, + yi + 
. a 
(Y-Y Eya 
ddw, v a a e A H , 《5-34) 
人 
wW . 
一 [PP 
Ua 
d Aè 
Au Fdo 
d Ab , 、 
dt ` sinf y Ant urneosd AË 
PAEA dAd 
= dw: -R 
di di 


WEEER EE 3 Lt HAAR nur: MAR (5734) lg E Au TF H) 
HIPER: 


An 
dAn 
di" Je. 
! Ad 
『 x” Xe 十 grosh y 0 一 eost y ~ 
| — y" — Mi Dain fl yl Uy ] ~ Geif sis | Jes | 
| oy ne (ygsingwdus) utrou negsingwi vn |i dos. 
0 0 1 0 St 
wf x? ` 
I 
SU y” TERE 
t D a P oof | 十 
uiz — Hy z HaT dk?) | A Po 
0 0 o’ 
一 x” 全 和 0 ` 
1 0 Wx 
HN "Hai = 
+ . o . wy (5-35) 
一 
Wys 
0 0 0 4 
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VW TE 


d 
RW x= Å, x, +B u, +409.: Lä 36) 
共 中 vz An daa Aw: ArI REKA 
u= Aë. AP] EHRE 
g Siwe wyp Wyrd EURE 


下 标 d RRA o 
有 了 了 方程 《5-35) ， 就 能 够 计算 飞行 前 对 大 气 案 流 的 纵 商 响应 特性 。 


4-4 横 侧 运动 方程 

取 第 四 音 中 线性 化 的 米 侧 小 摔 动 方程 {4-23)、 (4-24) 、 (4-25); Eh A E Zp EIE 
的 效应 Ba. xa Giaëix Oygg 的 公式 (5-24), (5-26) D (5-27). (5-28), 就 得 到 在 
ZS Sum ak, 

方程 《4-23) 成 为 


d AO 
di 





M Ug 一 mie Dës A on, T MU dl y 


C- Qat ZA) (ABa EE 
tZ Aag t (Z wss (Zaw t 


十 之 2 A4. mm enen Ay (3-37) 
方程 (4-24) 成 为 
d'Ae, d Aw, 
deer iwer 





=M${ 4 Ba Mit do, km, At Mi Aos + 
Da 


+ Minitom,, Mini, + Me As Mis Aë, (5-38) 
方程 《4-25) 成 为 
d Au Y _ Í d Jo, 
" di We di 


=Ms( 4 Bao Ea) Moe Ae, tr wy d+ My dod 
+MG w, My w +My AÖ +M? A, (5-39) 


除了 利用 式 〈4-64》 和 (4-69) BELEA MA DF RA BEA k 
MEU eyii M 














DE 3. ZETI, 
KE Me +t (Ts ys) My | 
fa 一 l, -IZ d, i | 
. {3-40} 
GET E . "" 
WT 
aya M+ /Te ) Me . 
Py 8 JI. cf Zait, ` 


» lij e 


e" p a Dm, PERA (537), (5-38), 
AERAN MR T AER e A E A R R EE 











(5-39) OG 


"bk, gd Jm 上 过 


Pa Cof oe 1 Drtaflar Us Pua 
d DO, HÊ in no 0 ox 
一 7 一 ， = ,本 
aE "oe, (nf wi Hp” 0 | wy | 
wu Jo 1 -tety IW Au. 
0 ze -f jug 0 0 D d 
E A å, — Bia Dy Gaul Gu? a 
+ “x tx Ge Ja lx fx fx ， b : (5-41) 
Dei Dei Lë 一 Ra Hy” Hoi Asil we 
0 D. 0 0 0 ü Waa / 
成 者 当成 
d — A B | C 5-42) 
di LÉI KÉ ell, "7 P LA 
其 中 vz Cla Ox o, Y) 状态 矢量 
nt, All EHRE 
Js= W, Wyz Warr Wya)! FAKE 
Egi F. o ARRE w, HR ppmizpnu"v, MPRRM L, 
Min, RR” SE (5-41》 的 右边 最 后 一 项 就 简化 为 
『 — 2f jug 0 ü 
I 
` wii air wëll . 
p: R a js (5-43) 
一 Hy Vy Hy” Hy” 
Wax 
0 0 D, 


4-5 运动 方程 的 运用 


RREA HEE (5-35) 和 (5-41) 可 以 用 来 计算 和 分 析 飞行 器 对 任何 性 质 的 宾 风 {包括 
随机 突 风 和 离散 突 风 ) HEM (AERAR bT KIRE) ， 但 主要 还 足 计算 和 分 析 飞 
行将 对 随机 突 风 〈 即 连续 大 气 案 流 ) 的 响应 特性 。 其 主要 步 县 如 下 : 

第 一 ， 按 照 运 动 方程 ， 挫 导出 状态 和 失 量 对 疾风 和 儿 级 的 传递 画 数 短 阵 : 

第 二 ， 根 据 设 计 规 范 ， 崩 定案 流 的 尺度 和 蝇 讼 ， 半 且 按 选 定 的 大 气 案 流 模型 ( 例 如 Yon 
Karman 模型 或 Dryden 模型 ) ， 俏 定 大 气 案 流 的 各 分 量 《 亦 是 方程 (5-35) 和 (5-11) RE 
风 人 欠 量 的 分 量 》 的 频谱 密度 和 交叉 闫 谱 密 庶 ， 亦 即 突 风 的 频谱 密度 矩阵 ; 

第 :， 接 巾 线 性 系统 的 随机 过 程 理论 方法 ， 由 以 .上 两 个 录 阵 ， 计 算 则 响应 时 《 杰 即 状态 
RO 的 频谱 密谋 征 阵 。 进 而 计算 响应 量 的 各 种 鬼 率 特性 。 

K 上 这 个 问题 的 详细 理 认 和 方法 ， 读 者 可 查阅 文献 [1]，[6] 等 。 
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第 六 章 ”飞行 器 运动 的 计算 和 分 析 方 法 


Si d 言 


以 上 备 章 中 已 经 建立 了 各 种 类 型 飞行 器 在 各 种 飞行 笨 件 下 的 运动 方程 。 所 是 为 广 社 出 对 
痢 行 强 这 计 和 使 用 有 意 文 的 结论 ， 还 必须 设法 求解 这 些 方 程 。 术 章 的 目的 就 是 讲 进 《行将 适 
动 方程 的 求解 方法 。 

运动 方程 的 性 质 不 辣 ， 则 求解 的 方法 也 不 相同 。 我 们 将 考虑 以 下 二 种 情 闹 : 

第 一 种 情况 是 线性 微分 方程 弓 。 主 要 是 常 柔 数 线性 微分 方程 组 ， 第 四 童 中 建 人 内 K 4 
葛 小 沈 动 线 化 片 程 组 就 是 杠 丁 这 种 情况 。 这 样 的 运动 方程 可 以 用 解析 法 求解 。 它 和 包括: RH 
飞行 器 运动 的 特征 杨 和 特征 矢量 ， 在 知 始 扰 动 后 飞行 器 的 运动 过 程 ; 飞行 蕉 对 井 锥 作用 的 嘛 
Kam, 飞行 器 运动 变 最 对 操 级 作用 的 传递 柄 数 以 及 相应 的 频率 特性 。 至 于 变 系 数 线 往 微 分 
方程 组 的 求解 万 法 具 作 简要 介绍 。 所 有 这 些 问题 都 属于 一 般 的 线性 系统 动力 学 问题 - 

党 二 种 情况 是 非 线 性 、 变 系数 的 微分 方程 组 。 第 一 、 三 章 及 第 五 章 《 部 分 ) Diere 
(BAH “28” JED) 即 属 于 这 种 情况 。 这 样 的 上 程 一 般 只 能 用 数值 法 求解 ， 册 在 给 定 的 
初 苔 条 什 和 给 定 的 摊 纵 输 大 作用 下 求 量 飞行 器 和 运动 的 数值 峙 果 。 为 此 要 介绍 微分 方程 组 的 数 
人 慎 求 解 方 法 。 

第 三 种 情况 足 非 线性 的 代数 或 超越 方程 组， 即 描述 飞行 器 平衡 状态 的 方 入 给。 这 样 的 广 
FALIE H RHE HEH A Ro 





SEN EI KATESH a RE A ti 


2-1 概 述 
ARE p A EEE RAER Eie EAA R E E S Du EAN TEA E 
如 下 一 般 形式 : 


dx 
EF =Ax+Bu {6-1} 


其 中 x jE n ERRARE: 
A Èn x n 维 系统 矩阵 (或 称 状态 录 阵 ); 
U Samu RS, 
BE nx u HPHHP. 
此 按照 承 数 冻结 法 ， 认 为 矩阵 上 和 8B 都 是 不 随时 间 而 变 的 什 阵 。 这 样 ， 以 方程 (6-1) 描 
述 的 动力 学 系统 侯 是 定常 线性 系统 。 这 样 的 条 统 可 以 用 解析 半 完 满 地 进行 研究 ， 并 生得 到 解 
HA FR 2-2 至 2-5 中 讲述 定常 线性 系统 的 动态 特性 分 析 ; 2-6 中 讲述 非 定常 (时空 ) SS 
性 系统 的 动态 特性 分 析 。 
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2-2 无 操纵 输 人 时 定常 线性 系统 的 动态 特性 


大 操纵 输 大 【 即 控制 矢量 4u=0) 时 作为 定常 线性 孙 统 的 飞行 器 的 运动 由 如 上 方程 来 措 
H: 


dx 
efi 


An (6-2) 
1. Jj Laplace 2282 RAR 
方程 (6-2) 进行 4aplace 变换 后 成 为 
sx(s)=Åx(s)tx(0) 
PxC) RA x(7) 的 Laplace 变换 式 ; x(0) 是 x(t) 的 初始 什 。 于 是 
Esl — AjxCs)=x(0) 


d ei Fal - AJ-ix(0) (6-3) 
RI ER EnH’, ESRR 
Gol. AT = ra rtl - A] (5-4) 


FRP ad ROR HERE, EX, nxn tEftpRa sE pE adEsl- AJ] KJ "re: 列 的 
RÆ nxn MEIER [31- AJ R JAT i ARR EA PR BIERE Tel- A] 
AEGA ARRETE n- 1 RIAR Im. meen 工 阶 多 项 式 ( 泊 7 王子 
Hie 2 Sm/A (8 Ts 时 ) 。 

A (6-4) 的 分 母 是 行列 式 。 命 

ADS jsl- A |= te, asl t estea (6-5) 
EE n META KOER A MESMA. 

MEKA (6-4), (6-5) RAX (6-3), 9 

adjlsl Al 

Aís) 
利川 Heaviside 公式 将 上 式 进行 Lanlacs 逆 变换 ， Mirah ath x) 引起 的 系统 过 e 
过 程 





x{s)= x(0) 








TS 一 Sadj[st 一 下 ] } sat E 
xi 一 之 上 二 (可 fon OE (6-6) 
其 中 s (k=l; 23 eo n) 是 特征 方程 式 
中 (ss 一 3 二 cr_isn lf E E A) (6-7) 


Dom. ARA RRR ERE A RE. F A EE BARRE AIE RETE -EE K i 
ERME. MER. A (6-6) 只 适用 下 特征 方程 式 (6-7) 为 互 蜡 的 二 个 特征 根 的 情况 。 
当 特 征 棚 存在 重 根 时 ， 需 参阅 有 关 文 献 〈 例 如 [1]) 。 
2. 特征 矢量 与 运动 模 态 
EEA HREH THEEDA BA AE sa 相对 应 的 特征 矢 Œ pao E 
EL BEIE ss 与 对 应 的 特征 矢量 pe 满足 如 下 方程 : 
页 pa 一 3xps (6-8) 
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或 Lal - Ain, =0 ( 
因为 特征 值 sa 已 满足 条 件 
IT Å= Afs) 50 
ai a ERRIRE E (6-9) BIF M 
JE (674) HÆ Esl- AJ, 4 
jsl- Aji s=[sl- Aiadifel A) 
当 seor bh Ltb, ed 
Csel Aidifszal AT L610) 
HEH, acjEsal- A] 的 任何 -: 列 均 渗 足 特 征 矢量 的 定 六 万 程 (6-9) DI. die: Al 
的 企 何 一 个 非 填 列 均 可 取 为 对 应 村 sa 的 特征 矢 En, GAEE EARI eD tT EAB 
ERA - 询 )。 这 就 是 计算 特征 矢量 的 实力 方法 。 
让 ”个 特征 天 晶 拨 列 组 成 的 矩阵 
P =[Pi Pepa] (6-11) 
称 为 畦 入 和 矩阵。 
如 果 取 
x (C1) =pae’ st 


则 有 有 


x= pase t" 
Ax= Apae'n’ 
考虑 到 条 件 (6-8) ， 可 见 x= Axo ZAA, HTE A k piet EARRA i REA 
特征 得 互 不 相等 《 助 无 重 特征 值 》 的 条 件 下 ,各 个 Pre 是 线性 无 关 的 。 所 以 包 合 ”个 任 
E tee 的 线性 组 合 


XC epee tësbse"tikrzekéefaeiel! (6-12) 
就 是 微分 方程 (6-2) 的 通 解 。 其 中 的 在意 常数 c 由 和 初 姑 条 件 决 定 
xfKO) 一 Clip 二 erps 十 …… 十 Cop， 《6--J3 


化 臣 标 莽 ， 则 上 上 式 成 为 4 个 线性 方程 ， 由 此 可 以 确定 # 个 常数 ci 。 
表达 式 6-12) 具有 重要 意义 ， 它 表明 ， 一 般 情 况 下 定常 线性 系统 的 运动 臣 田 王 个 形式 
为 pse'*' 的 基本 运动 组 成 的 。 每 -… 个 这 样 的 基本 运动 称 为 一 个 运动 模 态 。 
考虑 畦 征 值 的 两 种 情况 。 第 --， 当 某 个 8 为 实数 时 ，pke"z* 避 是 按 指数 规律 “生变 化 的 
非 周 期 模 态 ， 其 中 省 分 晶 保 持 着 癌 定 的 比例 关系 ; 第 二 ， 若 某 个 sx 是 复数 ， 则 必然 有 另 A 
Seri FZ 
Ska DNET 
E RER h E R 
Pr, ari SEER 
AEREI A Beti té de HE A 
Ca, #+1 At 
这 一 对 共 辊 复 模 态 之 和 为 
CPE t" +H Ceti Patie kri! 
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—(A+i H) tie Tint -有 Jr 一 05 
zc Zenit Ar = D ivoscc? - (44+ Br isinei (6-14) 
这 中 一 个 周期 模 态 〔 振 沪 模 态 )， 其 中 各 个 分 量 的 振幅 均 按 指数 规律 e wie. HIRR AEE 
的 比例 关系 ， 相 己 出 保持 固定 的 关系 。 
在 特殊 的 初始 策 件 x{0) 二 cop4 激 发 下 ， 采 统 的 运动 为 
Tt) 一 CpG 
KRHA, HE HUMR ERTE RE ， 系 统 的 运动 吉良 单个 模 态 组成 。 
2. PPIE IRE Ret N TE 
OG, HAREE PARER E x 直行 线性 变换 





x= Pz c6-15} 
IWSzk hm (6-2) WE 
Pz 一 APz 
成 z=P-:APz (6-16) 
D P-1AP 一 人 
Wl AP=PS 
Dy Ann, …p] 一 [pips…pv]S 
Zu EA (6-8) IEE ou SEHRE, (Sl, 290e n) RR luch 
Ši > 
s= â (6-17) 
s, d 
FR, SE (6-16) 成 为 
DEER (6-18) 
Bn Z= SiZ} (k=l. 2, =, n) (6-19) 


可 见 ， 经 过 线性 变换 (6-15) 后 形成 的 以 特征 坐标 矢 草 z RRIA A Ee n HTA 
(CHERDD 方程 。 显 然 ， 方 程 (6-19) 的 解 为 
za(1)—=e z(t0) (6-20) 
方程 组 【6-18) WRT EUH R 
setat 


g'az! | 


z(t)=' ` z0) 


=e"'z(0) (6-21) 
d. Esl-A' HHA 
ZEAE ARRE SMAA REREN mt GAES] Aimi) o a 


-1 adjlsl- Aj 
Esl- Ae TA] 
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Raast tR, as” 二 Ris 二 Rn 





To s poas li ie sie, ( 
首先 取 
KR. 一 0 ， cn 一】 
然后 按 如 下 递 拱 公 式 计算 ; 
RK. = orptil HÅR, gt 
1 
Cr- 一 E tr[ AR., Al ) { 


(k=l; 23 We HR 
Zä, 


2-3 定常 线性 系统 的 传递 负数 


6-229) 


6-23) 


HEBART AN ETA ETARE ESA E (6-1) o EM laplace TEA 


NET EIERE 
(el -Alx(s) = Bu(s) 
这 蛙 订 以 将 4 看 成 输入 ，x 表 成 输出 。 
输出 与 输入 的 Laplace 变换 式 之 比 称 为 系统 的 传递 醒 数 于 阵 


viel e — 
Gis) uts) =[sl -~1B 
=- aj[sl -AJB ( 
Ais) bi : 


落 x 为 上 维 伙 最 员 为 万 维 大 量 , 则 G(s) 为 ax 和 维 征 阵 ERII: s) EN a Ea H 


6-24) 





DS, ,对 


WAE, ERAR ERRAR E n RELIRA G) mAT z 0) MEE 


RERE adj[ sl - AJB ms, iban- 1 阶 多 项 式 。 


WELE (6-22) 、 (6-23) J 已 倪 绍 了 [sl - AJRIIANRAMEREH RER St, Dr 


BIRA (6-24) 计算 传递 画 数 矩阵 就 没有 什么 国难 了 。 


2-4 定常 线性 系统 对 揉 纵 输 入 的 响应 
1. EHAA Laplare 变换 


当 同 纵 输 大 ui 已 给 定时 ， 可 求 出 它 的 Lapac 变换 式 风 s) ， 然 后 根据 式 {6-24) $- 


viei-, Gisiuisi 
最 后 对 x(s) 作 Laplace 20 27 Ié. 有 即 可 求 出 飞行 器 的 及 应 过 程 x (1)。 
2. pkita he 


HAEA. iA du, A PARIR Bu G=). Amg, e)la He 
mD, BR MESE Aa KR, (1)。 这 时 有 
his (s)}= z: (s) = gy Cst; ales ee sl 
BELLI PpD žeh, O) TREA g: C) 的 Laplace 逆 变 换 式 
h(t)=L [gi;(s)] (8-25) 


s (ID: 


由 # mih, CHR Mn x mAH (z) 8: D DK Amar DE. Ié (6-25) JET, 就 得 
到 
H(t) =L'[GCs)] (6-26) 
2. Haha ra 
E -RAs (OWE TE, A 
Selz, Lëläcts) 
BI x, (5s) ApCs), (si 
DR, 可 应 利用 Laplace WHM ERER E E fi SE ni 


TOE hQ- rju, (rT)dr (6-27) 
ESTRATE K Jr 
xD 人 Èro- u(r)dr 
f-0f-i 
T ER SZ 


vuai-L HO- ur dr {6—28} 


2-5 定常 线性 系统 响应 的 其 它 问题 
除了 下 述 在 给 定 操 纵情 况 下 求 飞行 器 响应 时 间 过 程 的 问题 处， 还 有 两 类 化 频率 域内 研究 
的 问题 ， 简 单 地 介绍 一 个 。 
第 一 是 飞行 器 的 频率 响应 蛙 性 问题 。 设 操 约 量 以 频率 ao REIR 
uli) = aucos + pu) 
在 便 稳 状态 下 《行路 啊 应 量 x 亦 以 同样 频率 w EIRE Ra 
dE aolo EN 
EE SE AECH 
Alo)= H 
Plo) = — Fn 
MEHE o MHK. ATAWA CRIER) > # GCs) 是 x 对 的 传递 出 数 ， 则 
Geh, ass Gelz A(w)e'? (6-29) 
Gw) AS RN H io 
第 一 是 在 随机 输 人 作 几 下 飞行 囊 响 应 问题 。 设 输 大 草 u) 是 随机 过 程 ， 其 全 部 特性 以 
频谱 画 数 中 (ol 描述 。 输 市 对 输 基 的 传递 画 数 为 GC). AARU vir" 由 是 随机 过 程 ， 其 
Zi Le ro (on HAt. m th 可 接 如 下 关系 式 计算 : 
d. toi ltetioil Bulow) (6-30) 
LAE e A E AE LEEST Aa A FA A fa ai fE T Arta o EEA Linh T A 
KRAE DPHE 


- ill» 





Zb 非 定常 线性 系统 的 动态 响应 
于 定常 《时 变 ) RERA- ESEN 


2 x= Å(}x+B(7)u vr 21) 
ize IHARA A Jr 


有 TREIMER So, eh, OT eo ga Oe MEAR MR A 
Cga A) E ega Cr) te teq Il 
EFE (6-32) 的 通 解 ， 
Jä ie (6-31) 的 解 为 
TEDER QC YC) db 
Hif Q =g CE) gafr) e qa] 
yO S A E Holt) e OL 
将 式 【6-33) RA IR (6-31) ; 得 


OCOY QOY 一 ACTDQCDYCE BIODuCt) (6-34) 
HTE giel, D 11 ERRA 6-32) 
E QQ = ACFIQC) (6-35) 


因而 式 (6-34) 成 为 
Däi Bu) 


其 äi QC)B (u(r) (6-386) 
BTCA al Dë gem, KTTC) 存在 。 将 式 (5-36) 积分 ， 得 
y) “yG +f QB dr (6-37) 


gu (6-33) , # 
yita =Q etai) 
HERIR AA (6737) ， 然 后 代 久 式 (6-33); KI 
ECHTER Q-:(r)B(r)ugrjedr 
wT 
Di, r= QQ Cr) (6-38) 

则 最 后 荐 到 非 齐 次 方程 《6-31) 的 解 为 

JEU. BUSES I DU. t)Birjulr)dr (6-39) 
EREA Aa GO 引起 的 过 滤 过 程 ; 第 二 项 表示 出 aua) 8 引起 的 系 
Sdt? 


04 (6-38) 所 定义 的 Oo r) 称 为 状态 转移 矩阵 ， 式 中 的 OUO 为 基础 解 系 。 
TEM WEHE OC r) Afin TER: 


s 112 * 











dii, Dh 
D 加 2 加 | 
di, ri Së dt, ri Air, r) | (6-10) 
ils, t) if’ fr) 中 (r， £1) | 

®t p= (et dE ta) 


3 3 作为 非 线 性 系统 的 飞行 器 的 运动 方程 的 数值 来 解 

5-1 e $ 

TAS, T, ME (WA) iii Kiria A E EIRENE R R ia H Pia JER 
性 不 仅仅 枉 接 袁 现 在 运动 请 程 的 形式 上 ， 听 其 存在 于 飞行 器 的 气动 畦 性 和 发 动机 推力 特性 
中 ， 这 些 特 性 是 用 复杂 的 表达 式 或 实验 临 线 【着 转 化 成 表格 ) 来 描述 的 。 所 以 飞行器 运动 太 
释 组 在 大 多 数 情况 下 不 能 用 解析 法 求解 ， 只 能 用 数 介 法 求解 ， 即 在 具体 的 划 苔 条 件 和 具体 的 
拘 维 输入 条 件 下 求 由 飞行 器 运动 过 程 的 数值 解 。 

另 一 方面 ， 数 鼻 法 对 于 飞行 器 运动 方程 的 求解 几乎 是 力 能 的 ， 有 只 要 长 行 器 的 运动 能 用 方 
MRH. TERE FEAR EPRE RAER. ME KITARRE SA E H Bek h 
FIUK. HA ERER TARIA RE KARERE EARR TIA. 

RTE I HERRERA T RRF ARRE. EAA fi TI EHH 
Zap, EA Ti ATARA. BATRA A RRE (Pandi o AREFE 
JE, BRATIA iE T R A K aa D A E R E R ERE mA E o 

飞行 器 运动 方程 -: 般 地 可 以 归纳 成 如 下 形式 的 一 阶 常 微分 方程 组 ; 


i= filie nun Aan "ii: Xn) 


i= falio Xli Xar UG Xa) (6-41) 
a= falto TI? Aan Ta Xa) 7 
或 者 与 成 矢量 形式 
x=f(t, x) . (6-42) 
Jh X 一 [xi xad RER 


FEAH (6-41) 或 《6-42) 中 不 仅 真 接 反映 运动 方程 中 的 微分 方程 ， 市 并 在 面 数 户 ( 
Xia ans xa) PRR Ta A R p AE A A A EA APPR RA o 

微分 方程 组 (6-41) 或 (6-42) 的 数值 求解 方法 有 许多 种 。 基 本 上 可 以 分 为 多 步 法 和 一 
步 法 两 人 类 。 


5-2 ei 


多 步 法 的 特点 症 : 步 长 必须 为 常数 。 根 据 % Corinne, xats 
1’ Anne un 求 出 xx 在 后 一 步 的 数值 x 4+;。 多 步 法 双 可 以 进而 划分 为 外 插 法 和 预报 校正 
法 。 现 仅 忒 四 步 法 为 例 讲述 这 类 方法 的 实质 。 

DH Satir, EA Ri, RIRC pe Xans Enn un BI SI, Te 


-113 e 





Ins Fo Hee {laago fn-3); (tess faz), [taio fro fipo t) 作 
ETC, Stop DC EE OHE) ern, Sari) Mo DIO OC)MiE 
fst! 的 积分 第 作为 x 从 ia Ein WEE OL 6-1 Ai) : 


"kil 

XI "EA =| i F ldi 
D 
D 


按 上 水 患 路 经 过 推导 后 得 到 的 最 后 过 算 公 式 为 


55 59 37 9 
Xati E + 所 Ge Cedies fa “sa faa) (6-A4} 


HI, MEEMI- APREA, Ah 次 计算 导数 值 。 


T T 


Ta Tt É kej ik LEET t fgur ik- tk Fiat t 


SL 内 搬 法 


图 6-1 四 步 法 示意 图 

国 东 预报 核 上 正法 。 基 本 思路 是 : 首先 应 用 四 步 包 播 灶 的 公式 (6-44) ED tari aifi 
Rho AER, EI En, SES BE C aago Sa-a) Diane Jac) 
Ge fab’ Gano KE EENEG HSO GJA a B) t 的 积分 作为 
zc Mr, Sr, SE (HH 6-1 看 部 ) H HREH Xati? EAEE: 


f hatl 
se fd 
Wi 


由 十 在 这 全 区 间 内 周 的 是 内 插 哇 ， 收 它 比 预报 值 xir 更 为 精确 。 
所 以 预报 梳 正 法 的 实质 是 ; 先 用 录 播 法 预报 ， 然 后 用 内 插 法 校 止 。 计 算 公 式 册 纳 如 下 : 
站 S feat 9 WEE HEN ` 
firi =f Cata, xk+t1) ， (5-45) 
19 . 
Aki Tial d EEN 37 fe > f 11 yy faa) 
可 见 ， 攻 此 法 每 前 进 个 步兵 内， 需 啊 次 计算 导数 值 。 

总 之 ， 多 步 法 的 能 点 是 : L 由 于 在 任何 实际 问题 中 ， 已 知 初 始 条 件 上 只 有 II 
xo) ， 故 用 多 步 法 不 能 起 动 ， 必 须 借 助 于 一 步 法 【现下 述 ) 计算 出 5 个 点 以 后 才能 采用 多 步 
法 。 这 就 早教 计 算 程 序 复杂 化 : 2. 多 步 法 中 步 长 必须 为 常数 ， 因 而 不 能 灵活 地 改变 步 长 : 
3. 多 步 法 是 以 导数 1 的 连续 性 为 前 提 的 ， 因 此 当 肖 到 不 连续 点 (例如 艇 粳 阶 跃 仿 续 ， 发 动机 

es ]1l4，* 











突然 停车 的 情况 ) 时 此 方法 失效 ， 这 时 必须 重新 “起 动 ”。 册 于 以上 佣 点 ， 多 步 法 在 改行 器 
运动 方 稳 求解 中 使 用 不 多 。 


3-3 一 步 法 
一 步 法 的 畦 点 是 ， 内 要 有 了 一 点 的 数 秆 x 就 能 求 出 下 一 点 的 数值 x sr) .在 棚 于 一 步 法 的 
许多 方法 中 ， 这 里 仅 讲 述 标 维 四 阶 Runge~Kutta 靶 。 虽 然 它 并 不 是 最 先进 、 最 优越 的 方法 ， 
但 仍 是 虽 前 广泛 应 用 的 户 法 。 
PRIE T Runge- Kurta 中 的 计算 公式 如 下 : 
FaFa xa) 


Í= (ity h, xa HT 


D3 | = 
EN 
KE 


f= I(t +h, vat 2 WA (6-46) 
fo=f Uath, xat fch) 


H 
Stu AT e 由 (7 二 2870 :2fct+fp) 





公式 (6-46) 的 开 何 解 笠 是 ( 见 
4 6-2) : RA A C md d 
Ea dt 2. E, D ERER, 
JEH iR AE D IEA mmer HL 
Ip ABE Axo mixxat 
Ako 

一 步 法 的 优点 是 ， 本寺 能 
26 ah, Kiat rE ip, Sa 
REECH Dap, tut 
Jl RIRE ai, 图 6-2 Runge-Kutta 法 示意 转 





§ 4 飞行 器 平衡 状态 的 计算 
4-1 ES Ë 


改行 首 的 于 衡 状态 是 由 一些 非 线 性 方程 描述 的 。 几 于 许多 复杂 的 非 线性 关系 的 存在 ;这 
此 方 各 征 往 不 龙 川 星 搂 太 法 求解 ， 员 能 利用 尖 近 方法 求解 。 
这 些 非 线性 六 程 的 一 般 形 式 是 





f(x1s Aan e, Xa) =D ` 
Lal Xas Xas o Xa) =O | (EAF) 
Fatxo anren: Xa) =0 7 
KIK JER 
f(x}=0 (6-48) 


lls: 


在 众多 方法 由， 仅 讲述 此 种 ，Newron-Raphson 法 和 简单 选 代 法 。 
4-2 Newton-Raphson 法 


RECREATE rabu O, gf T -E e i e Ex ?为 此 将 方程 
(6-48) TEx 处 展开 成 Taylor 级 数 ， 仅 保留 线性 项 


f(x Cr J+ Zi dg vii viri j=9 
GEN 


yrth? —xir? -| 


Tina Ce? 、 4 
x (x 1 faxt?) (6-49) 
Ep THOR Ax Jacobi 短 阵 
afa fı AU. 
Axı dxa dir, j 


af Əfa fə Ais 
ax 一 SÉ 2x; PES (D bn) 


ai. fa Ai, 
OX] GER dr, 
HW AAEH s RPA (6-49) ， 得 到 x O, ni, "989 
SEMESTER 
l p kie At, ée, ERITA i $ ax (ci, B ovo m jais 
2o n) ， 丰 飞行 力学 问题 中 这 往往 是 困难 的 。 





4-5 简单 选 代 法 
沼 单 过 代 法 要 求 把 方 种 fx) 一 0 转化 戈 
x 三 中 Cx) (8-51) 
WR. CRATE r 次 近似 值 x “后 ， 即 可 求 出 个 -- 交 近似 位 
Xir+ = glx?) (5-52) 


MEEA. MEAR (6-52 0 RAA O x, es ABR 
KE RAA ko 
公式 (6-52) 称 为 按 全 部 坐标 透 代 公 式 。 另 一 个 更 有 效 也 更 简单 的 方式 称 为 接 部 分 坐标 
UHARIRI diodei 格式 ， 它 的 具体 公式 是 
wt sell, Ki, eg xf) 
KED px tO, Ki, wii? | EN 
AE E, x , wir? 
它 的 特点 是 : 尽 贡 利用 最 新 结果 ， 因 而 收 说 性 于 好 ; 而 且 计 算 程序 也 更 简单 ， 
Bigi MEO = 0 到 x 二 甲 (x) 的 转化 不 是 唯一 的 。 而 简单 选 代 法 的 收 敏 性 取决 了 这 
个 转化 :在 Siedel ERK dE 与 ,的 排列 磊 认 有关。 至 丁 如何 选 择 转 化 形式 ， 这 目 要 依 
第 经 验 积 累及 对 问题 的 物理 本 质 的 了 解 ; 还 可 以 通过 试验 来 选择 。 一 般 原则 是 ， 从 短 个 方程 
e JIG 。 


{uR H E EH fi A RAE o 
H AATED. EALA SERR ago jr 


1 2 
Menst zz yp- -5 PU S Cso tag +a ail: meint =0 

1 2 - 
Pail ryp): CS PUS eyot ega teže) — mgo- =0 15-54) 


i 
一 Pept a of baam, ot ano mez 6.) 0 


RK 中 KÉ bas EH Fp Ceos Ala Oga Cyto cko er, Mizos mgo mir, Pa RÉI 8 33 Z 已 
知 ， 答 求 平衡 状态 的 推力 只 ， 池 和 角 a 和 升降 舵 偏 角 56,。 

根据 经 验 和 对 物理 意义 的 了 解 ， 选择 的 转化 形式 基 : gro) Eër 
aH aM f 


"H 
L 2 pv’ S 


3-8 CORDAR) RAERD 





P =[ oS les Hara? toi mgin | EE 


出 第 三 个 方程 (DEAR) 求 扬 需 的 升降 舵 偏 角 
ô, -| x Pes "man" welt 
o Or Sha 


EDERLER, Keep P =0, L =0, Æ 用 Siedel RRR. Së 
{RTR Hb ii g EI E TAn F E PE o 

反之 ， 假 如 由 切 向 力 方 程 求 5 o BERJA EÆRP, MJETE a, WTA EER 
过 程 将 和 失败。 





附录 四 元 数 在 飞行 器 运动 方程 中 的 应 用 
i. SI 
在 本 书 中 {了 世间 迄今 大 多 数 飞 行 力 学 交 献 中 -一样 ) ， 我 们 利 用 Fuler fk e ERA E Er 


FRKA. Phin. 机体 华 标 对 相对 于 地 面 华 标 妹 的 方位 是 用 人 妃 航 角 罗 、 储 仰角 # m 
角 ? 来 描述 的 。 而 了 pnler 角 变 化 率 与 机 体 旋 转角 速度 的 关系 是 用 如 下 的 微分 方 称 求 建立 


中 


的 : 
D 1 -voy sin cost siny sing cot, ye 
| a l = | Oo cosy eosi 一 siny :cos L ,| ri) 
f 0 sinye cosy FER 


这 就 是 第 二 章 方 程 (2-46) o RAMNES R E T EITE Wiete cot 0p 
特殊 情况 下 in ERE KELN E ALa E tT P TREH RARE o Eii A rh gth 
现 奇 罕 点 ， 使 计算 开 法 进行 下 去 ， 也 使 模拟 试验 被 追 中国 。 避 如 这 种 奇异 情况 的 有 效 办 法 就 
是 存 飞 行 器 运动 方程 中 用 所 谓 “ 四 元 数 ” (guaternion) KRE RAI Euler 角 。 


2. 四 元 数 的 定义 和 性 质 
kuk TT 
Q =g t giit gaj tysk iZ} 
Sat, jo k ERRER WETER FIRE SH M: 
TER ET CEA 
ij=-ji, jk= -kj; ki= -ik | {3) 
ij =k, jk=i, ki=j 
由 定义 式 (2) 和 规定 的 运算 法 则 式 (3) 可 见 ， 对 于 四 元 数 而 绢 ， 琵 法 的 分 记 律 仍然 成 








下 
Q+ QQ tQ a 
{B si iR ETE r. 
Qr, 
Di ER O HERRE A 
B90- gii- gaj- gak EN 
不 难 导 出 


QB= DO af tai tgi tgi = NIQ) 
Hop v QORA EE EE 
如 果 NQS MORAER CURET EREN itto 
"e 118 e 


2. 用 四 元 数 表 示 坐 标 变换 


在 第 一 章 已 指出 ， 检 何 啊 个 坐标 和 《例如 5。 与 9) 之 闻 的 关系 可 以 "tie: H 
RS, IS ARA (Sa) SA Ami ZE — A Euler Di, AE ei db HE 
(Sa) Eo 

ER TEIEGA E Irina, 3 PITIRAA A SA >. ERË 
豆 找 蔓 个 旋转 畏 和 天 和 -个 角度 a ， 使 得 坐标 系 S i OR pA E a aep A 
"Wë, (nei, y Zik Br Bs OELH TEREA R ET 
"R a mgh OR MAREA 





D. cospe (i =i, 2, 3) (5) 
HN "rn 
a a 
yn co > 9= Pain g (i=l, 2, 3) (ei 
求 构 成 四 元 数 
Q$ -gat gii+gaj+ gk t7) 


A T R AE Sr bk Se, MS MRR. WES ELARA Sa Sa H Hi 
WRS o WOARE AR KiE ERR o 

这 里 订 能 会 提出 问题 : HARK Eue AFATE AEE E k iT ai: 
Bo mme Tt: 原因 在于 四 个 参数 之 间 还 在 在 一 个 学 化 条 件 





gs tgi rgd+g3=1 (8) 
或 彰 说 ， 出 于 旋转 织品 凡 的 -个 方向 余 蔬 之 间 存 在 约束 条 件 
Pitpatpi=l (8°) 


EER UTT Aiia o 
H ERIE séC) Aën DERO, Oxa, Oy, Otr, Oz, O2 A9 KA MSR D 
学 的 方法 表示 在 网 上 中。 
在 这 些 球 所 一 角形 中 ， 已 知 的 边 为 
XaYa S Yata Zaa T Xp > Yna = Xp T 2 
Rx = Rax, = Py ane eonp 
Ry, = PY, = P a= arc corpa 
Rz, = R2; = Êa Sare cos ps 


2 知 的 角 为 
ago 
= fZ Xaj 5 TIR 
Aas, xy 5 LY NYa = Lz R2 a 
m R aA 2 


AERKEAn Elfen CA Zu? 
HpXas Hen Ha 
Zh Za. Zeto 


在 球面 三 骨 形 x ,xs rh, HAI R WMR E A x Rach JO Ao Genf 
. 19> 








KEE 





vs Sa xa) = cost vi sin D pecs 


ER 
ERMA Rx nu. un, AAAG An, D acte, ls, EIS vr 
D-—vosff rof, tsin f sinf acolar ója) (10) 
由 此 可 解 出 68 


ERMS ARY hs WA (i, Po WRA G) 为 已 知 ， 故 可 求 上 出 第 三 边 
cosf xaya) = cos p cosh a tsin ein fi seosd]s 
SIT, Stu EI ka, DI enn, Rëm EF te ne Haas 
Dë, AE, Sal, Zs2s) 也 可 以 用 奖 伺 的 方法 求 出 。 而 这 万 个 边 《实际 上 是 岂 个 角 ) HR 
SECH. Sax, 2 bp E be 2 Së REN CS 
ER COSC Xia co xala) vëeizexë, ue Xa 
, ! = EZ ost HaHa) orza) IP | 
z 


EI Cosi za, oosfzng Ä Gosizz, 
GERS Sg 角形 大 经 过 许多 代数 运算 《推导 从 上 略 ) DE, Sea 
H Dër -gå Ztogs goga) 2(9i9s 十 dogs) 


MÉ 
Héi 2(cids 十 gogs) ` oi -dt+tos -a Boss ` gegi) 1, | (12) 
Z; Zones Goga} 2{gags+ goga? oë qt gtg "2 


Ek LG IR 8 T EEEa. Eit Ti FURRA: 
xai typj t zeK= (ga t gii tgaj + oak) (Cx, ity j+ 
+k) Cgo 一 ol -qaj ` gak) 


{ [1} 


Si 


或 ET Hdëfek (13) 
如 果 说 由 ,到 S ;的 转移 是 5S, 线 轴 0QR 转 过 角 a 而 实现 的 ， 则 反 过 来 由 3 到 5. 的 转移 
yS RAO RA -am KE. EAO BA (7) 相对 应 的 是 


QES 911 ~ gaj -gk (14) 
HALENA A an TER: 
v= ği (15) 
该 公式 (13) 亦 可 写成 
总 一 名 or (16) 
如 果 接 连 进行 两 次 变换 
rs—= Qs HEET det 
则 re=QiQ r EÜE 
Es re= Qir OO? 
而 V =Q? (i7) 
LORATA 
go = ge oi? gip qi?) -gf of o Zii oi? 
` ent At gig tgi oft?) | (an 
oa 一 gti gi 一 qi qi P Igi gin g gor?) ` 
DEE AAA -qig iS tgi gi tgi or 


Hi gét, HAN + e a AQ AIQ, W Ro 


d. 四 元 数 与 也 uler 角 的 关系 
按 第 AMURE AAE pR Sa 首先 绕 轴 ys 转 过 偏 航 和 前兆 , 然后 绕 轴 zz ead H Ar 
d. Er, Siss, BIRERE S Ee. TR, ESA ni Zem 
Qı? {cos(yi2), D, -sin(yi2), 0} 
Q= {tos{8/2); D, 0, -sin(#/2)} 
Qa = {eo0s(v/2), -sinCyi2), 0, 0} 
上 是 击 地 面 系 到 机 体系 的 坐标 变换 四 元 数 为 , 


QSO Wa 
展 于 后 得 到 怠 的 各 元 素 为 
go en pi) eos 2) eos(p/2) — siet /2)sin(d/2)sin( /2)} 
= coy hein A Diet /2) - sin(p/2)cos(d/2)cos(p/2) | (19) 


ga- cosl p/2)cost tio) sin( /2) -sin{ y; 2)sin(d/2)cos{Ņ/ 2) | 
ds = 一 rogfT 2 lais dl 2 )eost Y2) int ie /2) cot 2)sin(y 2) 
AR huler 角 求 四 元 数 分 量 的 公式 。 
另 一 方面 ， 车 把 等 式 (12) 中 的 xs， Hi, zu H Yin Zen 把 >x。、 Van >, 换 成 xy 
nu, za, SE SE Dm Desk TA 50, 的 坐标 变换 矩阵 Bi, mer Bi 已 在 第 二 
e 12} e 





章 的 式 〔2-2) hënt, TE- LEET LBE Eu e, DEE Sie Ei dun E 
数 束 Euler 角 的 公式 【以 组 会 形 式 给 出 更 为 方便 ) 
sind- Zou: - goga) 
codsiny — — St sgs 十 God 
onBrosy =q oi roigii ， (2i 
Zoe tgga) | 


也 :一 = ~ 7. 
s” gë tgi -gigs 


5 四 无 数 分 量 的 微分 方程 式 
d kiréé ag, KRAAD EE Ion KITAR IECH (aw 
其 地 说 时 出 增 面 邓 到 机 体系 的 变换 四 元 数 )》 的 分 量 也 不 断 地 变化 。 因 此 需要 用 微分 浪 程式 来 
28 AAR EE. 
H UAE RA 
[do j goney DBgo19g D909. ， 
$1| Zei Zeie Bq/ay D 
| da |7; Baay agay 20130 IS 
La J .Ogs/QF 99s713 Əqsa’ | 
Kirps, (äi, ue (2, RIDO Sé 《19) 推导 出 。 
"IS, amb (iTA Léna, REAT. Bana taR E h bar hm 


式 
#0 | d KÉ 
` | 8 8 | i s 
MAM > Tu o a) -| (21) 
| #2 | | gs “qe qi 
| i] | ] Ed as 
Wë 、 Dn "On 


Hea oe, vu. AE UE DEH DE. 


56. 利用 四 元 数 建立 的 飞行 器 运动 方程 组 
这 时 采用 的 力 方程 和 力矩 方程 都 是 在 机 体 坐 标 孙 中 建立 的 运 动 bie, ur 体系 
(大 看 第 二 章 ) 。 原 来 的 大 多 数 方程 都 可 以 保持 不 变 ， 仅 两 组 方程 需要 改变 。 第 一 是 利用 公 
K CO 把 包含 kle 角 的 力 方程 [ 见 第 二 吉方 程 (2-31)] 改 咸 包含 四 元 数 的 形式 ， 即 
el Lers Tou ai OUy )= Peosgpp — Qose 8 +Y sing - 
- Zveosasinp - 2mgigiga 一 Qaqa) 
gd H FaU T Upe ) = Psings + Osinacosp +F rosg + (22) 
|- Zsingsin f — mg oi — gz tgs q3) 


loas 
m e -二 (E ~ Qsinf + Zeos f - Zmgiosga Fäng) 


第 二 是 用 方程 《21) REREH (2-46) o 

用 于 航 迹 溪 标 系 中 的 力 方 称 如 何 改 造 ， 这 里 就 不 讲述 了 。 

Raati Fo A KIT ESA EAN, RAAE DICH koler 角 ， 和 有 了 两 个 优点 : 
第 -` ， 可 以 导 郭 出 现 奇 措 虚 ， 而 保证 在 任何 条 件 下 订 算 或 模拟 试验 不 会 因此 而 中 断 : 第 二 ， 
减少 了 - 些 三 角 图 数 的 计算 ， 困 而 可 以 节约 机 时 。 但 也 有 缺点 ， 四 元 数 的 物理 意义 不 如 
Balier 角 清 乱 ， 面 有 由 阿 元 数 与 飞行 仪表 的 读数 波 有 起 楼 联系 ， 而 Evler 角 则 可 以 在 仪表 中 前 示 
BEE 


e, J33 e, 


